N«RSAR 


NORSAR  Scientific  Report  No.  2-2002 

Semiannual  Technical  Summary 

1  January  -  30  June  2002 


Frode  Ringdal  (ed.) 


Kjeller,  August  2002 


DISTRIBUTION  STATEMENT  A  I 

Approved  for  Public  Release  I 
Distribution  Unlimited 


20030312  243 


SFCURITY  CLASSIFICATION  OF  THIS  PAGE  - — — — — 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 

Not  applicable _ 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 


Fom  Approved 
0MB  No.  0704-01 88 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified  _ 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 

Not  Applicable  _ 


2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER<S) 
Scientific  Rep.  2-2002 


Approved  for  public  release:  distribution  unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUM8ER(S) 

Scientific  Rep.  2-2002 


6.  HAMiOPPiUroRMlWOTGANIZOToS - Ub.  OFHCE  SVMWl.  ■'(  7..  NAME  OF  MONIIOKINS  0MANI2AT10N 

(If  appliceoie}  I 


NORSAR 


6c  ADDRESS  (Oty,  State,  and  ZIP  Code) 
Post  Box  53 

NO-2027  Kjeller,  Norway 


HQ/AFTAC/TTS 


7b.  ADDRESS  (Cfty,  State,  and  ZIP  Code) 
Patrick  AFB,  FL  32925-6001 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable)  ^ 

DTRA/NTPO  Contract  No.  F08650-01-C-0055 


Sc  ADDRESS  (Cfty,  State,  and  ZIP  Code) 


10.  SOURCE  OF  FUNDING  NUMBERS _ 

«  ipROJEa  I  TASK 

no^orsar 

Phase  3  Task  5.0 


I  WORK  UNIT 


■  ■  PE7mi^r  I  ■  lASK  |ves./r>F> 

1515  Wilson  Blvd.,  Suite  720 
Arlington,  VA  22209  _ 


11.  TITLE  (Include  Security  aassification) 

Semiannual  Technical  Summary,  1  January  -  30  June  2002  (Unclassified) 


12.  PERSONAL  AUTHOR(S) 

13a.  TYPE  OF  REPORT - |13b.  TIME  COVEReF  |14.  DATE  OF  REPORT  (Year,  Month,  Day)  115.  PAGE  COUNf 

Scientific  Summary  from  1  Jan  02  To30Jun02j^  2002  Aug  1 


16.  SUPPLEMENTARY  NOTATION 


18.  SUBJECT  TERMS  (Continue  on  reverse  If  necessary  and  identify  by  block  number) 

NORSAR,  Norwegian  Seismic  Array 


19.  ABSTRACT  {Continue  on  reverse  if  necessary  and  Identify  by  block  number) 


17.  COSATI  CODES _ 


FIELD  I  GROUP  I  SUB-GROUP 


8 


This  report  describes  the  research  activities  carried  out  at  NORSAR  under  Contract  No.  F08650-01  -C-0055  for  the  perio 
1  January  -30  June  2002.  In  addition,  it  provides  summary  information  on  operation  and  maintenance  (O&M)  activities  at 
the  Norwegian  National  Data  Center  (NDC)  during  the  same  period.  Research  activities  described  in  this  r^ort,  as  well 
as  transmission  of  selected  data  to  the  United  States  NDC,  are  funded  by  the  United  States  Department  of  Defense  The 
O&M  activities,  including  operation  of  transmission  links  within  Norway  and  to  Vienna,  Austria,  are  being  funded  jointly  by 
the  CTBTO/PTS  and  the  Norwegian  Government,  with  the  understanding  that  the  funding  of  IMS-related  activities  will 
gradually  be  transferred  to  the  CTBTO/PTS.  The  O&M  statistics  presented  in  this  report  are  included  for  the  purpose  of 
completeness,  and  in  order  to  maintain  consistency  with  earlier  reporting  practice. 

(cont.) 


20.  DISTRIBUTION /AVAIUBILITY  OF  ABSTRACT 
DuNCLASSIFIEDAJNLIMITEO  □  SAME  AS  RPT. 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Lt.  Col.  William  S.  Jones 


rtn  JUN  8S 


□  OTIC  USERS 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


22 b.  TELEPHONE  f/nctode  Area  Code)  22c.  OFFICE  SYMBOL 

(407)  494-7985  AFTAC/TTS 


NORSAR  Sci.  Rep.  2-2002 


August  2002 


Abstract  (cont.) 

The  seismic  arrays  operated  by  the  Norwegian  NDC  comprise  the  Norwegian  Seismic  Array 
(NOA),  the  Norwegian  Regional  Seismic  Array  (NORES),  the  Arctic  Regional  Seismic  Array 
(ARCES)  and  the  Spitsbergen  Regional  Array  (SPITS).  This  report  also  presents  statistics  for 
additional  seismic  stations  which  through  cooperative  agreements  with  institutions  in  the  host 
countries  provide  continuous  data  to  the  NORSAR  Data  Processing  Center  (NDPC).  These  sta¬ 
tions  comprise  the  Finnish  Regional  Seismic  Array  (FINES),  the  Hagfors  array  in  Sweden  and 
the  regional  seismic  array  in  Apatity,  Russia. 

The  NOA  Detection  Processing  system  has  been  operated  throughout  the  period  with  an  average 
uptime  of  100.00%.  A  total  of  1872  seismic  events  have  been  reported  in  the  NOA  monthly  seis¬ 
mic  bulletin  from  January  through  June  2002.  On-line  detection  processing  and  data  recording  at 
the  NDC  of  NORES,  ARCES  and  FINES  data  have  been  conducted  throughout  the  period.  Data 
from  two  small-aperture  arrays  at  sites  in  Spitsbergen  and  Apatity,  Kola  Peninsula,  as  well  as  the 
Hagfors  array  in  Sweden,  have  also  been  recorded  and  processed.  Processing  statistics  for  the 
arrays  for  the  reporting  period  are  given. 

A  summary  of  the  activities  related  to  the  GSETT-3  experiment  and  experience  gained  at  the 
Norwegian  NDC  during  the  reporting  period  is  provided  in  Section  4.  Norway  is  now  contribut¬ 
ing  primary  station  data  from  two  seismic  arrays:  ARCES  and  NOA  and  one  auxiliary  array 
(SPITS).  These  data  are  being  provided  to  the  IDC  via  the  global  communications  infrastructure 
(GCI).  Continuous  data  from  all  three  arrays  are  in  addition  being  transmitted  to  the  US  NDC. 
The  performance  of  the  data  transmission  to  the  US  NDC  has  been  satisfactory  during  the 
reporting  period. 

The  PrepCom  has  encouraged  states  that  operate  IMS-designated  stations  to  continue  to  do  so 
on  a  voluntary  basis  and  in  the  framework  of  the  GSETT-3  experiment  until  the  stations  have 
been  certified  for  formal  inclusion  in  IMS.  So  far  among  the  Norwegian  stations,  the  NOA  and 
the  ARCES  array  (PS27  and  PS28  respectively)  have  been  certified.  We  envisage  continuing  the 
provision  of  data  from  these  and  other  Norwegian  IMS-designated  stations  in  accordance  with 
current  procedures. 

Summaries  of  six  scientific  and  technical  contributions  are  presented  in  Chapter  6  of  this  report. 

Sectioii  6.1  contains  a  report  from  the  meeting  of  the  IDC  Technical  Experts  Group  on  Seismic 
Event  Location  in  Oslo,  Norway  on  22-26  April  2002.  This  was  the  fourth  annual  meeting  of  the 
Experts  Group  in  support  of  Working  Group  B  of  the  CTBTO  Preparatory  Commission.  The 
meeting  was  held  jointly  with  the  IDC  Technical  Experts  Group  on  Seismoacoustic  Event 
Screening.  The  purpose  of  the  meeting  was  to  support  the  ongoing  calibration  and  screening 
efforts  of  the  IDC  and  in  particular  to  review  progress  toward  developing  regionalized  travel 
times  to  improve  the  quality  of  location  estimates  of  seismic  events  reported  in  the  IDC  bulle¬ 
tins. 

The  meeting  was  attended  by  sixty-six  technical  experts,  coming  from  twelve  signatory  coun¬ 
tries  and  the  Provisional  Technical  Secretariat.  Dr.  Frode  Ringdal  of  Norway  chaired  the  meet- 
ing,  which  was  organized  into  four  sessions,  and  included  Working  Group  discussions  to  address 
the  technical  issues  in  detail.  Recommendations  from  the  meeting  have  been  presented  to  Work¬ 
ing  Group  B  in  Vienna  during  its  June  2002  session  (CTBTAVGB/TL-2/70). 
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Section  6.2  is  entitled  “Research  in  regional  seismic  monitoring”  and  contains  a  summary  paper 
presented  at  the  24th  Annual  Seismic  Research  Review,  describing  a  continued  effort  in  regional 
monitoring  research.  We  have  used  data  from  the  regional  networks  operated  by  NORSAR  and 
the  Kola  Regional  Seismological  Centre  (KRSC)  to  assess  the  seismicity  and  characteristics  of 
regional  phases  of  the  European  Arctic.  An  especially  interesting  recent  seismic  event,  which  has 
been  analyzed  in  detail,  occurred  on  Novaya  Zemlya  on  23  February  2002.  This  small  event 
(magnitude  about  3.0)  is  the  first  seismic  event  detected  in  or  near  Novaya  Zemlya  since  the  two 
Kara  Sea  events  of  16  August  1997.  The  event  was  recorded  with  particularly  high  SNR  by  the 
International  Monitoring  System  (IMS)  arrays  Spitsbergen  and  ARCES,  and  by  the  Amderma  sta¬ 
tion  (operated  by  KRSC)  south  of  Novaya  Zemlya.  We  have  located  the  event  at  the  eastern  coast 
of  Novaya  Zemlya,  approximately  100  km  NE  of  the  former  nuclear  test  site.  Our  attempts  to  esti¬ 
mate  the  depth  of  this  event  proved  inconclusive,  and  a  depth  of  0  could  not  be  ruled  out. 

We  have  begun  an  effort  to  develop  an  optimized  site-specific  threshold  monitoring  system  for  the 
Lop  Nor  test  site  in  China.  Initial  results  show  that,  among  the  available  stations,  the  IMS  Makan- 
chi  array  in  Kazakhstan  (at  an  epicentral  distance  of  6  degrees)  is  the  most  sensitive  station  for 
detecting  events  at  Lop  Nor,  followed  by  some  of  the  best  teleseismic  arrays,  for  example  the 
NORES  array  in  Norway.  Combining  the  data  by  using  the  threshold  monitoring  technique,  we 
have  found  that  the  Lop  Nor  site  can  be  monitored  down  to  a  threshold  close  to  3.2.  This 
threshold  is  about  0.8  magnitude  units  higher  than  the  corresponding  capability  for  the  Novaya 
Zemlya  test  site,  for  which  the  highly  sensitive  arrays  ARCES  and  Spitsbergen  are  both  at  a  dis¬ 
tance  of  only  10  degrees.  We  note  that  no  stations  inside  China  were  selected  for  the  Lop  Nor 
study. 

Section  6.3  is  a  design  study  of  a  possible  reconfiguration  of  the  Spitsbergen  array  (SPITS).  This 
array  needs  a  technical  refurbishment  before  it  can  be  certified  for  inclusion  in  the  IMS,  and  the 
motivation  for  the  present  study  is  to  determine  possible  options  for  reconfiguring  the  array  so  as 
to  improve  its  monitoring  capabihties.  The  SPITS  array  is  situated  in  a  geological/tectonic  envi¬ 
ronment  that  is  very  different  from  that  of  Fennoscandia.  This  results  in  a  number  of  characteristic 
features  that  must  be  taken  into  account  when  processing  SPITS  data. 

The  SPITS  array  often  reports  a  large  number  of  detections,  which  can  reach  thousands  per  day.  A 
detailed  analysis  shows  that  the  majority  of  these  detections  are  real  seismic  signals  caused  hy 
small  sources  located  at  close  distances.  The  numerous  local  events  typically  show  P  onsets,  no 
well  defined  S  onsets,  and  dominant  Rg  onsets.  Investigating  the  mean  backazimuths  of  these  Rg 
phases,  one  can  see  a  strong  correlation  with  the  direction  to  a  nearby  coal-mine  area  and  to 
nearby  glaciers.  To  avoid  spatial  abasing  of  these  signals  a  denser  inner  ring  of  seismometers 
would  be  helpful.  To  improve  the  separation  between  phase  types  based  on  apparent  velocities,  an 
increase  in  the  array  aperture  is  required.  This  would  also  provide  improved  resolution  in  the  azi¬ 
muth  estimates. 

At  the  SPITS  array  we  observe  a  relatively  large  number  of  high  SNR  P  detections  from  regional 
events  without  any  corresponding  detection  of  an  S  phase.  The  regional  S-phases  at  SPITS  in  gen¬ 
eral  have  the  largest  ampbtudes  and  SNR  on  the  horizontal  components.  The  current  SPITS  con¬ 
figuration  has  only  one  three-component  instrument  at  site  SPB4.  With  more  three-component 
sites,  the  detection  capabibty  of  S-phases  can  be  improved  by  creating  S-velocity  beams  from  the 
rotated  radial  and  transverse  components,  as  now  done  for  ARCES  and  NORES.  In  addition,  the 
F-K  analysis  can  be  performed  on  the  horizontal  components. 
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The  paper  suggests  specific  extensions  for  the  SPITS  array  to  increase  its  detection  and  monitor¬ 
ing  capabilities.  These  points  will  form  the  basis  for  our  proposed  Spitsbergen  upgrade/refur¬ 
bishment. 

Section  6.4  is  a  study  of  seismic  data  from  a  series  of  underground  cavity  decoupled  chemical 
explosions  in  the  Alvdalen  region  of  Sweden,  with  the  aim  of  quantitatively  predicting  the 
decoupling  effect  of  underground  cavities  of  known  dimensions  and  shape.  A  database  of  seis¬ 
mic  waveform  data  has  been  compiled  for  seven  such  cavity  explosions.  Only  three  factors  dif¬ 
ferentiate  these  events;  the  quantity  of  explosive,  the  composition  and  arrangement  of  the 
explosives,  and  the  size  and  configuration  of  the  explosion  chamber.  In  addition,  the  Swedish 
Armed  Forces  regularly  detonate  outdated  ammunition  at  ground  level  at  a  site  approximately 
15  km  from  the  underground  chambers.  Seismic  waveforms  from  12  such  events  were  added  to 
the  database  to  provide  a  comparison  with  the  decoupled  explosions. 

In  order  to  quantify  the  differences  between  these  events,  we  examine  the  power  density  spectra 
of  the  signals  recorded  by  NOA,  NORES  and  HFS.  The  recordings  from  explosions  in  the  cavi¬ 
ties  at  Alvdalen  are  characterized  by  a  large  amount  of  high-frequency  energy  (8-16  Hz)  com¬ 
pared  with  the  energy  at  lower  frequencies  (2-5  Hz).  Such  spectral  characteristics  are  not 
observed  for  the  surface  explosions.  The  effect  is  more  pronounced  for  P  than  for  S  waves. 

With  regard  to  decoupling  efficiency,  three  10  ton  explosions  in  a  lOOOm^  chamber  appear  to 
have  been  more  efficiently  decoupled  than  three  5  ton  explosions  in  chambers  of  200m^  and 
300m  .  The  signals  from  the  explosions  in  the  smaller  chambers  were  far  stronger,  in  spite  of 
significantly  lower  yields,  and  displayed  a  less  pronounced  high-to-low  frequency  energy  ratio 
than  the  explosions  in  the  larger  chambers. 

We  plan  to  assemble  NOA,  HFS  and  NORES  recordings  of  earthquakes  located  at  distances 
similar  to  Alvdalen  and  Mossibranden,  but  at  different  azimuths,  in  order  to  make  comparative 
studies  of  their  spectral  characteristics.  Comparisons  between  records  from  decoupled  explo¬ 
sions  and  ripple  fired  mining  explosions  are  also  recognized  as  an  important  future  target. 

Section  6.5  is  a  study  of  experimental  Site-Specific  Generalized  Beamforming  (SSGBF)  applied 
to  the  Lop  Nor  test  site  in  China.  The  SSGBF  approach  supplements  the  SSTM  system 
described  in  the  preceding  NORSAR  Semiannual  Technical  Summary.  We  also  comment  briefly 
on  the  relative  merits  of  the  two  approaches  (SSTM  and  SSGBF),  and  the  extent  to  which  they 
complement  each  other  in  this  particular  case  study. 

In  the  present  study  we  have  calibrated  all  available  high-quality  seismic  stations  within  regional 
distance  from  Lop  Nor  and  supplemented  them  with  the  best  IMS  arrays  at  teleseismic  dis¬ 
tances.  As  a  calibration  data  base,  we  use  several  past  nuclear  explosions  at  the  test  site,  com¬ 
bined  with  nearby,  well-recorded  earthquakes.  We  present  some  preliminary  results  in  applying 
site-specific  generalized  beamforming  to  the  Lop  Nor  test  site,  using  two  data  sets  for  perfor¬ 
mance  testing.  The  first  data  set  covers  one  day  00  September  2001),  and  in  these  data  the 
recordings  of  4  explosions  were  scaled  and  embedded  at  some  (but  not  all  of)  the  stations.  Sec¬ 
ondly,  a  10  day  test  period  with  data  from  August  2  through  1 1,  2001  has  been  used.  Examples 
of  the  data  analysis  for  these  data  sets  are  presented. 

The  paper  concludes  that  the  Site-Specific  Generalized  Beamforming  performance  tests  for  the 
Lop  Nor  test  site  were  successful.  The  combination  of  the  SSTM  and  the  SSGBF  techniques 
appears  to  be  a  veiy  promising  new  development.  Each  of  these  two  approaches  has  its  specific 
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strengths.  The  SSTM  technique  has  as  its  main  strength  the  abihty  to  display  the  real  seismic 
field,  without  regard  to  “station  detector  performance”.  On  the  other  hand,  the  SSGBF  tech¬ 
nique  takes  advantage  of  the  individual  station  detector  outputs,  and  uses  this  combined  infor¬ 
mation  to  narrow  down  the  number  of  possible  candidates  for  events  in  the  target  area. 

Section  6.(5  is  a  study  of  seismic  threshold  monitoring  apphed  to  the  Barents  Sea  region.  A 
database  comprising  a  total  of  45  events,  selected  to  provide  the  best  possible  ray  path  coverage 
of  the  Barents  Sea  and  adjacent  areas,  was  compiled  and  reanalyzed  in  a  consistent  manner. 
This  resulted  in  new  regional  attenuation  relations  for  Pn  and  Sn,  together  with  a  preferred 
average  velocity  model  to  he  used  for  predicting  the  travel  times  of  regional  phases.  We  have 
now  applied  these  attenuation  relations  to  investigate  a  regional  threshold  monitoring  scheme 
for  the  Barents  Sea  area. 

A  grid  system  with  an  approximately  100-km  grid  spacing  was  deployed  for  the  Barents  Sea 
region,  and  the  observations  at  the  arrays  ARCES,  SPITS,  FINES  and  NORES  were  then  used 
for  calculating  threshold  magnitudes  for  each  of  the  grid  points.  During  an  interval  without 
seismic  signals,  the  threshold  magnitudes  showed  large  variations  over  the  region,  and,  in  par¬ 
ticular,  in  the  vicinity  of  each  array.  However,  for  the  region  around  the  island  of  Novaya 
Zemlya  the  variations  are  modest,  varying  around  a  mean  of  magnitude  2. 1-2.2. 

We  have  studied  the  tradeoff  between  beam  area  coverage  (“beam  sharpness”)  and  monitoring 
performance  for  the  site-specific  threshold  monitoring  of  the  Novaya  Zemya  test  site.  For  this 
purpose,  we  have  deployed  a  very  dense  grid  system  covering  an  area  of  500  x  500  km  around 
the  test  site,  and  calculated  magnitude  thresholds  for  each  grid  point  for  a  two  hour  time  inter¬ 
val,  which  included  a  seismic  event  in  the  region  (23  February  2002, 01.12.21,  magnitude  3). 
Using  the  arrays  ARCES,  SPITS,  FINES  and  NORES,  we  find  that  even  for  a  target  region 
with  radius  as  large  as  100  km,  the  variations  in  threshold  magnitudes  are  all  within  0.2  magni¬ 
tude  units.  This  applies  both  for  the  time  interval  including  the  event  and  for  background  noise 
conditions.  For  the  investigated  station  geometry,  it  will  therefore  be  meaningful  to  represent 
the  monitoring  threshold  of  the  entire  Novaya  Zemlya  region  using  the  values  of  a  single  target 
point,  together  with  the  d priori  determined  uncertainty  bounds. 

For  areas  with  larger  variations  in  threshold  magnitudes,  like  in  the  vicinity  of  the  arrays,  a 
l(X)-km  radius  target  region  will  obviously  show  larger  differences  between  the  maximum  and 
minimum  values.  Examples  illustrating  this  point  are  shown. 
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1  Summary 

This  report  describes  the  research  activities  carried  out  at  NORSAR  under  Contract  No.  F08650- 
Ol-C-0055  for  the  period  1  January  -  30  June  2002.  In  addition,  it  provides  summary  informa¬ 
tion  on  operation  and  maintenance  (O&M)  activities  at  the  Norwegian  National  Data  Center 
(NDC)  during  the  same  period.  Research  activities  described  in  this  report,  as  well  as  transmis¬ 
sion  of  selected  data  to  the  United  States  NDC,  are  funded  by  the  United  States  Department  of 
Defense.  The  O&M  activities,  including  operation  of  transmission  hnks  within  Norway  and  to 
Vienna,  Austria  are  being  funded  jointly  by  the  CTBTO/PTS  and  the  Norwegian  Government, 
with  the  understanding  that  the  funding  of  all  IMS-related  activities  will  gradually  be  transferred 
to  the  CTBTO/PTS.  The  O&M  statistics  presented  in  this  report  are  included  for  the  purpose  of 
completeness,  and  in  order  to  maintain  consistency  with  earlier  reporting  practice. 

The  seismic  arrays  operated  by  the  Norwegian  NDC  comprise  the  Norwegian  Seismic  Array 
(NOA),  the  Norwegian  Regional  Seismic  Array  (NORESX  the  Arctic  Regional  Seismic  Array 
(ARCES)  and  the  Spitsbergen  Regional  Array  (SPITS).  This  report  also  presents  statistics  for 
additional  seismic  stations  which  through  cooperative  agreements  with  institutions  in  the  host 
countries  provide  continuous  data  to  the  NORSAR  Data  Processing  Center  (NDPC).  These  sta¬ 
tions  comprise  the  Finnish  Regional  Seismic  Array  (FINES),  the  Hagfors  array  in  Sweden  and 
the  regional  seismic  array  in  Apatity,  Russia. 

The  NOA  Detection  Processing  system  has  been  operated  throughout  the  period  with  an  average 
uptime  of  100.00%.  A  total  of  1872  seismic  events  have  been  reported  in  the  NOA  monthly  seis¬ 
mic  bulletin  from  January  through  June  2002.  On-line  detection  processing  and  data  recording  at 
the  NDC  of  NORES,  ARCES  and  FINES  data  have  been  conducted  throughout  the  period.  Data 
from  two  small-aperture  arrays  at  sites  in  Spitsbergen  and  Apatity,  Kola  Peninsula,  as  well  as  the 
Hagfors  array  in  Sweden,  have  also  been  recorded  and  processed.  Processing  statistics  for  the 
arrays  for  the  reporting  period  are  given. 

A  summary  of  the  activities  related  to  the  GSETT-3  experiment  and  experience  gained  at  the 
Norwegian  NDC  during  the  reporting  period  is  provided  in  Section  4.  Norway  is  now  contribut¬ 
ing  primary  station  data  from  two  seismic  arrays:  ARCES  and  NOA  and  one  auxiliary  array 
(SPITS).  These  data  are  being  provided  to  the  IDC  via  the  global  communications  infrastructure 
(GCI).  Continuous  data  from  all  three  arrays  are  in  addition  being  transmitted  to  the  US  NDC. 
The  performance  of  the  data  transmission  to  the  US  NDC  has  been  satisfactory  during  the 
reporting  period. 

The  PrepCom  has  encouraged  states  that  operate  IMS-designated  stations  to  continue  to  do  so 
on  a  voluntary  basis  and  in  the  framework  of  the  GSETT-3  experiment  until  the  stations  have 
been  certified  for  formal  inclusion  in  IMS.  So  far  among  the  Norwegian  stations,  the  NOA  and 
the  ARCES  array  (PS27  and  PS28  respectively)  have  been  certified.  We  envisage  continuing  the 
provision  of  data  from  these  and  other  Norwegian  IMS-designated  stations  in  accordance  with 
current  procedures. 

Summaries  of  six  scientific  and  technical  contributions  are  presented  in  Chapter  6  of  this  report. 

Section  6.1  contains  a  report  from  the  meeting  of  the  IDC  Technical  Experts  Group  on  Seismic 
Event  Location  in  Oslo,  Norway  on  22-26  April  2002.  This  was  the  fourth  annual  meeting  of  the 
Experts  Group  in  support  of  Working  Group  B  of  the  CTBTO  Preparatory  Commission.  The 
meeting  was  held  jointly  with  the  IDC  Technical  Experts  Group  on  Seismoacoustic  Event 
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Screening.  The  purpose  of  the  meeting  was  to  support  the  ongoing  calibration  and  screening 
efforts  of  the  IDC  and  in  particular  to  review  progress  toward  developing  regionalized  travel  times 
to  improve  the  quality  of  location  estimates  of  seismic  events  reported  in  the  IDC  bulletins. 

The  meeting  was  attended  by  sixty-six  technical  experts,  coming  from  twelve  signatory  countries 
and  the  Provisional  Technical  Secretariat.  Dr.  Frode  Ringdal  of  Norway  chaired  the  meeting, 
which  was  organized  into  four  sessions,  and  included  Working  Group  discussions  to  address  the 
technical  issues  in  detail.  Recommendations  from  the  meeting  have  been  presented  to  Working 
Group  B  in  Vienna  during  its  June  2002  session  (CTBTAVGB/TL-2/70). 

Section  6.2  is  entitled  “Research  in  regional  seismic  monitoring”  and  contains  a  summary  paper 
presented  at  the  24th  Annual  Seismic  Research  Review,  describing  a  continued  effort  in  regional 
monitoring  research.  We  have  used  data  from  the  regional  networks  operated  by  NORSAR  and 
the  Kola  Regional  Seismological  Centre  (KRSC)  to  assess  the  seismicity  and  characteristics  of 
regional  phases  of  the  European  Arctic.  An  especially  interesting  recent  seismic  event,  which  has 
been  analyzed  in  detail,  occurred  on  Novaya  Zemlya  on  23  February  2002.  This  small  event 
(magnitude  about  3.0)  is  the  first  seismic  event  detected  in  or  near  Novaya  Zemlya  since  the  two 
Kara  Sea  events  of  16  August  1997.  The  event  was  recorded  with  particularly  high  SNR  by  the 
Litemational  Monitoring  System  (IMS)  arrays  Spitsbergen  and  ARCES,  and  by  the  Amderma  sta¬ 
tion  (operated  by  KRSC)  south  of  Novaya  Zemlya.  We  have  located  the  event  at  the  eastern  coast 
of  Novaya  Zemlya,  approximately  100  km  NE  of  the  former  nuclear  test  site.  Our  attempts  to  esti¬ 
mate  the  depth  of  this  event  proved  inconclusive,  and  a  depth  of  0  could  not  be  ruled  out. 

We  have  begun  an  effort  to  develop  an  optimized  site-specific  threshold  monitoring  system  for  the 
Lop  Nor  test  site  in  China.  Initial  results  show  that,  among  the  available  stations,  the  IMS  Makan- 
chi  array  in  Kazakhstan  (at  an  epicentral  distance  of  6  degrees)  is  the  most  sensitive  station  for 
detecting  events  at  Lop  Nor,  followed  by  some  of  the  best  teleseismic  arrays,  for  example  the 
NORES  array  in  Norway.  Combining  the  data  by  using  the  threshold  monitoring  technique,  we 
have  found  that  the  Lop  Nor  site  can  be  monitored  down  to  a  threshold  close  to  m(j  3.2.  This 
threshold  is  about  0.8  magnitude  units  higher  than  the  corresponding  capability  for  the  Novaya 
Zemlya  test  site,  for  which  the  highly  sensitive  arrays  ARCES  and  Spitsbergen  are  both  at  a  dis¬ 
tance  of  only  10  degrees.  We  note  that  no  stations  inside  China  were  selected  for  the  Lop  Nor 
study. 

Section  6.3  is  a  design  study  of  a  possible  reconfiguration  of  the  Spitsbergen  array  (SPITS).  This 
array  needs  a  technical  refurbishment  before  it  can  be  certified  for  inclusion  in  the  IMS,  and  the 
motivation  for  the  present  study  is  to  determine  possible  options  for  reconfiguring  the  array  so  as 
to  improve  its  monitoring  capabilities.  The  SPITS  array  is  situated  in  a  geological/tectonic  envi¬ 
ronment  that  is  very  different  from  that  of  Fennoscandia.  This  results  in  a  number  of  characteristic 
features  that  must  be  taken  into  account  when  processing  SPITS  data. 

The  SPITS  array  often  reports  a  large  number  of  detections,  which  can  reach  thousands  per  day.  A 
detailed  analysis  shows  that  the  majority  of  these  detections  are  real  seismic  signals  caused  by 
small  sources  located  at  close  distances.  The  numerous  local  events  typically  show  P  onsets,  no 
well  defined  S  onsets,  and  dominant  Rg  onsets.  Investigating  the  mean  backazimuths  of  these  Rg 
phases,  one  can  see  a  strong  correlation  with  the  direction  to  a  nearby  coal-mine  area  and  to 
nearby  glaciers.  To  avoid  spatial  aliasing  of  these  signals  a  denser  inner  ring  of  seismometers 
would  be  helpful.  To  improve  the  separation  between  phase  types  based  on  apparent  velocities,  an 


2 


NORSAR  Sci.  Rep.  2-2002 


August  2002 


increase  in  the  array  aperture  is  required.  This  would  also  provide  improved  resolution  in  the 
azimuth  estimates. 

At  the  SPITS  array  we  observe  a  relatively  large  number  of  high  SNR  P  detections  from  regional 
events  without  any  corresponding  detection  of  an  S  phase.  The  regional  S-phases  at  SPITS  in 
general  have  the  largest  amplitudes  and  SNR  on  the  horizontal  components.  The  current  SPITS 
configuration  has  only  one  three-component  instrument  at  site  SPB4.  With  more  three-compo¬ 
nent  sites,  the  detection  capability  of  S-phases  can  be  improved  by  creating  S-velocity  beams 
from  the  rotated  radial  and  transverse  components,  as  now  done  for  ARCES  and  NORES.  In 
addition,  the  F-K  analysis  can  be  performed  on  the  horizontal  components. 

The  paper  suggests  specific  extensions  for  the  SPITS  array  to  increase  its  detection  and  monitor¬ 
ing  capabilities.  These  points  will  form  the  basis  for  our  proposed  Spitsbergen  upgrade/refur¬ 
bishment. 

Section  6.4  is  a  study  of  seismic  data  from  a  series  of  underground  cavity  decoupled  chemical 
explosions  in  the  Alvdalen  region  of  Sweden,  with  the  aim  of  quantitatively  predieting  the 
decouphng  effect  of  underground  cavities  of  known  dimensions  and  shape.  A  database  of  seis¬ 
mic  waveform  data  has  been  compiled  for  seven  such  cavity  explosions.  Only  three  factors  dif¬ 
ferentiate  these  events;  the  quantity  of  explosive,  the  composition  and  arrangement  of  the 
explosives,  and  the  size  and  configuration  of  the  explosion  chamber.  In  addition,  the  Swedish 
Armed  Forces  regularly  detonate  outdated  ammunition  at  ground  level  at  a  site  approximately 
15  km  from  the  underground  chambers.  Seismic  waveforms  from  12  such  events  were  added  to 
the  database  to  provide  a  comparison  with  the  decoupled  explosions. 

In  order  to  quantify  the  differences  between  these  events,  we  examine  the  power  density  spectra 
of  the  signals  recorded  by  NOA,  NORES  and  HFS.  The  recordings  from  explosions  in  the  cavi¬ 
ties  at  Alvdalen  are  characterized  by  a  large  amount  of  high-frequency  energy  (8-16  Hz)  com¬ 
pared  with  the  energy  at  lower  frequencies  (2-5  Hz).  Such  spectral  characteristics  are  not 
observed  for  the  surface  explosions.  The  effect  is  more  pronounced  for  P  than  for  S  waves. 

With  regard  to  decoupling  efficiency,  three  10  ton  explosions  in  a  lOOOm^  chamber  appear  to 
have  been  more  efficiently  decoupled  than  three  5  ton  explosions  in  chambers  of  200m^  and 
300m^.  The  signals  from  the  explosions  in  the  smaller  chambers  were  far  stronger,  in  spite  of 
significantly  lower  yields,  and  displayed  a  less  pronounced  high-to-low  frequency  energy  ratio 
than  the  explosions  in  the  larger  chambers. 

We  plan  to  assemble  NOA,  HFS  and  NORES  recordings  of  earthquakes  located  at  distances 
similar  to  Alvdalen  and  Mossibranden,  but  at  different  azimuths,  in  order  to  make  comparative 
studies  of  their  spectral  characteristics.  Comparisons  between  records  from  decoupled  explo¬ 
sions  and  ripple  fired  mining  explosions  are  also  recognized  as  an  important  future  target. 

Section  6.5  is  a  study  of  experimental  Site-Specific  Generalized  Beamforming  (SSGBF)  apphed 
to  the  Lop  Nor  test  site  in  China.  The  SSGBF  approach  supplements  the  SSTM  system 
described  in  the  preceding  NORSAR  Semiannual  Technical  Summary.  We  also  comment  briefly 
on  the  relative  merits  of  the  two  approaches  (SSTM  and  SSGBF),  and  the  extent  to  whieh  they 
complement  each  other  in  this  particular  case  study. 

In  the  present  study  we  have  eahbrated  all  available  high-quality  seismic  stations  within  regional 
distance  from  Lop  Nor  and  supplemented  them  with  the  best  IMS  arrays  at  teleseismic  dis¬ 
tances.  As  a  cahbration  data  base,  we  use  several  past  nuclear  explosions  at  the  test  site,  com- 
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bined  with  nearby,  well-recorded  earthquakes.  We  present  some  preliminary  results  in  applying 
site-specific  generalized  beamforming  to  the  Lop  Nor  test  site,  using  two  data  sets  for  perfor¬ 
mance  testing.  The  first  data  set  covers  one  day  (10  September  2001),  and  in  these  data  the 
recordings  of  4  explosions  were  scaled  and  embedded  at  some  (but  not  all  of)  the  stations.  Sec¬ 
ondly,  a  10  day  test  period  with  data  from  August  2  through  1 1,  2001  has  been  used.  Examples 
of  the  data  analysis  for  these  data  sets  are  presented. 

The  paper  concludes  that  the  Site-Specific  Generalized  Beamforming  performance  tests  for  the 
Lop  Nor  test  site  were  successful.  The  combination  of  the  SSTM  and  the  SSGBF  techniques 
appears  to  be  a  very  promising  new  development.  Each  of  these  two  approaches  has  its  specific 
strengths.  The  SSTM  technique  has  as  its  main  strength  the  ability  to  display  the  real  seismic 
field,  without  regard  to  “station  detector  performance”.  On  the  other  hand,  the  SSGBF  tech¬ 
nique  takes  advantage  of  the  individual  station  detector  outputs,  and  uses  this  combined  infor¬ 
mation  to  narrow  down  the  number  of  possible  candidates  for  events  in  the  target  area. 

Section  6.6  is  a  study  of  seismic  threshold  monitoring  applied  to  the  Barents  Sea  region.  A 
database  comprising  a  total  of  45  events,  selected  to  provide  the  best  possible  ray  path  coverage 
of  Ae  Barents  Sea  and  adjacent  areas,  was  compiled  and  reanalyzed  in  a  consistent  manner. 
This  resulted  in  new  regional  attenuation  relations  for  Pn  and  Sn,  together  with  a  preferred 
average  velocity  model  to  be  used  for  predicting  the  travel  times  of  regional  phases.  We  have 
now  applied  these  attenuation  relations  to  investigate  a  regional  threshold  monitoring  scheme 
for  the  Barents  Sea  area. 

A  grid  system  with  an  approximately  100-km  grid  spacing  was  deployed  for  the  Barents  Sea 
region,  and  the  observations  at  the  arrays  ARCES,  SPITS,  FEMES  and  NORES  were  then  used 
for  calculating  threshold  magnitudes  for  each  of  the  grid  points.  During  an  interval  without 
seismic  signals,  the  threshold  magnitudes  showed  large  variations  over  the  region,  and,  in  par¬ 
ticular,  in  die  vicinity  of  each  array.  However,  for  the  region  around  the  island  of  Novaya  Zem- 
lya  the  variations  are  modest,  varying  around  a  mean  of  magnitude  2. 1-2.2. 

We  have  studied  the  tradeoff  between  beam  area  coverage  (“beam  sharpness”)  and  monitoring 
performance  for  the  site-specific  threshold  monitoring  of  the  Novaya  Zemya  test  site.  For  this 
purpose,  we  have  deployed  a  very  dense  grid  system  covering  an  area  of  500  x  500  km  around 
the  test  site,  and  calculated  magnitude  thresholds  for  each  grid  point  for  a  two  hour  time  inter¬ 
val,  which  included  a  seismic  event  in  the  region  (23  February  2002, 01.12.21,  magnitude  3) 
Using  the  arrays  ARCES,  SPITS,  FINES  and  NORES,  we  find  that  even  for  a  target  region 
with  radius  as  large  as  100  km,  the  variations  in  threshold  magnitudes  are  all  within  0.2  magni¬ 
tude  units.  This  apphes  both  for  the  time  interval  including  the  event  and  for  background  noise 
conditions.  For  the  investigated  station  geometry,  it  will  therefore  be  meaningful  to  represent 
the  monitoring  threshold  of  the  entire  Novaya  Zemlya  region  using  the  values  of  a  single  target 
point,  together  with  the  d priori  determined  uncertainty  bounds. 

For  areas  with  larger  variations  in  threshold  magnitudes,  like  in  the  vicinity  of  the  arrays,  a 
100-km  radius  target  region  will  obviously  show  larger  differences  between  the  maximum  and 
minimum  values.  Examples  illustrating  this  point  are  shown. 


Erode  Ringdal 
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2  Operation  of  International  Monitoring  System  (IMS)  Stations 
in  Norway 

2.1  PS27  —  Primary  Seismic  Station  NOA 

The  average  recording  time  was  100%,  the  same  as  for  the  previous  reporting  period. 

Monthly  uptimes  for  the  NORSAR  on-hne  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  fine,  data  center  operation)  affecting  this  task  were  as 
follows: 

January  :  100% 

February  :  100% 

March  :  100% 

April  :  100% 

May  :  100% 

June  :  100% 


Fig.  2.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORSAR  data  in  our  tape  archive,  on  a  day-hy-day  basis  for  the  reporting  period. 

J.  Torstveit 
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Fig.  2.1.1.  The  figure  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availability  of 
NOA  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period.  (Page  1  of  2, 
Jan-Mar  2002). 
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NOA  Uptime  for  May.  2002  (This  month:  100.00  percent) 
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NOA  uptime  for  Jun.  2002  (This  month:  100.00  percent) 


Fig.  2.1.1.  (cont.)  (Page  2  of  2,  Apr-Jun  2002). 
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NOA  Event  Detection  Operation 

In  Table  2.1.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  hsts  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detector,  the 
total  number  of  detections  processed  by  the  automatic  event  processor  (EPX)  and  the  total' 
number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases  and  total). 


Total 

DPX 

Total 

EPX 

Accepted  Events 

Sum 

Daily 

P-phases 

Core 

Phases 

Jan 

12,660 

870 

184 

67 

251 

8.1 

Feb 

11,032 

721 

162 

48 

210 

^B9 

Mar 

12,508 

918 

266 

82 

348 

11.2 

Apr 

7,919 

853 

295 

67 

362 

12.1 

May 

5,144 

699 

286 

53 

339 

10.9 

Jun 

6,352 

830 

293 

69 

362 

12.1 

55,615 

4,891 

1,486 

386 

1,872 

10.3 

Table  2.1.1,  Detection  and  Event  Processor  statistics,  1  January  -  30  June  2002. 


NOA  detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  001,  2002, 
through  day  181, 2002,  was  55,615,  giving  an  average  of  307  detections  per  processed  day  (181 
days  processed). 

B.  Paulsen 
U.  Baadshaug 
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2.2  PS28  —  Primary  Seismic  Station  ARCES 

The  average  recording  time  was  99.19%  as  compared  to  99.20%  for  the  previous  period. 
Table  2.2.1  lists  the  reasons  for  and  time  periods  of  the  main  downtimes  in  the  reporting  period. 


Date 

Time 

Cause 

02  Feb 

0949 

- 

Power  failure 

03  Feb 

-  1800 

24  Apr 

1232 

-  1301 

Maintenance  in  array 

Table  2.2.1.  The  main  interruptions  in  recording  of  ARCES  data  at  NDPC,  1  January  - 
30  June  2002. 


Monthly  uptimes  for  the  ARCES  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmission  lines,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows; 


January 

99.99% 

February 

95.21% 

March 

100% 

April 

99.93% 

May 

100% 

June 

100% 

J.  Torstveit 
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Fig.  2.2.1.  The  figure  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availability 
of  ARCES  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period.  (Page  1 
of  2,  Jan-Mar  2002) 
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Event  Detection  Operation 
ARCES  detections 

The  number  of  detections  (phases)  reported  during  day  001, 2002,  through  day  181,  2002,  was 
109,509,  giving  an  average  of  605  detections  per  processed  day  (181  days  processed). 

Events  automatically  located  by  ARCES 

During  days  001,  2002,  through  181,  2002,  7,986  local  and  regional  events  were  located  by 
ARCES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of 
44.1  events  per  processed  day  (181  days  processed).  62%  of  these  events  are  within  300  km, 
and  85%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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2.3  AS72  —  Auxiliary  Seismic  Station  Spitsbergen 

The  average  recording  time  was  97.60%  as  compared  to  99.99%  for  the  previous  reporting 
period. 

Table  2.2.1  hsts  the  reasons  for  and  time  periods  of  the  main  downtimes  in  the  reporting  period. 


Date 

Time 

Cause 

20  Jan 

1801 

- 

Windmill  broken 

22  Jan 

-  1226 

Batteries  charged  by  petrol 
aggregate 

28  Jan 

2113 

- 

Empty  batteries 

29  Jan 

1335 

Batteries  charged  by  petrol 
aggregate 

06  Feb 

1155 

-  1840 

Empty  batteries  (charged  again) 

24  Feb 

0116 

- 

Empty  batteries 

25  Feb 

1304 

Batteries  charged  by  petrol 
aggregate 

18  Mar 

New  windmill  installed 

Table  2.2.1.  The  main  interruptions  in  recording  of  Spitsbergen  data  atNDPC,  1  January 
-  30  June  2002. 


Monthly  uptimes  for  the  Spitsbergen  on-hne  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  fine,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


January 

92.07% 

February 

93.63% 

March 

99.98% 

April 

99.94% 

May 

100% 

June 

99.98% 

J.  Torstveit 
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Fig.  23.1.  The  figure  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availability  of 
Spitsbergen  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period.  (Pase  1 
of  2,  Jan-Mar  2002). 
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Event  Detection  Operation 

Spitsbergen  array  detections 

The  number  of  detections  (phases)  reported  from  day  001,  2002,  through  day  181,  2002,  was 
295,419,  giving  an  average  of  1641  detections  per  processed  day  (180  days  processed).  ' 

Events  automatically  located  by  the  Spitsbergen  array 

During  days  001,  2002,  through  181,  2002,  36,467  local  and  regional  events  were  located  by 
the  Spitsbergen  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an 
average  of  202.6  events  per  processed  day  (180  days  processed).  68%  of  these  events  are 
within  300  km,  and  84%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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2.4  AS73  —  Auxiliary  Seismic  Station  at  Jan  Mayen 

The  IMS  auxiliary  seismic  network  will  include  a  three-component  station  on  the  Norwegian 
island  of  Jan  Mayen.  The  station  location  given  in  the  protocol  to  the  Comprehensive  Nuclear- 
Test-Ban  Treaty  is  70.9°N,  8.7°W. 

The  University  of  Bergen  has  operated  a  seismic  station  at  this  location  since  1970.  An  invest¬ 
ment  in  the  new  station  at  Jan  Mayen  wiU  be  made  in  due  course  and  in  accordance  with  Prep- 
Com  program  and  budget  decisions.  A  so-called  Parent  Network  Station  Assessment  for  AS73 
was  completed  in  April  2002.  Work  is  now  underway  to  prepare  for  the  installation  of  a  vault  at 
a  new  location  (71.0°N,  8.5°W)  recently  approved  by  the  PrepCom.  In  the  meantime,  data  from 
the  existing  seismic  station  on  Jan  Mayen  are  being  transmitted  to  the  NDC  at  Kjeller  and  to 
the  University  of  Bergen  via  a  VSAT  link  installed  in  April  2000. 

J.  Fyen 

2.5  IS37  —  Infrasound  Station  at  Karasjok 

The  IMS  infrasound  network  wiU  include  a  station  at  Karasjok  in  northern  Norway.  The  coor¬ 
dinates  given  for  this  station  are  69.5°N,  25.5°E.  These  coordinates  coincide  with  those  of  the 
primary  seismic  station  PS28. 

A  site  survey  for  this  station  was  carried  out  during  June/July  1998  as  a  cooperative  effort 
between  the  Provisional  Technical  Secretariat  of  the  CTBTO  and  NORSAR.  Analysis  of  the 
data  collected  at  several  potential  locations  for  this  station  in  and  around  Karasjok  has  been 
completed.  The  results  of  this  analysis  have  led  to  a  recommendation  on  the  exact  location  of 
the  inffasound  station.  This  location  needs  to  be  surveyed  in  detail.  The  next  step  wiU  be  to 
approach  the  local  authorities  to  obtain  the  permission  necessary  to  establish  the  station.  Sta¬ 
tion  installation  is  now  expected  to  take  place  in  the  year  2003. 

S.  Mykkeltveit 

2.6  RN49  —  Radionuclide  Station  on  Spitsbergen 

The  IMS  radionuclide  network  will  include  a  station  at  Longyearbyen  on  the  island  of  Spitsber¬ 
gen,  with  location  78.2°N,  16.4°E,  as  given  in  the  protocol  to  the  Comprehensive  Nuclear-Test- 
Ban  Treaty.  These  coordinates  coincide  with  those  of  the  auxihary  seismic  station  AS72. 
According  to  PrepCom  decision,  this  station  will  also  be  among  those  IMS  radionuclide  sta¬ 
tions  that  will  have  a  capability  of  monitoring  for  the  presence  of  relevant  noble  gases  upon 
entry  into  force  of  the  CTBT. 

A  site  survey  for  this  station  was  carried  out  in  August  of  1999  by  NORSAR,  in  cooperation 
with  the  Norwegian  Radiation  Protection  Authority.  The  site  survey  report  to  the  PTS  con¬ 
tained  a  recommendation  to  establish  this  station  at  Plataberget,  some  20  km  away  from  the 
Treaty  location.  The  PrepCom  approved  the  corresponding  coordinate  change  in  its  meeting  in 
May  2000.  The  station  installation  was  part  of  PrepCom’s  work  program  and  budget  for  the 
year  2000.  The  infrastructure  for  housing  the  station  equipment  has  been  established,  and  a 
noble  gas  detection  system,  based  on  the  Swedish  “SAUNA”  design,  was  instaUed  at  this  site  in 
May  2001,  as  part  of  PrepCom’s  noble  gas  experiment.  A  particulate  station  (“ARAME” 
design)  was  installed  at  the  same  location  in  September  2001.  Currently,  the  two  systems  are 
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undergoing  testing  and  evaluation.  It  is  expected  that  the  particulate  station  will  be  certified  by 
the  end  of  2002. 

S.  Mykkeltveit 
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3  Contributing  Regional  Seismic  Arrays 
3.1  NORES 

Average  recording  time  was  88.92%  as  compared  to  70.32%  for  the  previous  period. 

The  outage  that  starts  on  1 1  June  12.04  and  continues  throughout  the  month  was  due  to  a  thun¬ 
derstorm  that  caused  major  damage  to  the  equipment. 


Date 

Time 

Cause 

10  Mar 

1300  -  1630 

UPS  failure  at  HUB 

15  Mar 

0246  -  0950 

UPS  failure  at  HUB 

11  Jun 

1204  - 

HUB  damage  due  to  thunderstorm 

Table  3.1.1.  The  main  interruptions  in  recording  of  NORES  data  at  the  Norwe¬ 
gian  NDC,  1  January  -  30  June  2002. 


Monthly  uptimes  for  the  NORES  on-line  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  fine,  data  center  operation)  affecting  this  task  were  as 
follows: 


January 

100% 

February 

99.95% 

March 

98.58% 

April 

99.98% 

May 

99.99% 

June 

35.00% 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availabihty  of 
NORES  data  in  our  tape  archive  on  a  day-by-day  basis  for  the  reporting  period. 


J.  Torstveit 
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Fig.  3.1.1.  The  figure  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availability  of 
MORES  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1  of  2 
Jan-Mar  2002). 
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NORES  Event  Detection  Operation 
NORES  detections 

The  number  of  detections  (phases)  reported  from  day  001,  2002,  through  day  181,  2002,  was 
139,063,  giving  an  average  of  858  detections  per  processed  day  (162  days  processed). 

Events  automatically  located  by  NORES 

During  days  (X)l,  2002,  through  181,  2002,  3,646  local  and  regional  events  were  located  by 
NORES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of 
22.5  events  per  processed  day  (162  days  processed).  43%  of  these  events  are  within  300  km, 
and  76%  of  these  events  are  within  1000  km. 

U.  Baadshaug 


22 


NORSAR  Sci.  Rep.  2-2002 


August  2002 


3.2  Hagfors  (IMS  Station  ASlOl) 

The  average  recording  time  was  99.97%  as  compared  to  99.91%  for  the  previous  reporting 
period. 

Monthly  uptimes  for  the  Hagfors  on-hne  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as 
follows: 

January  :  99.84% 

February  :  100% 

March  :  100% 

April  :  99.99% 

May  :  100% 

June  :  100% 

J.  Torstveit 
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Fig.  3.2.1  (cont.)  (Page  2  of  2,  Apr-Jun  2002). 
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Hagfors  Event  Detection  Operation 

Hagfors  array  detections 

The  number  of  detections  (phases)  reported  from  day  001,  2002,  through  day  181,  2002,  was 
102,667,  giving  an  average  of  567  detections  per  processed  day  (181  days  processed). 

Events  automatically  located  by  the  Hagfors  array 

During  days  001, 2002,  through  181, 2002,  3,016  local  and  regional  events  were  located  by  the 
Hagfors  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  16.7  events  per  processed  day  (181  days  processed).  56%  of  these  events  are  within  300  km, 
and  83%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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3.3  FINES  (IMS  station  PS17) 

The  average  recording  time  was  100%  as  compared  to  99.99%  for  the  previous  reporting 
period. 

Monthly  uptimes  for  the  FINES  on-hne  data  recording  task,  taking  into  account  all  factors 
(field  installations,  transmissions  fine,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 


January 

100% 

February 

100% 

March 

100% 

April 

100% 

May 

100% 

June 

100% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
FINES  data  in  our  tape  archive  on  a  day-by-day  basis  for  the  reporting  period. 

J.  Torstveit 
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Fig.  33.1.  The  figure  shows  the  uptime  for  the  data  recording  task  or,  equivalently,  the  availability 
of  FINES  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1  of 
2,  Jan-Mar  2002 ). 
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FINES  Event  Detection  Operation 
FINES  detections 

The  number  of  detections  (phases)  reported  during  day  001,  2002,  through  day  181, 2002,  was 
44,060,  giving  an  average  of  243  detections  per  processed  day  (181  days  processed). 


Events  automatically  located  by  FINES 

During  days  001,  2002,  through  181,  2002,  2,220  local  and  regional  events  were  located  by 
FINES,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average  of  12.9 
events  per  processed  day  (181  days  processed).  77%  of  these  events  are  within  300  km,  and 
87%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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3.4  Apatity 

The  average  recording  time  was  99.94%  in  the  reporting  period  compared  to  98.31%  during  the 
previous  period. 

Monthly  uptimes  for  the  Apatity  on-line  data  recording  task,  taking  into  account  aU  factors 
(field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were  as  fol¬ 
lows: 

January  :  100% 

February  :  99.98% 

March  :  100% 

April  :  99.64% 

May  :  100% 

June  :  100% 

J.  Torstveit 
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Fig.  3.4.1.  The  figure  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability 
of  Apatity  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period  (Page  1  of 
2,  Jan-Mar  2002 }. 
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APA  uptime  for  Apr.  2002  (This  month:  99. 64  percent) 


APA  uptime  for  May.  2002  (This  month.  100.00  percent) 


1  2  3  4  5  6  7  8  9  10  1 1  12  13 


APA  Uptime  for  Jun.  2002  (This  month:  100.00  percent) 

Fig.  3.4.1  (com.)  (Page  2  of  2,  Apr-Jun  2002). 
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Apatity  Event  Detection  Operation 


Apatity  array  detections 

The  number  of  detections  (phases)  reported  from  day  001,  2002,  through  day  181,  2002,  was 
1 64,442,  giving  an  average  of  904  detections  per  processed  day  ( 1 8 1  days  processed). 

As  described  in  earUer  reports,  the  data  from  the  Apatity  array  is  transferred  by  one-way  (sim¬ 
plex)  radio  links  to  Apatity  city.  The  transmission  suffers  from  radio  disturbances  that  occa¬ 
sionally  result  in  a  large  number  of  small  data  gaps  and  spikes  in  the  data.  In  order  for  the 
communication  protocol  to  correct  such  errors  by  requesting  retransmission  of  data,  a  two-way 
radio  link  would  be  needed  (duplex  radio).  However,  it  should  be  noted  that  noise  from  culmral 
activities  and  from  the  nearby  lakes  cause  most  of  the  unwanted  detections.  These  unwanted 
detections  are  filtered  in  the  signal  processing,  as  they  give  seismic  velocities  that  are  outside 
accepted  limits  for  regional  and  teleseismic  phase  velocities. 


Events  automatically  located  by  the  Apatity  array 

During  days  001, 2002,  through  181, 2002, 1,872  local  and  regional  events  were  located  by  the 
Apatity  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  10.3  events  per  processed  day  (181  days  processed).  42%  of  these  events  are  within  300  km, 
and  76%  of  these  events  are  within  1000  km. 

U.  Baadshaug 
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3.5  Regional  Monitoring  System  Operation  and  Analysis 

The  Regional  Monitoring  System  (RMS)  was  installed  at  NORSAR  in  December  1989  and  has 
been  operated  at  NORSAR  from  1  January  1990  for  automatic  processing  of  data  from  ARCES 
and  NORES.  A  second  version  of  RMS  that  accepts  data  from  an  arbitrary  number  of  arrays 
and  single  3-component  stations  was  installed  at  NORSAR  in  October  1991,  and  regular  oper¬ 
ation  of  the  system  comprising  analysis  of  data  from  the  4  arrays  ARCES,  NORES,  FINES  and 
GERES  started  on  15  October  1991.  As  opposed  to  the  first  version  of  RMS,  the  one  in  current 
operation  also  has  the  capability  of  locating  events  at  teleseismic  distances. 

Data  from  the  Apatity  array  was  included  on  14  December  1992,  and  from  the  Spitsbergen 
array  on  12  January  1994.  Detections  from  the  Hagfors  array  were  available  to  the  analysts  and 
could  be  added  manually  during  analysis  fi-om  6  December  1994.  After  2  February  1995,  Hag¬ 
fors  detections  were  also  used  in  the  automatic  phase  association. 

Since  24  April  1999,  RMS  has  processed  data  from  all  the  seven  regional  arrays  ARCES, 
NORES,  FINES,  GERES  (until  January  2000),  Apatity,  Spitsbergen,  and  Hagfors.  Starting 
19  September  1999,  waveforms  and  detections  from  the  NORSAR  array  have  also  been  avail¬ 
able  to  the  analyst. 

Phase  and  event  statistics 

Table  3.5.1  gives  a  summary  of  phase  detections  and  events  declared  by  RMS.  From  top  to  bot¬ 
tom  the  table  gives  the  total  number  of  detections  by  the  RMS,  the  number  of  detections  that 
are  associated  with  events  automatically  declared  by  the  RMS,  the  number  of  detections  that 
are  not  associated  with  any  events,  the  number  of  events  automatically  declared  by  the  RMS, 
and  finally  the  total  number  of  events  worked  on  interactively  (in  accordance  with  criteria  that 
vary  over  time;  see  below)  and  defined  by  the  analyst. 

New  criteria  for  interactive  event  analysis  were  introduced  from  1  January  1994.  Since  that 
date,  only  regional  events  in  areas  of  special  interest  (e.g,  Spitsbergen,  since  it  is  necessary  to 
acquire  new  knowledge  in  this  region)  or  other  significant  events  (e.g,  felt  earthquakes  and 
large  industrial  explosions)  were  thoroughly  analyzed.  Teleseismic  events  of  special  interest 
are  also  analyzed. 

To  further  reduce  the  workload  on  the  analysts  and  to  focus  on  regional  events  in  preparation 
for  Gamma-data  submission  during  GSETT-3,  a  new  processing  scheme  was  introduced  on  2 
February  1995.  The  GBF  (Generalized  Beamforming)  program  is  used  as  a  pre-processor  to 
RMS,  and  only  phases  associated  with  selected  events  in  northern  Europe  are  considered  in  the 
automatic  RMS  phase  association.  All  detections,  however,  are  still  available  to  the  analysts 
and  can  be  added  manually  during  analysis. 
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Jan 

02 

Feb 

02 

Mar 

02 

Apr 

02 

May 

02 

Jun 

02 

Total 

Phase  detections 

115917 

161223 

930267 

-  Associated  phases 

4318 

3513 

6828 

5245 

6576 

4226 

30706 

-  Unassociated  phases 

111599 

155978 

208509 

166893 

899561 

Events  automatically 
declared  by  RMS 

783 

646 

1388 

976 

1287 

855 

5935 

No.  of  events  defined  by 
the  analyst 

76 

76 

82 

68 

66 

66 

434 

Xflblo  3.S>1.  RMS  phase  detections  and  event  summary. 

U.  Baadshaug 
B.  Paulsen 
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4  NDC  and  Field  Activities 

4.1  NDC  Activitities 

NORSAR  functions  as  the  Norwegian  National  Data  Center  (NDC)  for  CTBT  verification.  Six 
monitoring  stations,  comprising  altogether  119  field  instruments,  will  be  located  on  Norwegian 
territory  as  part  of  the  future  IMS  as  described  elsewhere  in  this  report.  The  four  seismic  IMS 
stations  are  all  in  operation  today,  and  three  of  them  are  currently  providing  data  to  the  IDC. 
The  radionuclide  station  at  Spitsbergen  is  currently  operating  in  a  testing  mode,  whereas  the 
infrasound  station  in  northern  Norway  will  need  to  be  established  within  the  next  few  years. 
Data  recorded  by  the  Norwegian  stations  is  being  transmitted  in  real  time  to  the  Norwegian 
NDC,  and  provided  to  the  IDC  through  the  Global  Communications  Infrastructure  (GCI).  Nor¬ 
way  is  connected  to  the  GCI  with  a  frame  relay  link  to  Vienna. 

Operating  the  Norwegian  IMS  stations  continues  to  require  increased  resources  and  additional 
personnel  both  at  the  NDC  and  in  the  field.  The  PTS  has  estabhshed  new  and  strictly  defined 
procedures  as  well  as  increased  emphasis  on  regularity  of  data  recording  and  timely  data  trans¬ 
mission  to  the  IDC  in  Vienna.  This  has  led  to  increased  reporting  activities  and  implementation 
of  new  procedures  for  the  NDC  operators.  The  NDC  carries  out  all  the  technical  tasks  required 
in  support  of  Norway’s  treaty  obligations.  NORSAR  will  also  carry  out  assessments  of  events 
of  special  interest,  and  advise  the  Norwegian  authorities  in  technical  matters  relating  to  treaty 
comphance. 

Verification  functions;  information  received  from  the  IDC 

After  the  CTBT  enters  into  force,  the  IDC  will  provide  data  for  a  large  number  of  events  each 
day,  but  will  not  assess  whether  any  of  them  are  likely  to  be  nuclear  explosions.  Such  assess¬ 
ments  will  be  the  task  of  the  States  Parties,  and  it  is  important  to  develop  the  necessary  national 
expertise  in  the  participating  countries.  An  important  task  for  the  Norwegian  NDC  will  thus  be 
to  make  independent  assessments  of  events  of  particular  interest  to  Norway,  and  to  communi¬ 
cate  the  results  of  these  analyses  to  the  Norwegian  Ministry  of  Foreign  Affairs. 

Monitoring  the  Arctic  region 

Norway  will  have  monitoring  stations  of  key  importance  for  covering  the  Arctic,  including 
Novaya  Zemlya,  and  Norwegian  experts  have  a  unique  competence  in  assessing  events  in  this 
region.  On  several  occasions  in  the  past,  seismic  events  near  Novaya  Zemlya  have  caused  polit¬ 
ical  concern,  and  NORSAR  specialists  have  contributed  to  clarifying  these  issues. 

International  cooperation 

After  entry  into  force  of  the  treaty,  a  number  of  countries  are  expected  to  estabhsh  national 
expertise  to  contribute  to  the  treaty  verification  on  a  global  basis.  Norwegian  experts  have  been 
in  contact  with  experts  from  several  countries  with  the  aim  of  establishing  bilateral  or  multi¬ 
lateral  cooperation  in  this  field.  One  interesting  possibihty  for  the  future  is  to  establish 
NORSAR  as  a  regional  center  for  European  cooperation  in  the  CTBT  verification  activities. 
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NORSAR  event  processing 

The  automatic  routine  processing  of  NORSAR  events  as  described  in  NORSAR  Sci.  Rep.  No. 
2-93/94,  has  been  running  satisfactorily.  The  analyst  tools  for  reviewing  and  updating  the  solu¬ 
tions  have  been  continually  modified  to  simplify  operations  and  improve  results.  NORSAR  is 
currently  applying  teleseismic  detection  and  event  processing  using  the  large-aperture 
NORSAR  array  as  well  as  regional  monitoring  using  the  network  of  small-aperture  arrays  in 
Fennoscandia  and  adjacent  areas. 

Certification  of  PS28 

On  8  November  2001  the  IMS  station  PS28-ARCES  was  formally  certified.  PTS  personnel  vis¬ 
ited  the  station  during  the  winter  of  2000  and  carried  out  a  detailed  technical  evaluation.  As  a 
result  of  this  inspection  and  subsequent  discussions  between  NORSAR  and  the  PTS,  it  was 
concluded  that  PS28  needed  only  one  enhancement  in  order  to  be  certified:  to  install  a  central¬ 
ized  authentication  process  at  the  central  array  recording  facility.  After  this  was  done  during 
the  fall  of  2001  and  the  subsequent  verification  by  the  PTS,  station  certification  was  granted. 

Communication  topology 

Norway  has  implemented  an  independent  subnetwork,  which  connects  the  IMS  stations  AS72, 
AS73,  PS28,  and  RN49  operated  by  NORSAR  to  the  GCI  at  NOR_NDC.  A  contract  has  been 
concluded  and  VS  AT  antennas  have  been  installed  at  each  station  in  the  network.  Under  the 
same  contract,  VS  AT  antennas  for  6  of  the  PS27  subarrays  have  been  installed  for  intra-array 
communication.  The  seventh  subarray  is  connected  to  the  central  recording  facility  via  a  leased 
land  line.  The  central  recording  facility  for  PS27  is  connected  directly  to  the  GCI  (Basic 
Topology).  All  the  VS  AT  communication  is  functioning  satisfactorily. 

The  Norwegian  NDC  has  been  cooperating  with  institutions  in  other  countries  for  transmission 
of  IMS  data  to  the  Prototype  IDC  during  GSETT-3.  Details  on  this  can  be  found  in  Section  4.2. 

Jan  Fyen 


4.2  Status  Report:  Norway’s  Participation  in  GSETT-3 

Introduction 

This  contribution  is  a  report  for  the  period  January  -  June  2002  on  activities  associated  with 
Norway’s  participation  in  the  GSETT-3  experiment,  which  provides  data  to  the  International 
Data  Centre  (IDC)  in  Vienna  on  an  experimental  basis  until  the  participating  stations  have  been 
commissioned  as  part  of  the  International  Monitoring  System  (IMS)  network  defined  in  the 
protocol  to  the  Comprehensive  Nuclear-Test-Ban  Treaty.  This  report  represents  an  update  of 
contributions  that  can  be  found  in  previous  editions  of  NORSAR’s  Semiannual  Technical 
Summary.  It  is  noted  that  as  of  30  June  2002,  two  out  of  the  three  Norwegian  seismic  stations 
providing  data  to  the  IDC  have  been  formally  certified  and  thus  commissioned  as  part  of  the 
IMS  network. 
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Norwegian  GSETT-3  stations  and  communications  arrangements 

During  the  reporting  interval  1  January  -  30  June  2002,  Norway  has  provided  data  to  the 
GSETT-3  experiment  from  the  three  seismic  stations  shown  in  Fig.  4.2.1.  The  NORSAR  array 
(PS27,  station  code  NOA)  is  a  60  km  aperture  teleseismic  array,  comprised  of  7  subarrays,  each 
containing  six  vertical  short  period  sensors  and  a  three-component  broadband  instrument. 
ARCES  is  a  25-element  regional  array  with  an  aperture  of  3  km,  whereas  the  Spitsbergen  array 
(station  code  SPITS)  has  9  elements  within  a  1-km  aperture.  ARCES  and  SPITS  both  have  a 
broadband  three-component  seismometer  at  the  array  center. 

The  intra-array  communication  for  NOA  utihzes  a  land  line  for  subarray  NC6  and  VS  AT  links 
based  on  TDMA  technology  for  the  other  6  subarrays.  The  central  recording  facility  for  NOA 
isatNOR_NDC. 

Continuous  ARCES  data  has  been  transmitted  from  the  ARCES  site  to  NOR_NDC  using  a  64 
kbits/s  VS  AT  satellite  link,  based  on  BOD  technology. 

Continuous  SPITS  data  has  been  transmitted  to  NOR_NDC  via  a  VS  AT  terminal  located  at 
Plataberget  in  Longyearbyen  (which  is  the  site  of  the  IMS  radionuchde  monitoring  station 
RN49  installed  during  2001). 

Seven-day  station  buffers  have  been  established  at  the  ARCES  and  SPITS  sites  and  at  all  NOA 
subarray  sites,  as  well  as  at  NOR_NDC  for  ARCES,  SPITS  and  NOA. 

The  NOA  and  ARCES  arrays  are  primary  stations  in  the  GSETT-3  network  and  the  IMS,  which 
implies  that  data  from  these  stations  is  transmitted  continuously  to  the  receiving  international 
data  center.  Since  October  1999,  this  data  has  been  transmitted  (from  NOR_NDC)  via  the  Glo¬ 
bal  Communications  Infrastructure  (GCI)  to  the  IDC  in  Vienna,  whereas  transmission  of  the 
same  data  to  the  Prototype  International  Data  Center  (IDC)  in  Arlington,  VA,  was  discontinued 
on  7  February  2000.  The  SPITS  array  is  an  auxihary  station  in  GSETT-3  and  the  IMS,  and  the 
SPITS  data  have  been  available  to  both  the  IDC  and  the  IDC  throughout  the  reporting  period  on 
a  request  basis  via  use  of  the  AutoDRM  protocol  (Kradolfer,  1993;  Kradolfer,  1996).  The  Nor¬ 
wegian  stations  are  thus  participating  in  GSETT-3  with  the  same  status  (primary/auxiliary  seis¬ 
mic  stations)  they  have  in  the  IMS  defined  in  the  protocol  to  the  Comprehensive  Nuclear-Test- 
Ban  Treaty.  In  addition,  continuous  data  from  all  three  arrays  is  transmitted  to  the  US  NDC. 

Starting  with  this  report,  data  availabihties  for  primary  stations  and  the  number  of  requests 
made  for  auxiliary  station  data  are  reported  with  reference  to  the  IDC,  as  opposed  to  the  PIDC, 
as  was  done  in  previous  reports. 

Uptimes  and  data  availability 

Figs.  4.2.2  -  4.2.3  show  the  monthly  uptimes  for  the  Norwegian  GSETT-3  primary  stations 
ARCES  and  NOA,  respectively,  for  the  period  1  January  -  30  June  2002,  given  as  the  hatched 
(taUer)  bars  in  these  figures.  These  barplots  reflect  the  percentage  of  the  waveform  data  that  are 
available  in  the  NOR_NDC  tape  archives  for  these  two  arrays.  The  downtimes  inferred  from 
these  figures  thus  represent  the  cumulative  effect  of  field  equipment  outages,  station  site  to 
NOR_NDC  communication  outage,  and  NOR_NDC  data  acquisition  outages. 

Figs.  4.2.2-4.2.3  also  give  the  data  availabihty  for  these  two  stations  as  reported  by  the  IDC  in 
the  IDC  Station  Status  reports.  The  main  reason  for  the  discrepancies  between  the  NOR_NDC 
and  IDC  data  availabilities  as  observed  from  these  figures  is  the  difference  in  the  ways  the  two 
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data  centers  report  data  availability  for  arrays:  Whereas  NOR_NDC  reports  an  array  station  to 
be  up  and  available  if  at  least  one  channel  produces  useful  data,  the  IDC  uses  weights  where  the 
reported  availability  (capabihty)  is  based  on  the  number  of  actually  operating  channels. 

Use  of  the  AutoDRM  protocol 

NOR_NDC’s  AutoDRM  has  been  operational  since  November  1995  (Mykkeltveit  &  Baads- 
haug,  1996).  The  monthly  number  of  requests  by  the  IDC  for  SPITS  data  for  the  period  January 
-  June  2002  is  shown  in  Fig.  4.2.4. 

NDC  automatic  processing  and  data  analysis 

These  tasks  have  proceeded  in  accordance  with  the  descriptions  given  in  Mykkeltveit  and 
Baadshaug  (1996).  For  the  period  January  -  June  2002,  NOR_NDC  derived  information  on  464 
supplementary  events  in  northern  Europe  and  submitted  this  information  to  the  Finnish  NDC  as 
the  NOR_NDC  contribution  to  the  joint  Nordic  Supplementary  (Gamma)  Bulletin,  which  in 
turn  is  forwarded  to  the  IDC.  These  events  are  plotted  in  Fig.  4.2.5. 

Data  forwarding  for  GSETT-3  stations  in  other  countries 

NOR—NDC  continued  to  provide  communications  for  the  GSETT-3  auxiliary  station  at  Nilore, 
Pakistan,  through  a  VS  AT  satellite  link  between  NOR_NDC  and  Pakistan’s  NDC  in  Nilore. 
The  IDC  obtains  data  from  the  Hagfors  array  (HFS)  in  Sweden  through  requests  to  the  Auto¬ 
DRM  server  at  NOR_NDC  (in  the  same  way  requests  for  Spitsbergen  array  data  are  handled, 
see  above).  Fig.  4.2.6  shows  the  monthly  number  of  requests  for  HFS  data  from  the  two  IDC 
accounts  “pipeline”  and  “testbed”. 

Current  developments  and  future  plans 

is  continuing  the  efforts  towards  improving  and  hardening  all  critical  data  acquisi¬ 
tion  and  data  forwarding  hardware  and  software  components,  so  as  to  meet  future  requirements 
related  to  operation  of  IMS  stations  to  the  maximum  extent  possible. 

The  PrepCom  has  tasked  its  Working  Group  B  with  overseeing,  coordinating,  and  evaluating 
the  GSETT-3  experiment.  The  PrepCom  has  also  encouraged  states  that  operate  IMS- 
designated  stations  to  continue  to  do  so  on  a  voluntary  basis  and  in  the  framework  of  the 
GSETT-experiment  until  such  time  that  the  stations  have  been  certified  for  formal  inclusion  in 
IMS.  The  NOA  array  was  formally  certified  by  the  PTS  on  28  July  2000,  and  a  contract  with 
the  PTS  in  \fienna  currently  provides  partial  funding  for  operation  and  maintenance  of  this  sta¬ 
tion.  The  ARCES  array  was  formally  certified  by  the  PTS  on  8  November  2001.  A  contract  has 
been  signed  with  the  PTS  for  operation  and  maintenance  of  this  station,  with  an  effective  date 
of  1  January  2002.  Provided  that  adequate  funding  continues  to  be  made  available  (from  the 
PTS  and  the  Norwegian  Ministry  of  Foreign  Affairs),  we  envisage  continuing  the  provision  of 
data  from  all  Norwegian  IMS-designated  seismic  stations  without  interruption  to  the  IDC  in 
Vienna. 

U.  Baadshaug 
S.  Mykkeltveit 
J.  Fyen 
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Fig.  4.2.1.  The  figure  shows  the  locations  and  configurations  of  the  three  Norwegian  seismic  array 
stations  that  have  provided  data  to  the  GSETT-3  experiment  during  the  period  1  January  - 
30  June  2002.  The  data  from  these  stations  are  transmitted  continuously  and  in  real  time  to 
the  Norwegian  NDC  (NOR_NDC).  The  stations  NOA  and  ARCES  have  participated  in 
GSETT-3  as  primary  stations,  whereas  SPITS  has  contributed  as  an  auxiliary  station. 
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ARCES  data  availability  at  NDC  and  IDC 


Jan  Feb  Mar  Apr  May  Jun 
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Fig.  4.2.2.  The  figure  shows  the  monthly  availability  of  ARCES  array  data  for  the  period  January  - 
June  2002  atNOR_NDC  and  the  IDC.  See  the  text  for  explanation  of  differences  in  definition 
of  the  term  “data  availability”  between  the  two  centers.  The  higher  values  (hatched  bars)  rep¬ 
resent  the  NOR_NDC  data  availability. 


NOA  data  availability  at  NDC  and  IDC 


o 

-  =□ 

, 

WmmM 

-t  ^ 

^ V"  ? 

*  i  '  ^ 

4 

r  >  >  » 

o  . 

-  .  •  . 

1 

CO 

- 

■  t-; 

,* 

i  *4 

?  V  • 

f  1  1 

,  .  ■ 

i  r  i  -• 

l  ,  y  1 

t  f 

V  4  ' 

°  - 

■ 

i  V-  ;• 

t  4  t  i 

a 

■ 

• 

■ 

§ 

■ 

• 

-  .1 

• 

■ 

c 

u 

■ 

e 

/ 

■ 

®  p  _ 

»  *  r  5 

I 

'  i  ^  j 

}  i.  t  i 

t  1  , 
y  \  V  ' 

lilt 

I 

o 

■ 

■>  • 

%  V  ( 

4^  a 

>  »  i 

y..'  i 

■ 

cvj 

*  "  ^  1 
?  V  1  ^ 

>  ■-  r 

^  f  ' 

n  1  i 

\  ^ 
i  i  A  ! 

V  V 

U  1’ 

ii 

\  H  ' 

?  5-  i 

'r  t  5  ; 

L 

o  - 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

02 

Fig.  4.2.3.  The  figure  shows  the  monthly  availability  of  NORSAR  array  data  for  the  period  January  - 
June  2002  at  NOR_NDC  and  the  IDC.  See  the  text  for  explanation  of  differences  in  definition 
of  the  term  “data  availability”  between  the  two  centers.  The  higher  values  (hatched  bars)  rep¬ 
resent  the  NOR_NDC  data  availability. 


Fig.  4.2.4.  The  figure  shows  the  monthly  number  of  requests  received  by  NOR_NDCfrom  the  IDC for 
SPITS  waveform  segments  during  January  -  June  2002. 
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Reviewed  Supplementary  events 


Fig.  4.2.5.  The  map  shows  the  464  events  in  and  around  Norway  contributed  by  NOR_NDC  during 
January  -  June  2002  as  supplementary  (Gamma)  events  to  the  IDC,  as  part  of  the  Nordic  sup¬ 
plementary  data  compiled  by  the  Finnish  NDC.  The  map  also  shows  the  seismic  stations  used 
in  the  data  analysis  to  define  these  events. 


Fig.  4.2.6.  The  figure  shows  the  monthly  number  of  requests  received  by  NOR_NDCfrom  the  IDC for 
HFS  waveform  segments  during  January  -  June  2002. 
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4.3  Field  Activities 

The  activities  at  the  NORSAR  Maintenance  Center  (NMC)  at  Hamar  currently  include  work 
related  to  operation  and  maintenance  of  the  following  IMS  seismic  stations:  the  NOA  teleseis- 
mic  array  (PS27),  the  ARCES  array  (PS28)  and  the  Spitsbergen  array  (AS72).  Some  prepara¬ 
tory  work  has  also  been  carried  out  in  connection  with  the  seismic  station  on  Jan  Mayen 
(AS73),  the  inffasound  station  at  Karasjok  (IS37)  and  the  radionuclide  station  at  Spitsbergen 
(RN49).  NORSAR  also  acts  as  a  consultant  for  the  operation  and  maintenance  of  the  Hagfors 
array  in  Sweden  (AS  101). 

In  addition  to  the  above  activities,  which  are  directly  related  to  the  International  Monitoring 
System,  NORSAR’s  field  staff  are  continuing,  within  available  resources,  to  maintain  the 
small-aperture  NORES  array,  which  is  co-located  with  NOA  subarray  06C.  These  efforts  are 
given  low  priority,  since  there  is  no  requirement  for  specific  uptimes  at  NORES. 

NORSAR  carries  out  the  field  activities  relating  to  IMS  stations  in  a  manner  generally  consis¬ 
tent  with  the  requirements  specified  in  the  appropriate  IMS  Operational  Manuals,  which  are 
currently  being  developed  by  Working  Group  B  of  the  Preparatory  Commission.  For  seismic 
stations  these  specifications  are  contained  in  the  Operational  Manual  for  Seismological  Moni¬ 
toring  and  the  International  Exchange  of  Seismological  Data  (CTBTAVGB/TL-1 1/2),  currently 
available  in  a  draft  version. 

All  regular  maintenance  on  the  NORSAR  field  systems  is  conducted  on  a  one-shift-per-day, 
five-day-per-week  basis.  The  maintenance  tasks  include: 

•  Operating  and  maintaining  the  seismic  sensors  and  the  associated  digitizers,  authentication 
devices  and  other  electronics  components. 

•  Maintaining  the  power  supply  to  the  field  sites  as  well  as  backup  power  supplies. 

•  Operating  and  maintaining  the  VSATs,  the  data  acquisition  systems  and  the  intra-array  data 
transmission  systems. 

•  Assisting  the  NDC  in  evaluating  the  data  quality  and  making  the  necessary  changes  in  gain 
settings,  frequency  response  and  other  operating  characteristics  as  required. 

•  Carrying  out  preventive,  routine  and  emergency  maintenance  to  ensure  that  all  field  sys¬ 
tems  operate  properly. 

•  Maintaining  a  computerized  record  of  the  utilization,  status,  and  maintenance  history  of  all 
site  equipment. 

•  Providing  appropriate  security  measures  to  protect  against  incidents  such  as  intrusion,  theft 
and  vandalism  at  the  field  installations. 

Details  of  the  daily  maintenance  activities  are  kept  locally.  As  part  of  its  contract  with  CTBTO/ 
PTS  NORSAR  submits,  when  applicable,  problem  reports,  outage  notification  reports  and 
equipment  status  reports.  The  contents  of  these  reports  and  the  circumstances  under  which  they 
will  be  submitted  are  specified  in  the  draft  Operational  Manual. 

P.W.  Larsen 
K.A.  L0ken 
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6  Summary  of  Technical  Reports  /  Papers  Published 

6.1  Seismic  Event  Location  Calibration 

Report  from  the  IDC  Technical  Experts  Meeting  in  Oslo,  Norway  22-26  April  2002 

6.1.1  Introduction 

The  International  Data  Centre  (IDC)  Technical  Experts  Group  on  Seismic  Event  Location  held 
its  fourth  annual  meeting  in  Oslo,  Norway  on  22-26  April  2002.  The  meeting  was  held  jointly 
with  the  IDC  Technical  Experts  Group  on  Seismoacoustic  Event  Screening.  The  purpose  of  the 
meeting  was  to  support  the  ongoing  calibration  and  screening  efforts  of  the  IDC  and  in  particu¬ 
lar  to  review  progress  toward  developing  regionalized  travel  times  to  improve  the  quality  of 
location  estimates  of  seismic  events  reported  in  the  IDC  bulletins. 

Sixty-six  technical  experts,  coming  from  twelve  signatory  countries  and  the  Provisional  Tech¬ 
nical  Secretariat,  participated  in  the  meeting.  In  accordance  with  previous  recommendations, 
the  focus  of  the  discussions  was  on  the  following  geographical  regions:  North  America,  Eur¬ 
asia,  Northern  Africa  and  Australia.  Dr.  Erode  Ringdal  of  Norway  chaired  the  meeting. 

6.1.2  Background  and  technical  objectives 

Working  Group  B  has  repeatedly  encouraged  States  Signatories  to  support  the  location 
improvement  efforts  by  supplying  relevant  location  calibration  information  for  their  own  terri¬ 
tories  as  well  as  for  other  regions  where  they  have  such  information  available.  The  following 
types  of  calibration  information  were  proposed  in  the  document  CTBTAVGB-6/CRP.26: 

•  Precise  information  on  location,  depth,  and  origin  time  of  previous  nuclear  explosions  or 
large  chemical  explosions 

•  Similar  information  on  other  seismic  events  that  have  been  located  by  regional  networks 
with  sufficient  precision 

•  Data  as  appropriate  on  seismic  travel-time  models 

•  Any  other  information  (e.g.,  geologic  or  tectonic  maps)  that  would  be  useful 

•  Ground  truth  data  from  chemical  explosions. 

At  its  first  meeting  in  January  1999,  the  IDC  Technical  Experts  Group  on  Seismic  Event  Loca¬ 
tion  developed  plans  and  recommendations  for  a  global  calibration  program,  and  presented  its 
report  to  Working  Group  B  in  February  1999  (CTBTAVGB/TL-2/18).  This  work  was  reviewed 
and  updated  during  the  second  and  third  meetings  of  the  Experts  Group  in  March  2000  and 
April  2001,  and  the  results  were  subsequently  presented  to  Working  Group  B  (CTBTAV GB/ 
TL-2/49  and  CTBTAVGB/TL-2/61).  The  fourth  meeting  of  the  Experts  Group  (22-26  April 
2002)  had  the  following  objectives: 

•  To  review  new  developments  of  common  interest  to  location  calibration  and  event  screen¬ 
ing,  with  special  emphasis  on  event  depth  determination  and  computation  of  location  con¬ 
fidence  ellipses 

•  To  report  on  and  review  progress  of  ongoing  research  work  on  location  calibration,  includ¬ 
ing  calibration  consortia  and  PTS  Phase  1  calibration  contracts 
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To  review  proposals  for  detailed  station-specific  regional  location  corrections,  with  partic¬ 
ular  emphasis  on  IMS  stations  in  North  America,  Europe,  North  Africa,  Asia  and  Australia 

•  To  recommend  a  set  of  such  corrections,  including  appropriate  model  errors,  for  incorpora¬ 
tion  into  the  operational  IDC  software 

To  develop  a  plan  for  future  extensions  and  improvements  of  this  regional  correction  data 
base,  to  be  incorporated  into  future  IDC  software  releases 

•  To  review  progress  in  the  general  recommendations  from  the  first  and  second  meetings, 
and  make  adjustments  and  updates  to  these  recommendations  as  required. 

The  primary  task  of  the  meeting  was  to  assess  the  status  and  availability  of  such  calibration 
information  for  the  regions  being  considered,  and  to  plan  for  implementing  regional  location 
calibration  at  the  IDC  as  well  as  discuss  the  need  for  future  research  and  development. 


6.13  Technical  Issues 


Presentations  during  the  meeting 

A  number  of  papers  relating  to  the  collection,  application  and  validation  of  calibration  informa¬ 
tion  were  presented  by  participants.  Models  for  regionalization  on  a  global  basis  were  pre¬ 
sented  and  discussed.  Specific  presentations  were  made  by  several  experts  describing  regional 
velocity  models  and  calibration  data  for  the  general  geographic  regions  being  considered  ini- 
tiaUy.  Information  was  provided  about  the  current  CTBTO  Calibration  Programme.  Progress 
was  reported  at  the  workshop  by  U.S  sponsored  consortia,  by  CTBTO  sponsored  contractors 
and  by  several  other  research  groups. 


It  was  noted  that  for  some  regions,  information  was  incomplete  or  lacking,  and  the  use  of 
default  “generic”  velocity  models  for  various  tectonic  regions  was  discussed  in  some  detail. 
Valuable  new  data  on  ground  truth  (GT)  information  for  seismic  events  was  presented.  These 
data  wiU  be  organized  and  made  available  to  the  IDC  and  interested  States  Signatories.  Coun¬ 
tries  were  encouraged  to  continue  to  provide  relevant  calibration  data  for  the  purpose  of  devel¬ 
oping  accurate  seismic  travel-time  curves  for  various  geographical  regions. 


Reports  were  presented  on  a  number  of  modelling  studies,  some  of  which  showed  significant 
improvement  in  location  precision  when  applied  to  test  sets  of  seismic  events.  Three-dimen¬ 
sional  models  were  introduced  for  several  regions  and  were  found  to  provide  considerable 
improvements  in  location  accuracy  compared  to  standard  (IASPEI-91)  models. 

Techniques  for  improved  regional  processing  using  sparse  seismic  networks  as  well  as 
improved  azimuth  determination  for  regional  arrays  were  presented  and  discussed.  The  appli¬ 
cation  of  special  location  and  depth  estimation  techniques  was  also  addressed. 


Working  Group  Discussions 

Three  Working  Groups,  each  focusing  on  specific  regions  of  the  world,  were  established  to 
discuss  technical  issues  in  detail  during  the  workshop: 

Working  Group  1:  Northern  Eurasia  and  East  Asia 

Working  Group  2:  Southwestern  Asia  and  the  African/Mediterranean  area 

Working  Group  3:  North  America,  Australia,  Global  models 
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The  Working  Groups  were  given  a  mandate  with  a  list  of  specific  questions  addressing  the  fol¬ 
lowing  topics; 

Topic  1;  Validation  and  Implementation  of  Regional  Calibration  Information 
Topic  2:  Collection  of  Regional  Calibration  Information 
Topic  3:  Application  of  Regional  Calibration  Information 
Topic  4:  Future  work  of  the  Experts  Group 

The  results  of  the  Working  Groups  were  presented  and  discussed  in  a  plenary  session.  In  some 
cases,  previous  recommendations  were  reiterated  or  amplified.  These  presentations  and  discus¬ 
sions  provided  the  basis  for  the  recommendations  presented  below.  The  detailed  reports  of 
these  Working  Groups  are  available  on  request  from  the  Chairman  of  the  Experts  Group,  Dr. 
Frode  Ringdal,  Norway. 

6.1.4  Results  and  recommendations 
Main  results 

Since  the  2001  workshop,  participants  reported  considerable  progress  in  reference  event  data 
collection,  GT  criteria,  and  regional  calibration.  Reference  event  fists  have  significantly 
increased  over  the  last  year.  New  GT  category  criteria  have  been  proposed.  Regional  calibra¬ 
tion  has  demonstrated  reduced  bias  (absolute  errors)  and  decreased  uncertainty  (smaller  error 
eUipses)  in  accordance  with  the  goals  of  the  IMS  calibration  effort. 

Regional  corrections  have  been  implemented  at  the  IDC  for  IMS  stations  in  northwestern  Eur¬ 
asia  and  northern  America.  Work  is  continuing  on  developing  such  corrections  for  the  remain¬ 
ing  priority  regions,  and  encouraging  results  have  been  achieved. 

For  example,  for  much  of  Eastern  Asia,  preliminary  work  on  source-specific  station  corrections 
(SSSCs)  for  Pn  arrivals  has  given  very  promising  results  suggesting  that  it  will  be  possible  to 
achieve  or  exceed  the  uncertainty  goal  of  1000  km^.  For  example,  it  was  reported  that  reloca¬ 
tion  studies  using  Soviet  explosions  recorded  by  about  90  regional  stations  have  succeeded  in 
reducing  the  median  90%  error  ellipse  to  200  km^  for  these  explosions,  using  model-based 
analysis  combined  with  an  interpolation  approach  to  provide  preliminary  SSSCs.  More  limited 
validation  studies  conducted  using  a  sparse  IMS/surrogate  network  have  documented  improve¬ 
ments  in  location  accuracy  to  significantly  better  than  10  km.  Some  cautions  needs  to  be  exer¬ 
cised  in  extrapolating  these  results  to  smaller  events,  but  aU  indications  are  that  significant 
improvements  in  location  wiU  be  obtained  by  applying  SSSC  corrections  to  the  IMS  stations  in 
this  area. 

Validation  data  sets 

During  the  course  of  the  experts’  work  the  importance  of  validation  of  ground  truth  data  has 
become  increasingly  clear.  Validation  is  a  far  more  complex  and  time-consuming  undertaking 
than  initially  expected.  The  experts  recommend  that  a  systematic  development  of  a  validation 
data  set  be  undertaken.  This  validation  set  should  include  GT  events  of  various  categories,  and 
should  contain  all  relevant  information  (metadata)  about  the  events  in  the  data  base.  Before 
delivering  such  datasets  to  the  IDC  archive,  the  data  should  be  carefully  quality  controlled  by 
the  organization  providing  the  data.  Information  on  the  quality  control  of  origins  and  arrival 
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times  should  be  provided  to  the  IDC  along  with  the  data.  WGB  should  consider  arrangements 
to  make  these  data  available  from  the  IDC  without  restrictions  to  the  research  community. 

The  Expert  Group  re-emphasizes  the  need  for  a  formal  procedure  for  validation.  In  addition, 
there  should  be  standards  for  implementation  and  periodic  checking  of  performance.  The 
Group  recommends  that  global  models  (large-scale  regional  models)  continue  to  be  tested,  and 
points  out  the  need  for  a  high  resolution  crustal  model.  It  is  also  noted  that  there  are  currently 
no  teleseismic  SSSCs  available,  and  the  experts  recommend  that  such  corrections  be  devel¬ 
oped. 

Validation  metrics 

The  initial  set  of  validation  metrics  suggested  by  the  1999  Workshop  (CTBTAVGB/TL-2/18) 
should  be  updated.  The  experts  recommend  that  the  metrics  below  to  be  used  as  a  common 
ground  to  evaluate  the  performance  of  SSSCs: 

•  Average  and  median  improvement  in  epicentre  mislocation  for  all  relocated  events  with 
depth  fixed  to  the  GT  depth. 

•  Average  and  median  reduction  in  the  size  of  90%  error  ellipses  where  these  error  ellipses 
are  adjusted  to  GT  accuracy 

•  Actual  coverage  of  the  90%  confidence  region  (adjusted  to  GT  accuracy). 

•  For  events  where  accurate  origin  time  is  available: 

•  average  and  median  improvement  in  origin  time; 

•  origin  time  error  should  provide  90%  coverage 

•  Percentage  of  events  improved/deteriorated  with  respect  to  GT  epicentres 

•  Percentage  of  events  improved/deteriorated  by  more  than  20%  w.r.t.  GT  epicentres;  aver¬ 
age/median  improvement  and  deterioration  in  km 

The  experts  expect  that  additional  metrics  will  be  developed  over  the  next  few  years  which  may 
provide  more  robust  comparisons/evaluation  capability.  In  order  to  compare  corrections  that 
have  overlapping  geographical  coverage,  a  central  archive  should  be  established  at  the  IDC  to 
encourage  the  international  community  to  contribute  their  insights  and  evaluation. 

It  would  be  desirable  to  develop  a  simpler  “measure  of  performance”  that  would  indicate  which 
set  of  models/calibrations  are  statistically  “better”.  It  was  recognized  that  the  problem  has 
multiple  dimensions  and  location  performance  is  more  complicated  than  what  can  be  described 
by  a  single  number,  but  it  will  be  easier  to  communicate  to  non-experts  if  the  advantages/disad¬ 
vantages  of  multiple  overlapping  SSSCs  can  be  ranked  in  some  sense. 

Implementing  preliminary  corrections 

The  experts  recommend  that  the  preliminary  sets  of  SSSCs  that  are  now  available  should  be 
considered  for  testing  at  the  IDC  as  soon  as  possible  along  with  existing  available  sets  of  offline 
validation  test  sets.  Furthermore,  such  testing  should  include  computations  of  multiple  loca¬ 
tions  of  REB  events  for  open  evaluation  by  the  expert  community.  Multiple  models  should  be 
tested.  Input  and  results  should  be  available  in  an  open  and  transparent  database  with  results 
compared  to  current  practice.  All  experts  are  encouraged  to  provide  their  GT  events,  velocity 
models  and  SSSC  corrections,  including  supporting  information  (meta-data),  to  the  IDC  as 
soon  as  possible. 
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Transparency 

The  experts  again  expressed  concern  about  the  current  restrictions  on  obtaining  IMS  data  and 
IDC  products,  and  recommended  that  the  IDC  make  openly  available  to  the  scientists  involved 
in  the  IDC  location  calibration  effort  aU  of  the  waveform  data  and  associated  IDC  products  that 
are  needed  in  order  to  successfully  carry  out  the  caUbration  program.  The  experts  re-emphasize 
concerns  over  lack  of  transparency  in  IDC  bulletin  products  that  reflect  corrections  apphed  to 
travel  times,  azimuths,  and  slownesses  for  location.  A  mechanism  for  access  to  and  distribu¬ 
tion  of  corrections  stored  in  databases  and  flat-files  should  be  developed. 

Need  for  wider  participation 

The  experts  are  concerned  with  the  unfortunate  low  level  of  participation  in  calibration  activi¬ 
ties  in  under-represented  areas  such  as  Africa.  The  IDC  and  concerned  states  may  wish  to 
engage  in  programs  to  encourage  participation  in  such  areas.  Such  activities  might  include 
professional  exchanges  of  personnel  with  groups  actively  engaged  in  calibration  to  promote 
exchange  of  data  and  expertise.  In  recognition  of  the  importance  of  aftershock  surveys  in  the 
generation  of  valuable  reference  events,  other  activities  might  include  support  of  temporary 
aftershock  recordings  (instruments  and  personnel)  and  a  clearing-house  to  collect  aftershock 
data  and  maintain  an  open  database  of  aftershock  metadata  that  can  be  used  for  calibration  and 
reference  event  collection. 

Future  work 

The  IDC  contractors,  U.S.-sponsored  consortia  and  other  research  groups  are  expected  to 
dehver  to  the  IDC  a  variety  of  proposed  SSSCs  in  the  coming  year.  Evaluation  and  validation 
of  these  corrections  will  require  substantial  additional  research  study.  The  experts  consider  that 
merging  SSSCs  from  different  study  groups  with  different  methodologies  in  an  operational  sys¬ 
tem  presents  a  significant  challenge.  It  is  recommended  that  the  IDC  begin  planning  for  such  a 
continuing  research  program. 

Encouraging  results  were  presented  for  depth  phase  detection  and  identification.  Methods  for 
confident  detection  and  identification  of  depth  phases  remain  an  important  problem,  and 
research  in  this  area  should  continue.  Focused  discussion  of  selected  topics  such  as  depth  esti¬ 
mation  in  the  full  assembly  of  experts  is  encouraged. 

The  experts  consider  that  investigation  is  needed  into  the  possibihty  that  analyst  time  picks  (as 
well  as  automatic  onset  estimates)  are  late  for  low  signal-to-noise  ratio  recordings.  Such  biases 
would  result  in  different  apparent  basehnes  for  regional  and  teleseismic  phases.  More  gener¬ 
ally,  presentation  and  discussion  of  IDC  analyst  procedures  and  of  the  problems  facing  the  ana¬ 
lysts  should  be  considered  for  the  future. 

Collection  of  a  set  of  reference  events  will  continue  to  be  a  priority  with  emphasis  on  precise 
hypocenters  and  origin  times  and  a  good  geographic  coverage.  The  reference  events  should  be 
chosen  so  as  to  keep  usage  of  surrogate  (non-IMS)  stations  to  a  minimum.  For  clusters  of 
events  it  is  recommended  that  a  few,  well-constrained  events  be  chosen  (using  the  smallest  con¬ 
fidence  eUipse  as  a  selection  criterion).  Reference  events  should  be  recorded  regionally  and 
should  comprise  a  range  of  magnitudes. 

These  recommendations  wiU  be  considered  before  the  next  meeting  of  the  Experts  Group. 

F.  Ringdal 
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6.2  Research  in  regional  seismic  monitoring 

(Paper  presented  at  the  24th  Annual  Seismic  Research  Review) 

Abstract 

This  project  represents  a  continuing  effort  aiming  at  three  main  topics:  (a)  to  carry  out  research 
in  regional  monitoring  of  tiie  European  Arctic,  (b)  to  apply  experimental  methods  such  as  site- 
specific  threshold  monitoring  to  target  areas  of  interest  and  assess  the  results  and  (c)  to  contrib¬ 
ute  to  the  global  location  calibration  effort  currently  being  undertaken  by  the  Preparatory  Com¬ 
mission  (PrepCom)  of  the  Comprehensive  Nuclear  Test-Ban-Treaty  Organization’s  Working 
Group  B  in  Vienna,  Austria. 

We  have  used  data  from  the  regional  networks  operated  by  NORSAR  and  the  Kola  Regional 
Seismological  Centre  (KRSC)  to  assess  the  seismicity  and  characteristics  of  regional  phases  of 
tile  European  Arctic.  An  especially  interesting  recent  seismic  event,  which  has  been  analyzed 
in  detail,  occurred  on  Novaya  Zemlya  on  23  February  2002.  This  small  event  (magnitude  about 
3.0)  is  the  first  seismic  event  detected  in  or  near  Novaya  Zemlya  since  the  two  Kara  Sea  events 
of  16  August  1997.  The  event  was  recorded  with  particularly  high  SNR  by  the  International 
Monitoring  System  (IMS)  arrays  Spitsbergen  and  ARCES,  and  by  the  Amderma  station  (oper¬ 
ated  by  KRSC)  south  of  Novaya  Zemlya.  We  have  located  the  event  at  the  eastern  coast  of 
Novaya  Zemlya,  approximately  100  km  NE  of  the  former  nuclear  test  site.  Our  attempts  to  esti¬ 
mate  the  depth  of  this  event  proved  inconclusive,  and  a  depth  of  0  could  not  be  ruled  out.  We 
present  comparisons  of  the  recordings  of  this  event  to  recordings  of  previous  small  events  in 
the  Novaya  Zemlya  region. 

We  have  begun  an  effort  to  develop  an  optimized  site-specific  threshold  monitoring  system  for 
the  Lop  Nor  test  site  in  China.  Our  aim  is  to  calibrate  all  available  high-quality  seismic  stations 
within  reponal  distance,  and  supplement  with  the  best  IMS  arrays  at  teleseismic  distances.  As 
a  calibration  data  base,  we  have  used  several  past  nuclear  explosions  at  the  test  site,  combined 
with  nearby,  well-recorded  earthquakes.  Even  though  most  events  were  recorded  by  only  a  few 
stations,  this  data  set  was  sufficient  to  obtain  a  reasonable  amplitude  calibration  for  the  pur¬ 
poses  of  this  study. 

Initial  results  show  that,  among  the  available  stations,  the  IMS  Makanchi  array  in  Kazakhstan 
(at  an  epicentral  distance  of  6  degrees)  is  the  most  sensitive  station  for  detecting  events  at  Lop 
Nor,  followed  by  some  of  the  best  teleseismic  arrays,  for  example  the  NORES  array  in  Norway. 
Combining  the  data  by  using  the  threshold  monitoring  technique,  we  have  found  that  the  Lop 
Nor  site  can  be  monitored  down  to  a  threshold  close  to  mb  3.2.  This  threshold  is  about  0.8  mag¬ 
nitude  unit  higher  than  the  corresponding  capability  for  the  Novaya  Zemlya  test  site,  for  which 
the  highly  sensitive  arrays  ARCES  and  Spitsbergen  are  both  at  a  distance  of  only  10  degrees. 
We  note  that  no  stations  inside  China  were  selected  for  the  Lop  Nor  study. 

A  workshop  was  held  in  Oslo,  Norway,  during  22-26  April  2002  in  support  of  the  global  seis¬ 
mic  event  location  calibration  effort  currently  being  undertaken  by  PrepCom’s  Working  Group 
B  in  Vienna.  The  workshop,  which  was  chaired  by  Dr.  Erode  Ringdal,  was  attended  by  66  sci¬ 
entists  from  12  countries  and  the  Provisional  Technical  Secretariat  of  the  CTBTO.  The  work¬ 
shop  was  held  jointly  by  the  IDC  Experts  Groups  on  location  calibration  and  event  screening, 
and  the  results  and  recommendations  were  reported  to  Working  Group  B. 
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6.2.1  Objective 

This  work  represents  a  continued  effort  in  seismic  monitoring,  with  emphasis  on  studying  earth¬ 
quakes  and  explosions  in  the  Barents/Kara  Sea  region,  which  includes  the  Russian  nuclear  test 
site  at  Novaya  Zemlya.  The  overall  objective  is  to  characterize  the  seismicity  of  this  region,  to 
investigate  the  detection  and  location  capability  of  regional  seismic  networks  and  to  study  var¬ 
ious  methods  for  screening  and  identifying  seismic  events  in  order  to  improve  monitoring  of  the 
Comprehensive  Test  Ban  Treaty.  Another  objective  is  to  apply  advanced  site-specific  seismic 
monitoring  methods  to  other  sites  of  special  interest,  in  particular  known  nuclear  test  sites.  A 
third  objective  is  to  support  the  international  effort  to  provide  regional  location  calibration  of 
the  International  Monitoring  System. 

6.2.2  Research  Accomplished 
Introduction 

NORSAR  and  Kola  Regional  Seismological  Centre  (KRSC)  of  the  Russian  Academy  of  Sci¬ 
ences  have  for  many  years  cooperated  in  the  continuous  monitoring  of  seismic  events  in  North- 
West  Russia  and  adjacent  sea  areas.  The  research  has  been  based  on  data  from  a  network  of 
sensitive  regional  arrays  which  has  been  installed  in  northern  Europe  during  the  last  decade  in 
preparation  for  the  CTBT  monitoring  network.  This  regional  network,  which  comprises  sta¬ 
tions  in  Fennoscandia,  Spitsbergen  and  NW  Russia  provides  a  detection  capabihty  for  the  Bar¬ 
ents/Kara  Sea  region  that  is  close  to  m^  =  2.5  (Ringdal,  1997). 

The  research  carried  out  during  this  effort  is  documented  in  detail  in  several  contributions  con¬ 
tained  in  the  NORSAR  Semiannual  Technical  Summaries.  In  the  present  paper  we  will  limit 
the  discussions  to  some  recent  results  of  interest  in  the  general  context  of  regional  monitoring 
of  seismic  events  in  the  European  Arctic,  a  presentation  of  an  initial  threshold  monitoring  study 
for  the  Lop  Nor  test  site  in  China  and  a  brief  review  the  location  cahbration  effort  currently 
underway  for  the  International  Monitoring  System  (IMS). 

Low-magnitude  seismic  events  in  or  near  Novaya  Zemlya 

We  have  updated  our  previous  smdies  of  low-magnitude  events  in  the  Novaya  Zemlya  region 
recorded  by  the  NORSAR  regional  seismic  network.  Table  6.2.1  and  Fig.  6.2.1  summarize  the 
information  on  such  events  detected  since  1980.  In  spite  of  the  low  threshold  for  detection, 
only  8  seismic  events  with  estimated  location  outside  the  test  site  have  been  recorded  during 
these  more  than  20  years.  Furthermore,  we  have  since  1992  apphed  the  site-specific  threshold 
monitoring  network  to  the  test  site,  and  have  thereby  been  able  to  confirm  that  no  other  seismic 
event  of  magnitude  2.5  or  greater  has  occurred  in  this  region  since  then. 

An  especially  interesting  recent  seismic  event,  which  has  been  analyzed  in  detail,  occurred  on 
Novaya  Zemlya  on  23  February  2002.  This  small  event  (magnimde  about  3.0)  is  the  first  seis¬ 
mic  event  detected  in  or  near  Novaya  Zemlya  since  the  two  Kara  Sea  events  of  16  August  1997 . 
The  event  was  recorded  with  particularly  high  SNR  by  the  International  Monitoring  System 
(IMS)  arrays  Spitsbergen  (SPITS)  and  ARCES,  and  by  the  Amderma  station  (AMD),  operated 
by  KRSC,  south  of  Novaya  Zemlya.  We  have  located  the  event  at  the  eastern  coast  of  Novaya 
Zemlya,  approximately  100  km  NE  of  the  former  nuclear  test  site. 
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For  this  recent  event,  as  well  as  other  events  recorded  by  the  regional  network,  we  have  fre¬ 
quently  noted  uncertainties  in  onset  measurements  amounting  to  1-2  seconds.’ This  is  in  fact  a 
general  problem,  and  is  due  partly  to  the  low  SNR  for  some  phases,  partly  to  the  emergent 
nature  of  some  signal  recordings.  We  have  conducted  a  study  to  assess  the  effects  of  such 
uncertainty  on  location  and  depth  estimates,  using  a  number  of  seismic  events  in  the  European 
Arctic  (Novaya  Zemlya,  Kara  Sea,  Northern  Norway,  Kola  Peninsula).  For  each  event,  and 
each  detected  phase,  we  selected  up  to  three  onset  time  picks,  spanning  the  estimated  uncer¬ 
tainty  (typically  +/- 1-2  seconds).  We  then  relocated  each  event  using  all  possible  combinations 
of  such  time  picks,  and  compiled  all  solutions  except  those  which  showed  clear  inconsisten¬ 
cies.  In  this  way,  we  obtained  a  large  number  of  locations  and  depths  for  each  event,  and  made 
plots  to  indicate  the  distribution  of  both  the  epicenters  and  the  estimated  depths. 

The  results  for  the  23  February  2002  event  near  Novaya  Zemlya  are  illustrated  in  Fig.  6.2.2.  We 
note  that  the  epicentral  error  area  is  similar  to  the  90%  confidence  ellipse  in  Fig.  6.2.1.  The 
depth  distribution  is  inconclusive,  although  most  solutions  have  an  estimated  depth  near  0. 
From  the  general  study,  we  note  that  for  all  the  analyzed  events,  except  for  one  known  earth¬ 
quake  in  Northern  Norway,  a  depth  of  zero  cannot  be  ruled  out.  We  consider  that  this  should 
serve  as  a  caution  against  using  estimated  depth  as  a  screening  criterion  in  this  region.  In  par¬ 
ticular,  we  note  that  depth  estimation  using  regional  velocity  models  depend  very  strongly  on 
the  assumed  Yp/Vs  ratio. 

Waveforms  for  the  23  February  2002  event  from  a  number  of  seismic  stations  at  regional  dis¬ 
tances  from  Novaya  Zemlya  are  shown  in  Fig.  6.2.3.  The  high  signal-to-noise  ratios  indicates 
that  this  event  is  well  above  the  detection  threshold,  especially  for  the  Amderma,  SPITS  and 
ARCES  stations.  Note  that  by  array  beamforming,  the  SNR  at  ARCES  and  SPITS  could  be 
improved  by  more  than  0.5  magnitude  units.  This  again  reconfirms  the  capability  of  this 
regional  network  to  monitor  seismic  events  at  Novaya  Zemlya  down  to  a  magnitude  level  of 
2.0-2.5,  depending  on  the  actual  background  noise  levels. 

In  the  context  of  event  screening,  which  is  an  important  topic  for  the  CTBT  International  Data 
Center,  it  is  clear  from  these  results  that  no  screening  should  take  place  for  a  region  with  as  low 
seismic  activity  as  Novaya  Zemlya.  In  fact,  any  seismic  event  in  this  region  is  of  potential  mon¬ 
itoring  interest,  and  with  a  frquency  of  occurrence  of  magnitude  2.5  or  larger  at  less  than  one 
event  per  year,  the  actual  analysis  work  would  be  quite  modest.  Similarly,  as  noted  by  Ringdal 
et.  al.  (20{X)),  the  much  larger  geographical  area  of  western  Russia  and  the  Barents  Sea  is  also 
characterized  by  modest  natural  seismic  activity,  and  could  be  exempted  from  event  screening, 
at  least  at  the  current  magnitude  threshold  of  3.5  for  screening  purposes. 
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Table  6.2.1.  Seismic  events  in  or  near  Novaya  Zemlya  (1980-2002)  located  outside  the 

test  site. 


Date/time 

Location 

mb 

Comment 

01.08.86/ 

13.56.38 

72.945  N,  56.549  E 

4.3 

Located  by  Marshall  et.al.  (1989) 

31.12.92/ 

09.29.24 

73.600  N  55.200  E 

2.7 

Located  by  NORSAR 

23.02.95/ 

21.50.00 

71.856  N,  55.685  E 

2.5 

Located  by  NORSAR 

13.06.95/ 

19.22.38 

75.170  N,  56.740  E 

3.5 

Located  by  NORSAR 

13.01.96/ 

17.17.23 

75.240  N,  56.660  E 

2.4 

Approximately  co-located  with  preced¬ 
ing  event 

16.08.97/ 

02.11.00 

72.510  N,  57.550  E 

3.5 

Located  by  NORSAR 

16.08.97/ 

06.19.10 

72.510  N,  57.550  E 

2.6 

Co-located  with  preceding  event 

23.02.02/ 

01.21.14 

74.047  N,  57.671  E 

3.0 

Located  by  NORSAR 
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Fig.  6.2.1.  Seismic  events  detected  by  the  NORSAR  regional  array  system  in  or  near  Novaya  Zemlya 
during  J 980-2002.  Only  events  with  estimated  location  outside  the  former  nuclear  test  site  are 
included.  Estimated  90%  error  ellipses  are  indicated. 


Fig.  6.2.2.  Estimates  of  the  location  of  the  23.02.2002  seismic  event,  as  explained  in  the  text.  Note 
that  the  epicenters  corre.spond  well  with  the  estimated  uncertainty  ellipse  of  Fig.  6.2.1,  and 
that  the  depth  of  the  event  cannot  he  conclusively  determined. 
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23,GS.200S  e^ent  recordings*  filtered  4-8  Hz 


Fig.  6.2.3.  Recorded  SPZ  waveforms  for  the  23.02.2002  seismic  event.  The  data  are  filtered  in  the  4- 
8  Hz  band.  Note  the  high  P -phase  SNR  in  paricular  at  the  Spitsbergen  array  (sensor  SPAO) 
and  the  Ainderma  station  (AMDD).  Note  also  the  relatively  weak  S-phase  for  the  two  stations 
at  Spitsbergen  (SPAO  and  KBS). 

Khibiny  mines 

We  have  continued  our  research  on  rockbursts  and  mining  explosions  in  the  mining  areas  of 
NW  Russia,  in  particular  the  Khibiny  Massif.  Recently,  seismic  instrumentation  has  been 
installed  inside  the  mines  in  the  Khibiny  Massif  of  the  Kola  peninsula  in  order  to  provide  origin 
times  of  the  seismic  events  as  well  as  to  contribute  to  additional  vahdation  of  the  location  accu¬ 
racy.  Fig.  6.2.4  shows  an  example  of  recordings  for  a  mining  explosion  of  about  magnitude  2.0. 

We  are  also  cooperating  with  Lawrence  Livermore  National  Laboratory  in  a  DOE-funded 
project  to  carry  out  more  detailed  studies  of  the  characteristics  of  recordings  from  mining 
events  in  nothem  Fennoscandia  and  Western  Russia.  That  project  includes  the  installation  of 
additional  seismometers  along  profiles  in  Norway,  Finland  and  the  Kola  Peninsula,  for  record¬ 
ing  over  a  period  of  one  year. 
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5. 19:08  4 


Explosion  in  Kirovsk  mine 


Fig.  6.2.4.  Example  of  recordings  of  a  mining  explosion  in  the  Khibiny  Massif  Kola  Peninsula,  using 
a  network  of  recently  installed  local  stations  as  well  as  stations  at  regional  distances. 

Threshold  monitoring  of  the  Lop  Nor,  China  test  site 

As  part  of  the  DoD  Advanced  Concept  Demonstration,  we  have  applied  the  Site-Specific 
Threshold  Monitoring  (SSTM)  technique  to  the  Lop  Nor  test  site  in  China  (Lindholm  et.  al., 
2002).  The  emphasis  has  been  on  detection  (and  location)  of  small  seismic  events  with  mb  < 
4.0,  and  the  purpose  is  to  evaluate  the  SSTM  method  as  a  potential  monitoring  tool  for  this  site. 
In  contrast  to  most  previous  case  studies,  which  have  been  based  on  recordings  by  seismic 
arrays  at  regional  distances,  we  have  in  this  study  applied  a  combination  of  3-component  sta¬ 
tions  and  arrays,  at  both  regional  and  teleseismic  ranges. 

Our  efforts  so  far,  as  reported  in  this  contribution,  comprises  mainly  a  study  of  available  seis¬ 
mic  stations,  selection  of  those  stations  which  are  most  sensitive  to  seismic  events  in  the  Lop 
Nor  general  area,  and  tuning  of  the  signal  parameters  of  these  stations  so  as  to  prepare  process¬ 
ing  recipes  for  the  application  of  the  threshold  monitoring  tool. 

Development  of  processing  recipes 

For  the  successful  implementation  of  SSTM,  beams  filtered  in  optimal  frequency  bands  must 
be  steered  from  the  individual  arrays  towards  the  Lop  Nor  test  site,  and  amplitude  calibration 
constants  developed  from  older  events  are  applied  to  facilitate  the  calculation  of  continuous 
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magnitude  thresholds  for  the  site.  For  single  (or  3-component)  stations,  the  vertical  component 
is  filtered  in  the  optimal  frequency  band. 

For  the  purposes  of  this  study,  23  single  stations  and  arrays  were  initially  chosen  as  shown  in 
Fig.  6.2.5.  For  each  of  the  stations  data  were  collected  from  known  nuclear  tests  at  Lop  Nor, 
and  supplemented  with  data  for  low-magnitude  earhquakes  near  the  test  site.  We  note  that  for 
most  stations  only  a  few  data  sets  with  explosions  were  available  for  the  calibration. 


Fig.  6.2.5.  Network  selected  for  threshold  monitoring  of  the  Lop  Nor,  China  test  site.  The  test  site  is 
marked  with  a  star  symbol. 


6.2.3  Preliminary  Results 

Two  performance  tests  were  carried  out.  First,  data  for  one  day  (September  10,  2001)  was  col¬ 
lected,  and  in  these  data  the  recordings  of  4  explosions  were  scaled  and  embedded  at  some  (but 
not  all  of  the)  stations.  The  explosions  were  scaled  to  magnitude  3.5  and  3.0,  so  that  a  total  of  8 
embedded  explosions  were  included.  Secondly,  a  10  day  test  period  with  data  from  August  2 
through  1 1,  2001  was  selected  and  data  were  collected  for  all  available  stations  in  our  network. 
Fig.  6.2.6  shows  examples  of  network  results  of  these  two  performance  tests. 

From  these  tests  the  following  preliminary  observations  can  be  stated: 

•  Out  of  the  eight  embedded  explosions  on  day  253  six  were  detected  automatically  on  the 
threshold  trace,  and  the  remaining  two  could  be  clearly  seen  on  the  trace,  even  if  the  auto¬ 
matic  detector  did  not  trigger. 

•  The  quiet  day  (214)  did  not  have  any  events  flagged. 
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•  On  day  217  the  SSTM  method  triggered  on  two  large  teleseismic  events. 

•  The  network  threshold  is  close  to  or  better  than  3.0  (based  on  both  arrays  and  single  sta¬ 
tions).  The  trigger  level  (shown  as  a  continuous,  smooth  line  in  Fig.  6.2.6)  is  at  magnitude 
3.2  or  lower. 

•  In  a  monitoring  situation,  all  peaks  exceeding  the  threshold  could  be  analyzed.  The  maxi¬ 
mum  number  of  such  peaks  during  one  day  was  10-15  for  the  days  processed.  If  such  peaks 
were  analyzed,  all  of  the  scaled  events  on  day  2001-253  would  have  been  found. 

The  above  observations  should  be  evaluated  under  the  perspective  that  the  calibration  is  a  pre¬ 
liminary  one,  where  the  majority  of  the  stations  could  only  be  calibrated  with  one  or  two 
events/explosions,  and  where  only  some  of  the  calibrated  stations  were  included  in  the  perfor¬ 
mance  test.  Furthermore,  the  day  with  embedded  events  did  not  include  embedded  data  for  all 
the  stations,  and  this  did  reduce  the  performance  of  the  threshold  monitoring. 

In  the  paper  by  Lindholm  et.  al.  (2002),  there  is  also  some  discussion  of  the  relative  perfor¬ 
mance  of  individual  stations  in  the  network.  The  SSTM  threshold  traces  for  single  stations  and 
arrays  can  in  fact  be  used  to  assess  this  performance.  We  have  observed  that  the  arrays  gener¬ 
ally  have  thresholds  between  magnitude  3.2  and  3.8,  and  that  the  single  stations  vary  consider¬ 
ably,  from  about  3.5  for  the  best  (closest)  stations  to  almost  4.5  for  some  of  the  lower 
performing  stations.  Some  stations  in  the  data  set  feature  frequent  data  problems. 

As  the  IMS  system  is  further  developed,  there  is  a  clear  potential  for  imrovement  in  the  perfor¬ 
mance  illustrated  above.  For  example,  the  sensitive  array  in  Makanchi,  Kazakhstan,  was  not 
available  to  us  for  the  test  period  in  August/September  2001 .  Studies  of  the  performance  of  this 
array  has  indicated  that  it  has  a  superior  capability  for  detecting  events  at  Lop  Nor  compared  to 
the  other  stations  in  the  network,  and  therefore  would  contribute  to  significantly  improved 
thresholds,  once  included. 

We  also  note  that  the  current  study  has  been  done  without  including  any  Chinese  IMS  stations 
in  our  experimental  network.  Including  such  stations  would  no  doubt  contribute  to  improved 
thresholds,  but  there  are  other  considerations  that  may  make  it  useful  to  evaluate  how  the  test 
site  can  be  monitored  using  only  station  outside  the  country. 
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Fig.  6.2.6.  Threshold  monitoring  network  traces  from  three  different  days.  The  upper  trace  is  from 
day  253  with  eight  embedded  Lop  Nor  explosion  recordings  scaled  to  magnitudes  3.0  and  3.5. 
The  ellipses  indicate  the  times  of  the  embedded  explosions.  For  six  of  the  eight  embedded 
explosions  the  SSTM  automatic  analysis  picked  the  events,  while  two  events  (No.  1  and  No.  7) 
were  not  picked.  The  lower  traces  are  from  days  214  and  217.  Day  214  is  quiet  with  no  detec¬ 
tions,  while  day  217  show  two  detections  resulting  from  “interference”  by  large  earthquakes 
(located  far  from  Lop  Nor).  The  continuous,  smooth  line  above  each  threshold  trace  is  an  esti¬ 
mated  “trigger  level”  for  identifying  significant  peaks  on  the  threshold  trace. 

6.2.4  Location  Calibration 

Oslo  Workshop  on  location  calibration 

A  workshop  was  held  in  Oslo,  Norway  during  22-26  April  2002  in  support  of  the  global  seis¬ 
mic  event  location  calibration  effort  currently  being  undertaken  by  Working  Group  B  of  the 
CTBTO  Preparatory  Commission  in  Vienna,  Austria.  Dr.  Frode  Ringdal  chaired  the  meeting, 
which  was  attended  by  66  experts  from  12  countries  and  the  Provisional  Technical  Secretariat 
of  the  CTBTO  The  recommendations  from  this  workshop  have  been  provided  in  the  paper 
CTBTAVGB/TL-2/70,  issued  by  Working  Group  B  of  the  CTBTO  Preparatory  Commission. 
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6.2.5  Conclusions  and  Recommendations 

The  analysis  of  the  23  February  2002  Novaya  Zemlya  seismic  event  has  reconfirmed  our  previ¬ 
ous  estimates  of  the  detection  and  location  capability  of  the  regional  network  in  northern  Eu¬ 
rope.  This  experimental  network  operation  will  continue,  and  will  be  used  to  evaluate  and 
supplement  the  IMS  network  in  this  region.  Further  research  is  required  in  regional  event  iden¬ 
tification  criteria,  especially  depth  estimation.  In  the  case  of  depth  estimates,  a  detailed  study  of 
the  Vp/Vs  ratio  as  well  as  its  uncertainty  limits  should  be  undertaken. 

Our  initial  performance  tests  of  Site-Specific  Threshold  Monitoring  for  the  Lop  Nor  test  site 
were  successful.  We  recommend  that  these  initial  tests  be  followed  by  more  detailed  studies, 
where  in  particular  the  calibration  parameters  should  be  more  firmly  established.  It  is  also  of 
interest  to  include  additional  sensitive  IMS  arrays  in  the  monitoring  network,  and  evaluate  the 
ensuing  improvement  in  capability. 

The  location  calibration  effort  will  continue  to  be  an  important  part  of  our  work.  The  recommen¬ 
dations  provided  in  the  paper  CTBTAVGBArL-2/70  should  be  followed  up  by  the  international 
community,  and  the  progress  of  this  work  will  be  reviewed  in  a  planned  workshop  in  Oslo  during 
4-9  May  2003.  ^ 

F.  Ringdal 
T.  Kvaerna 
E.  Kremenetskaya 
V.  Asming 
S.  MykkeItveit 
C.  Lindholm 
J.  Schweitzer 
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6.3  Design  Study  for  the  Refurbishment  of  the  SPITS  Array  (AS  72) 

6.3.1  Introduction 

The  SPITS  array  is  located  on  Jansonhaugen  in  Adventdalen  on  Spitsbergen,  Svalbard,  approx¬ 
imately  15  km  east-southeast  of  Longyearbyen.  Jansonhaugen  is  a  hill  in  the  middle  of  a  valley 
(Adventdalen),  and  the  largest  part  of  the  array  is  deployed  on  the  plateau  of  this  hill.  The  rocks 
at  the  site  are  of  Cretaceous  age,  partly  covered  by  thin  top  soil  or  by  moraine  deposits  of  vari¬ 
able  depth.  The  twelve  sensors  of  the  array  are  placed  at  the  bottom  of  nine  6  m  deep  cased 
boreholes.  The  bottom  of  the  boreholes  are  either  in  Cretaceous  rock  or  in  moraine  material  in 
stable  permafrost  conditions  (temperature  approximately  -5°C  all  year  around  at  a  depth  of  6 
m),  such  that  there  is  no  melting/freezing  taking  place  at  this  depth.  The  9  array  sites  are 
deployed  in  two  rings  around  a  center  element  (see  Fig.  6.3.1).  The  radii  of  the  rings  are 
approximately  250  m  for  the  A-ring  and  500  m  for  the  B-ring.  More  detailed  descriptions  of  the 
array  and  an  initial  data  processing  were  presented  by  Mykkeltveit  et  al.  (1992)  and  Fyen  and 
Ringdal  (1993). 

As  of  summer  2002  the  SPITS  array  consists  of  the  following  equipment: 

•  9  vertical  component  short-period  seismometers,  Guralp  CMG-3V-borehole,  placed  at 
the  bottom  of  6  m  deep  cased  boreholes. 

•  1  three  component  broad-band  seismometer  of  type  Guralp  CMG-3T-borehole,  also 
placed  at  the  bottom  of  the  borehole  at  array  site  SPB4.  This  instrument  has,  however, 
been  out  of  operation  since  March  2001. 

•  2  Nanometrics  RD-6  digitizers,  16-bit  resolution,  40  Hz  sampling  rate. 

The  transfer  function  of  the  SPITS  array  is  shown  in  Fig.  6.3.2.  The  yellow  dashed  circle  indi¬ 
cates  the  highest  wavenumber  that  is  likely  to  be  considered  for  any  signals  at  SPITS,  i.e.,  high 
frequency  Rg-phases  (apparent  velocity:  1.3  km/s,  dominant  frequency:  6  Hz).  As  seen  from 
Fig.  6.3.2,  the  array  exhibits  significant  sidelobes  that  are  likely  to  influence  the  data  processing 
of  Rg-phases  due  to  aliasing. 

6.3.2  Problematic  Signals 

We  will  in  this  section  show  some  typical  examples  of  signals  that  are  difficult  to  process  prop¬ 
erly  with  the  current  SPITS  configuration.  The  problems  are  related  to  both  signal  detection, 
and  the  estimation  of  the  slowness  vector  (i.e.,  apparent  velocity  and  backazimuth)  of  the 
detected  phases. 

The  Spitsbergen  array  (SPITS)  often  reports  a  large  number  of  detections  (Schweitzer,  1998), 
which  can  reach  thousands  per  day.  A  detailed  analysis  shows  that  the  majority  of  these  detec¬ 
tions  are  real  seismic  signals  caused  by  small  sources  located  at  close  distances.  The  numerous 
local  events  (see  e.g..  Fig.  6.3.3)  typically  show  P  onsets,  no  weU  defined  S  onsets,  and  domi¬ 
nant  Rg  onsets.  Investigating  the  mean  backazimuths  of  these  Rg  phases  one  can  see  a  strong 
correlation  with  the  direction  to  the  coal-mine  area  at  about  8  km  distance  and  directions  to 
nearby  glaciers  (shortest  distance  about  3  km),  for  more  details  see  Schweitzer  (1998). 

At  the  SPITS  array  we  also  observe  a  relatively  large  number  of  high  SNR  P  detections  from 
regional  events  without  any  corresponding  detection  of  an  S  phase.  Fig.  6.3.4  shows  seismo¬ 
grams  observed  at  the  three-component  site  SPB4  from  an  earthquake  located  at  the  Kmpovich 
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Ridge  west  of  Svalbard  at  an  epicentral  distance  of  about  3°.  In  addition  to  the  original  horizon¬ 
tal  components,  the  rotated  radial  and  transverse  components  are  shown  (at  the  bottom).  The  S 
phase  has  a  low  SNR  and  the  signal  is  particularly  weak  in  SV  energy,  which  should  be  visible 
on  the  radial  and  vertical  components.  The  S  phase  can,  however,  be  observed  on  the  transverse 
component  SPB4_bt. 

Due  to  the  presence  of  low-velocity  materials  in  the  uppermost  crust,  Sn  arrivals  from  regional 
events  at  larger  distances  also  show  the  largest  amplitudes  and  SNR  on  the  radial  or  transverse 
components.  As  an  example,  we  show  in  Fig.  6.3.5  SPITS  recordings  of  the  Kara  Sea  event  on 
16  August  1997.  The  epicentral  distance  to  this  event  was  1 1.6°.  The  upper  trace  shows  the  ver¬ 
tical-component  array  beam  steered  with  the  S-phase  apparent  velocity  of  4.0  km/s.  The  three 
lower  traces  show  the  vertical,  radial  and  transverse  components  for  the  three-component  sen¬ 
sor  SPB4.  All  traces  of  Fig.  6.3.5  are  bandpass  filtered  between  4  and  8  Hz.  Notice  the 
improved  SNR  on  the  radial  and  transverse  components  as  compared  to  the  vertical-component 
S-beam. 

633  Spatial  aliasing  of  high  frequency  Rg-phases 

As  seen  from  Fig.  6.3.2,  the  transfer  function  of  the  SPITS  array  exhibits  significant  side  lobes 
that  are  likely  to  influence  the  data  processing  of  the  numerous  Rg-phases,  all  having  low 
apparent  velocities.  The  low  apparent  velocities  of  the  Rg  phases  are  caused  by  the  low  seismic 
velocities  of  the  Cretaceous  sediments. 

Schweitzer  (1998)  has  introduced  a  rather  complex  signal  processing  scheme  to  reduce  abasing 
problems  when  processing  the  Rg  onsets.  There  are,  however,  still  cases  when  Rg  phases  are 
mis-identified  as  P-onsets,  and  vice  versa.  In  such  cases  the  F-K  analysis  picks  the  energy  max¬ 
imum  on  a  sidelobe,  in  particular  for  onsets  with  low  SNR. 

The  numerous  local  events  often  contain  a  mixture  of  Pg,  Sg,  and  Rg  within  a  time  window  of 
2-5  seconds  (see  Fig.  6.3.3).  In  order  to  estimate  the  slowness  and  backazimuth  of  such  onsets, 
short  time  windows  are  needed  for  the  F-K  analysis.  With  the  current  SPITS  configuration  we 
will  need  at  least  a  2  s  long  time  window  for  processing  of  the  low-velocity  Rg-phases.  For 
these  local  events  an  automatically  positioned  2  s  long  time  window  often  contains  all  types  of 
phases  (Pg,  Sg,  Rg). 

In  order  to  remove  the  side  lobes  of  the  array  transfer  function  and  to  give  us  the  possibihty  to 
analyze  the  low-velocity  onsets  in  shorter  time  windows,  we  propose  to  install  three  additional 
sensors  close  to  the  center  element  SPAO.  The  geometry  of  an  array  cannot  be  inverted  from  a 
given  transfer  function.  It  must  be  found  by  a  trial-and-error  search  trough  different  geometries. 
Therefore,  to  find  a  position  for  the  three  new  sensors,  which  gives  the  best  suppression  of  the 
side  lobes  in  Fig.  6.3.2,  we  applied  a  systematic  search  by  adding  one  station  after  the  other  at 
different  locations  and  calculating  the  array-transfer  function  (see  e.g.,  Schweitzer  et  al.,  2002). 
The  searched  area  was  a  square  of  600*600  m  centered  at  the  central  site  SPAO.  In  this  area  aU 
possible  site  locations  on  a  5*5  m  grid  were  tested  as  locations  for  the  new  array  sensors;  the 
criterion  for  optimizing  the  transfer  function  was  to  minimize  this  function  in  the  slowness 
range  between  the  main  lobe  and  the  yellow  broken  line  of  Fig.  6.3.2.  The  best  solution  found 
is  shown  in  Fig.  6.3.6  and  the  location  of  the  three  newly  proposed  sensors  SPEl,  SPE2,  and 
SPE3  to  reach  this  solution  is  seen  on  Fig.  6.3.7.  Notice  the  reduction  of  the  side  lobes  as  com¬ 
pared  to  Fig.  6.3.2. 
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63.4  Low  resolution  of  the  F-K  estimates 

As  seen  from  Fig.  6.3.6,  the  transfer  function  of  the  SPITS  array,  including  the  three  new  sites 
in  the  center,  has  a  relatively  wide  main  lobe.  This  feature  hmits  the  resolution  capabihty  of  the 
array,  which  is  critical  for  classifying  signals  as  local,  regional  or  teleseismic  based  on  the 
apparent  velocity  estimates  from  the  F-K  analysis.  In  addition,  the  backazimuth  estimates  from 
the  F-K  analysis  have  large  uncertainties. 

In  order  to  narrow  the  main  lobe  we  have  to  extend  the  aperture  of  the  array.  Based  on  experi¬ 
ence  with  other  regional  arrays  in  Fennoscandia  (FINES,  ARCES  and  NORES)  we  propose  to 
install  new  instruments  in  a  C-ring  hke  configuration.  The  plateau  on  the  hill  on  which  SPITS 
is  located  is  surrounded  by  steep  slopes  and  rivers  both  to  the  north  and  south,  and  the  array  can 
therefore  not  be  extended  in  these  directions.  From  analysis  of  topography  maps,  sateUite  pho¬ 
tos,  and  from  inspections  during  visits  to  the  SPITS  array  we  propose  a  new  configuration  as 
shown  in  Fig.  6.3.7:  three  C-ring  sites  towards  east  (SPCl,  SPC2,  SPC3)  and  three  C-ring  sites 
towards  west-southwest  (SPC4,  SPC5,  SPC6).  The  exact  positions  of  these  additional  six  new 
C-ring  sites  were  also  found  by  a  systematic  trial-and-error  search  on  a  10*10  m  grid.  The 
array-transfer  function,  for  which  the  main  lobe  became  the  smallest,  was  chosen  as  the  best 
one. 

The  transfer  function  of  the  new  configuration  of  the  SPITS  array,  i.e.  the  original  nine  sites 
extended  by  three  sites  in  the  center  (SPEl,  SPE2,  SPE3)  and  six  sites  in  the  C-ring  (SPCl, 
SPC2,  SPC3,SPC4,  SPC5,  SPC6)  is  shown  in  Fig.  6.3.8.  Notice  that  the  main  lobe  has  become 
narrower  as  compared  to  Figs.  6.6  and  6.2,  and  that  it  is  elongated  in  the  north-south  direction 
due  to  the  non-circular  array  geometry.  The  relative  coordinates  of  all  sites  with  respect  to  the 
center  site  SPAO  are  listed  in  6.4.1. 

6.3.5  Detection  and  F-K  analysis  of  S-phases 

As  shown  in  Figs.  6.4  and  6.5,  regional  S-phases  at  SPITS  in  general  have  the  largest  amph- 
tudes  and  SNR  on  the  horizontal  components.  The  current  SPITS  configuration  has  only  one 
three-component  instrument  at  site  SPB4.  This  instrument  has,  however,  been  quite  unstable  in 
operation.  In  order  to  detect  S-phases  more  efficiently,  more  three  component  instruments  are 
needed.  While  one  three-component  instrument  already  improves  the  detection  capability  of  S- 
phases  (see  Figs.  6.4  and  6.5),  the  subsequent  F-K  analysis  would  suffer  from  low  SNR  as  the 
vertical  component  sensors  have  to  be  used  for  this  purpose. 

With  more  three-component  sites,  the  detection  capahihty  of  S-phases  can  be  further  improved 
by  creating  S-velocity  beams  from  the  rotated  radial  and  transverse  components,  as  now  done 
for  ARCES  and  NORES  (Kvsma  et  al.,  1999).  In  addition,  the  F-K  analysis  can  now  be  run  on 
the  horizontal  components,  as  e.g.,  demonstrated  by  Schweitzer  (1994). 

As  a  first  attempt,  we  calculate  the  array  response  of  4  three-component  sensors,  which  is  a 
minimum  number  of  traces  for  running  an  F-K  analysis.  After  testing  several  combinations  of 
the  SPITS  geometry,  our  best  result  was  obtained  by  using  SPAl,  SPC3,  SPC4  and  the  already 
existing  three-component  sensor  SPB4  (see  Fig.  6.3.7).  With  only  four  sensors  we  were  unable 
to  remove  aU  side  lobes  from  the  S-phase  wavenumber  range,  but  we  can  obtain  a  relatively 
narrow  main  lobe.  However,  with  four  sensors  only,  we  will  lose  the  capability  to  run  an  F-K 
analysis  on  the  horizontal  components  in  the  case  that  one  of  the  horizontal  sensors  is  down. 
The  best  “reahstic”  three-component  configuration  was  obtained  by  having  three-component 
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instruments  at  all  five  B-ring  sites  and  at  SPAO  in  the  center.  Such  a  configuration  can  improve 
the  SNR  for  S  onsets  by  a  factor  of  about  2.5;  the  corresponding  transfer  function  of  Fig.  6.3.9 
shows  no  side  lobes  within  the  wavenumber  range  for  regional  Sn  phases  (i.e.,  between  the  red 
and  the  yellow  rings). 

6  J.6  Sampling  rate,  high-frequency  spectra,  and  24  bits  digitizers 

With  the  current  sampling  rate  of  40  Hz  at  SPITS,  we  can  conduct  array  processing  for  signal 
frequencies  up  to  16-20  Hz.  For  higher  frequencies  the  coherency  among  the  signals  at  the  dif¬ 
ferent  sensors  is  heavily  reduced  such  that  beamforming  and  F-K  analysis  have  little  effect. 
Therefore,  we  see  no  need  to  increase  the  general  sampling  rate  for  SPITS. 

However,  it  is  well  known  that  we  at  SPITS  observe  high  signal  frequencies  (e.g..  Bowers  et 
al.,  2001).  Therefore  we  suggest  that  one  of  the  three-component  sensors  (e.g.,  SPAO)  is  sam¬ 
pled  at  120  Hz,  or  alternatively  80  Hz. 

In  order  to  meet  the  standard  of  the  CTBT  International  Monitoring  System  an  upgrade  of  the 
digitizers  to  24  bit  words  is  needed. 

63.7  Conclusion 

In  conclusion,  the  following  extensions  are  suggested  for  the  SPITS  array  to  increase  its  detec¬ 
tion  and  monitoring  capabilities: 

•  Three  new  sites  in  the  center  of  the  array  (SPEl,  SPE2,  SPE3)  equipped  with  vertical 
component  instruments. 

•  Six  new  sites  in  a  C-ring  equipped  with  vertical  component  instruments. 

•  Add  5  three-component  instruments  to  SPAl,  SPBl,  SPB2,  SPB3  and  SPB5.  The 
existing  broad-band  three-component  instrument  at  SPB4  can  be  used  for  array  pro¬ 
cessing  after  digitally  reshaping  the  response. 

•  Install  24  bits  digitizer  for  all  channels. 

•  Sample  one  of  the  three-component  sensors  at  120  Hz,  e.g.,  SPAO. 


Johannes  Schweitzer 
Tormod  Kvaema 
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Table  6J.L  The  relative  coordinates  of  the  old  and  nevt'ly  proposed  sites  of  the  SPITS  array  with 
respect  to  the  central  element  SPAO,  which  is  located  at  78.1777°  N  and  16J6998°  E. 


Site 

XCE-W) 

[m] 

Y(N-S) 

[m] 

SPAO 

0 

0 

SPAl 

125 

227 

SPA2 

150 

-202 

SPA3 

-255 

-44 

SPBl 

470 

228 

SPB2 

333 

-381 

SPB3 

-255 

-451 

SPB4 

-501 

123 

SPB5 

-49 

513 

SPEl 

150 

-30 

SPE2 

-30 

-160 

SPE3 

-230 

-240 

SPCl 

1030 

500 

SPC2 

1100 

-200 

SPC3 

745 

-240 

SPC4 

-800 

-835 

SPC5 

-710 

-460 

SPC6 

-1140 

-40 
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2002-102:12 

Site  X 

SPLA0-S2 

SPLAI^z 

SPLA2.S2 

SPLA3.S2 

SPLB1-S2 

SPLB2-SZ 

SPLB3.S2 

SPI.B4.sz 

SPL84.bz 

SPI.B4.bn 

SPLB4.be 

SPLB5.SZ 
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-501 

-501 

-501 

-49 


218 

-381 

-451 


Z  height 
323 
320 
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200 

234 

340 

340 

340 

340 

295 


Fig.  6.3.1.  Current  geometry  of  the  SPITS  array. 


Fig.  6.3.2.  The  transfer  function  of  the  current  SPITS  array  ( vertical  components). 
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Fig.  6.3.3.  Example  of  a  local  event  observed  at  SPITS,  located  in  the  glacier  Gl0ttjjellhreen 

(hackazimuth  =  141°,  D  =  4  km).  The  traces  are  the  hand  pass  filtered  (3-8  Hz)  vertical-com¬ 
ponent  seismograms. 
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Fig.  6.3.4.  Filtered  seismograms  (3-6  Hz)  at  the  three-component  site  SPB4from  an  earthquake 
located  at  tlw  Knipovich  Ridge  northwest  of  Svalbard  at  an  epicentral  distance  of  about  3°. 
The  rotated  radial  and  transverse  components  are  shown  (at  the  bottom). 
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Fig.  6.3.5.  SPITS  observations  of  the  Kara  Sea  event  on  16.  August  1997.  The  upper  trace  shows 
the  vertical-component  array  beam  steered  with  the  S-phase  apparent  velocity  of  4.0  km/s. 
The  three  lower  traces  show  the  vertical,  radial  and  transverse  components  for  the  three- 
component  sensor  SPB4.  All  traces  are  bandpass  filtered  between  4  and  8  Hz. 
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Fig.  6.3.6.  The  transfer  function  of  the  SPITS  array  extended  hy  three  sites  (SPEI,  SPE2.  SPE3)  in 
the  center,  note  the  side-lohe  reduction  inside  the  yellow  ring  with  respect  to  Fig.  6.3.2. 
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Fig.  6.3.7.  The  suggested  extended  configuration  of  the  SPITS  array  plotted  on  a  satellite  photo  of 
Jansonhaugen.  The  original  configuration  is  extended  by  three  sites  (El,  E2,  E3)  in  the  cen¬ 
ter  and  six  sites  in  the  C-ring  (Cl,  C2,  C3,  C4,  C5,  C6).  Green  symbols:  sites  unchanged 
from  the  old  configuration.  White  circles:  new  sites  equipped  with  vertical  component  sen¬ 
sors.  Green  triangles:  Sites  of  the  old  configuration  upgraded  to  three-component  sensors. 
High-frequency  sensor:  the  three-component  instrument  at  AO  should  be  sampled  at  e.g.,  120 
Hz. 
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Fig.  6.3.8.  The  transfer  function  of  the  new  SPITS  array  configuration  for  vertical  sensors  at  all 
sites  as  shown  on  satellite  picture  in  Fig.  6.3.7. 
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Fig.  6.3.9.  The  transfer  function  of  the  6  SPITS  array  sites  SPAO,  SPBl,  SPB2,  SPB3,  SPB4  and 
SPB5.  The  yellow  and  red  circles  represent  the  wavenumber  range  for  regional  Sn  phases. 
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6.4  Analysis  of  cavity-decoupled  chemical  explosions 

This  research  is  conducted  under  Defense  Threat  Reduction  Agency  contract  no. 
Ol-OI-C-0069 

Introduction 

Cavity  decoupling  has  long  been  acknowledged  to  be  the  most  effective  means  of  evading 
detection  of  an  underground  nuclear  test,  and  is  thus  one  of  the  most  probable  scenarios  for  a 
violation  of  the  Comprehensive  Nuclear  Test  Ban  Treaty  (CTBT).  NORSAR  has  recorded  and 
analysed  seismic  data  from  a  series  of  underground  cavity  explosions  in  the  Kopparberg  region 
of  Sweden,  with  the  aim  of  quantitatively  predicting  the  decoupling  effect  of  underground  cav¬ 
ities  of  known  dimensions  and  shape. 

A  database  of  seismic  waveform  data  has  been  compiled  for  seven  such  cavity  explosions.  A 
series  of  four  explosions  has  taken  place  between  December  2000  and  June  2002  using  varying 
quantities  and  compositions  of  explosives  within  a  chamber  of  size  lOOOm"*  with  an  overburden 
of  rock  of  approximately  100m.  In  addition,  NORSAR  also  recorded  a  series  of  explosions 
between  1987  and  1989  which  took  place  within  smaller  chambers  on  the  same  site.  Only  three 
factors  differentiate  these  events;  the  quantity  of  explosive,  the  composition  and  arrangement 
of  the  explosives,  and  the  size  and  configuration  of  the  explosion  chamber. 
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Fig.  6.4.1.  The  location  of  the  explosion  sites  relative  to  the  HFS  and  NRS  arrays  and  the  NORSAR 
sub-arrays.  The  red  line  is  the  Norway-Sweden  national  boundary. 

The  cavity  explosion  data  has  been  complemented  by  records  from  a  series  of  surface  explo¬ 
sions  at  Mossibranden,  a  site  approximately  15  km  from  the  underground  cavities.  The  distance 
between  the  surface  and  cavity  explosions  is  small  compared  with  the  distance  between  the 
explosion  sites  and  the  NORSAR,  NORES  and  Hagfors  seismic  arrays  (see  Fig.  6.4.1).  This 
will  largely  eliminate  path  effects  and  allow  a  comparison  of  source  properties. 
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Cavity  decoupled  explosions 

The  cavity  explosions,  which  comprise  two  separate  series  of  experiments,  are  listed  in  Table 
6.4.1.  The  first  six  of  these  were  recorded  by  the  NORES  array,  which  was  unfortunately  put 
out  of  action  by  lightning  shortly  before  the  June  2002  event;  waveforms  from  the  NRAO 
instrument  are  displayed  in  Fig.  6.4.2.  The  signals  resulting  from  the  three  events  from  1987 
and  1989  have  far  higher  amplitudes  than  the  later  events,  which  probably  reflects  the  differ¬ 
ences  in  charge/cavity  volume  as  displayed  in  Table  6.4.1.  The  signal  from  event  2001  Cl 50 
(2500kg  TNT  in  the  lOOOm^  chamber)  is  indistinguishable  from  the  noise  without  filtering  of 
the  data.  The  explosions  were  set  off  in  chambers  open  to  the  access  tunnels  (unsealed  explo¬ 
sions). 


Table  6.4.1.  Cavity  decoupled  explosions  at  the  Alvdalen  site.  Origin  times  of  the  1987-1989  events 
were  estimated  from  arrival  times  at  NORES.  Origin  times  of  the  2000-2002 
events  were  determined  from  a  station  at  the  explosion  site. 


Origin  ID 

Explosion  origin  time 

Explosion 

charge 

(kg) 

Explosive 

Chamber 

volume 

(in^)“ 

Charge/ 

Chamber 

volume 

(kg/m^) 

1987C146 

1987-146:10.47.38.2 

5000 

ANFO^’ 

300 

16.7 

1987C259 

1987-259:10.36.13.0 

5000 

ANFO 

200 

25.0 

1989C263 

1989-263:10.06.03.5 

5000 

ANFO 

300 

16.7 

2000C348 

2000-348:10.03.02.0 

10000 

TNT 

1000 

10.0 

2001C150 

2001-150:10.03.56.2 

2500 

TNT 

1000 

2.5 

2001C186 

2001-186:10.41.23.5 

10000 

TNT  and 
ammuni¬ 
tion  shells 

1000 

10.0 

2002C164 

2002-164:08.59.25.1 

10000 

TNT/ 

powder 

1000 

10.0 

a.  Not  including  access  tunnel. 

b.  The  ANFO  explosive  used  in  the  1987  and  1989  events  had  an  explosion  equivalent  of 
0.82  of  TNT 


In  order  to  quantify  the  differences  between  these  six  events,  we  examine  the  power  density 
spectra  (PDS).  For  each  component,  three  time  windows  were  defined  corresponding  to  the  P 
arrival,  the  S  arrival  and  pre-event  noise.  The  P-window  was  defined  beginning  at  the  P-wave 
onset  and  ending  5  seconds  later.  The  S-window  consists  of  the  10  seconds  immediately  fol¬ 
lowing  the  S-wave  onset.  A  third  10  second  window,  ending  2  seconds  before  the  P-wave 
onset,  was  defined  to  measure  the  background  noise  level. 

For  each  event,  the  PDS  was  calculated  for  the  three  time  windows  for  each  vertical  component 
trace.  The  PDS,  averaged  over  all  vertical  component  traces  for  the  NORES  array,  is  shown  for 
the  first  six  decoupled  explosions  in  Fig.  6.4.3.  In  order  to  provide  a  common  basis  for  compar- 


79 


NORSAR  ScL  Rep.  2-2002 


August  2002 


ison,  we  have  converted  all  spectra  to  the  characteristics  of  the  HFS  response  which,  unlike  the 
NORES  array  response,  is  flat  to  velocity  over  the  bandpass. 


19 S 7— 146:1  0.47, 4.?, S21  19B7C146  i 

Yield  5000  kg.  Chamber  300  ra*3.  | 


2001  -  1  5>0-  1  0.04.00.826  200  1  C  1  50 

Yield  2500  kg.  Chamber  1000  m*3. 


200 1-1861041,28137  200lClS6i  1 

Yield  10000  kg.  Chamber  1000  m''3. 


Fig.  6.4.2.  Unfiltered  waveforms  for  the  6  decoupled  explosions  prior  to  June  2002  from  the  3-com¬ 
ponent  NRAO  instrument  of  the  NORES  array.  North  and  East  waveforms  have  been  rotated 
into  radial  and  transverse  components.  Maximum  values  of  traces  (counts)  are  indicated  in 
red. 

The  results  shown  in  Fig.  6.4.3  clearly  show  the  effect  of  substantially  greater  decoupling  for 
the  larger  chamber.  The  spectra  for  events  1987C146  and  1987C263  show  great  similarity, 
which  is  to  be  expected  given  that  both  occured  in  the  same  chamber  using  the  same  quantiw 
and  type  of  explosives.  Event  1987C259,  which  occured  in  the  200m^  rather  than  the  300m^ 
chamber,  has  a  spectrum  which  is  similar  to  the  1987C146  and  1987C263  spectra  at  high  fre¬ 
quencies  but  which  differs  significantly  in  the  3  to  8  Hz  frequency  band.  Events  2000C348  and 
2001  Cl  86  both  took  place  within  the  lOOOm^  chamber  using  the  same  yield  of  explosives. 
However,  the  PDS  for  event  2000C348  is  greater  than  that  for  200 1C  186  by  an  approximate 
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factor  of  2  for  all  frequencies.  The  only  difference  between  these  two  events  was  that  pure  TNT 
was  used  as  the  explosive  for  2000C348  whereas  TNT  mixed  with  ammunition  shells  was  used 
for  event  2001  Cl 86. 


Fig.  6.4.3.  Power  density  spectra,  averaged  over  all  vertical  component  traces  from  the  MORES 
array  for  six  cavity  explosions  as  indicated.  Solid  red  lines  indicate  PDS  for  P -window  (5 
seconds  following  P),  dashed  blue  lines  indicate  PDS  for  the  S-window  (10  seconds  follow¬ 
ing  S)  and  the  dotted  black  line  indicates  PDS  for  the  noise  window  (12  to  2  seconds  before 
P).  At  frequencies  for  which  the  PDS  for  P  or  S  window  is  less  than  4.0  times  that  for  the 
noise  window,  the  signal  spectrum  is  shown  as  a  thin  line. 

Using  the  average  of  vertieal  component  traees  from  the  NRS,  NOA  and  HFS  arrays,  we  com¬ 
pare  the  power  density  spectra  from  the  P-arrival  windows  from  all  of  the  seven  cavity  explo¬ 
sions  listed  in  6.4.1.  Performing  a  simple  arithmetic  mean  on  these  spectra  recorded  at  different 
locations  is  justified  by  the  comparable  source-receiver  distances  between  the  explosion  sites 
and  seismic  arrays  (see  Fig.  6.4.1).  Figure  6.4.4  demonstrates  the  greater  PDS  in  the  signal 
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resulting  from  explosions  in  the  smaller  chambers.  The  10000  kg  of  explosives  used  for  the 
2002C1 64  event  consisted  of  6600  kg  powder  and  3400  kg  TNT  in  shells. 

All  of  the  spectra  in  Fig.  6.4.4  exhibit  low  power  at  low  frequencies,  although  the  ratio  of  low- 
frequency  power  density  to  high-frequency  power  density  is  lower  for  the  events  within  the 
1000m  chamber,  indicating  a  higher  degree  of  decoupling.  The  signal  to  noise  ratio  for  event 
2001  Cl 50  was  very  low  for  frequencies  below  10  Hz  and  a  low-frequency  power  to  high-fre¬ 
quency  power  ratio  cannot  be  estimated. 


KEY 

1987C146  — 

1987C259  — 

1989CPG3 

2001C150  — 

2001C186  -~ 

2002C164  — 


Fig.  6.4.4.  Power  density  spectra  for  the  P  time  window  averaged  over  all  vertical  component 

instruments  for  the  arrays  NOA,  NRS  and  HFS.  Thin  lines  indicate  that  the  ratio  PDS(sig- 

0  for  a  given  frequency.  Events  J987CJ46,  1987C259and 
1989C263  only  include  instruments  from  the  NORES  array.  Event  2002CI64  only  includes 
instruments  from  the  HFS  and  NOA  arrays. 


A  comparison  between  the  underground  cavity  explosions  and  surface  explosions  at 
Mossibrdnden 

The  Swedish  Armed  Forces  regularly  detonate  outdated  ammunition  at  ground  level  at  a  site 
approximately  15  km  from  the  underground  chambers  (see  Fig.  6.4.1).  Seismic  waveforms 
from  12  such  events  were  added  to  the  database  to  provide  a  comparison  between  underground 
explosions  (where  the  explosive  is  completely  surrounded  by  air)  and  detonations  at  the  surface 
in  the  open  air.  The  surface  explosions  at  the  Mossibranden  site  are  listed  in  Table  6.4.2. 
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Table  6.4.2.  Explosions  at  the  Mossibranden  site  in  Alvdalen,  June  and  July  2001. 

Charges  are  TNT  weights. 


Origin  ID 

Date 

Charge/ 

kg 

Origin  time  as  determined  by 
NDC,  NORSAR 

2001S176 

25  June  2001 

18831 

2001-176:13.46.17.7 

200 IS 177a 

26  June  2001 

20327 

2001-177:07.15.31.2 

2001S177b 

26  June  2001 

20408 

2001-177:13.00.11.1 

2001 S 178a 

27  June  2001 

20101 

2001-178:09.16.04.5 

200 IS 178b 

27  June  2001 

20432 

2001-178:13.40.12.4 

2001S179a 

28  June  2001 

20793 

2001-179:09.40.55.4 

2001 S 179b 

28  June  2001 

20242 

2001-179:13.50.31.6 

2001S183 

2  July  2001 

21024 

2001-183:11.36.00.6 

2001S184a 

3  July  2001 

19044 

2001-184:07.31.01.9 

2001S184b 

3  July  2001 

20540 

2001-184:13.01.00.0 

2001S185 

4  July  2001 

19118 

2001-185:09.36.05.6 

2001S186 

5  July  2001 

20102 

2001-186:08.46.05.1 

In  Fig.  6.4.5  we  compare  the  waveforms  resulting  from  one  of  the  underground  cavity  explo¬ 
sions  (2002C164)  and  one  of  the  surface  explosions  (2001  SI 86)  filtered  through  several  differ¬ 
ent  frequency  bands.  The  amplitudes  resulting  from  the  surface  explosion  are  somewhat  greater 
than  those  resulting  from  the  decoupled  explosion.  This  is  particularly  pronounced  for  low  fre¬ 
quencies;  in  the  1.5-3.0  Hz  frequency  band,  the  signal  from  2002C164  barely  exceeds  the 
background  noise  level  whereas  there  is  a  clear  signal  at  these  low  frequencies  from  the  surface 
explosion.  In  the  3-6  Hz  band  there  is  a  weak  S-wave  signal  from  the  cavity  explosion  but  no 
discernible  P-wave  arrival.  In  the  higher  frequency  bands,  the  cavity  explosion  has  a  strong  P- 
signal  which  is  comparable  to  that  from  the  surface  explosion. 

The  spectral  properties  of  the  energy  resulting  from  these  two  events  is  seen  more  clearly  in  the 
power  density  spectra;  see  Fig.  6.4.6.  The  clear  distinction  between  the  two  events  is  the  energy 
in  the  P  spectrum  at  low  frequencies.  The  underground  cavity  explosion  on  the  left  has  a  peak 
of  P-energy  between  8  Hz  and  16  Hz,  with  approximately  an  order  of  magnitude  less  energy 
below  8  Hz.  By  contrast,  the  surface  explosion  on  the  right  has  slightly  more  P-wave  energy 
between  3  Hz  and  8  Hz  than  at  higher  frequencies.  There  is  less  difference  between  the  S-wave 
spectra  for  the  two  events. 
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Fig.  6.4.5.  Waveform  data  from  the  3-component  HFSC2  instrument  for  events  2002C164  (cavity 
explo.sion:  10000kg  powder  and  TNT  in  shells  in  a  lOOOm^  cavity)  and  2001 SI 86  (surface 
explosion  of 201 02kg  out-of-date  ammunition).  Vertical  component  traces  are  in  black; 
rotated  radial  and  transverse  components  are  in  red  and  blue  respectively.  All  traces  are 
scaled  identically  with  the  largest  amplitude  as  indicated.  Frequency  hands  for  filterins  are 
shown  on  the  right  hand  side.  '  j  j 
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Event  2002C164  Event  2001S186 


Fig.  6.4.6.  Power  density  spectra  (PDS)  of  the  radial  (top),  transverse  (center)  and  vertical  (bot¬ 
tom)  components  of  the  signals  from  events  2002C164  (decoupled,  left)  and  2001S186 
(Mossibranden,  right)  from  the  HFSC2  instrument.  The  solid  red  line  is  the  PDS  for  a  5  sec¬ 
ond  interval  starting  at  the  P  arrival,  the  dashed  blue  line  is  the  PDS  for  a  10  second  interval 
starting  at  the  S  arrival  and  the  dotted  line  is  the  PDS  for  a  10  second  inteiwal  ending  2  sec¬ 
onds  before  the  P  arrival  (noise  window).  The  dashed  and  solid  lines  are  displayed  with  dou¬ 
ble  thickness  for  frequencies  where  the  signal  to  noise  ratio  is  greater  than  4.0 . 

A  discriminant  for  cavity -decoupled  explosions 

The  comparison  between  events  2002C164  and  2001S186  (Fig.  6.4.5  and  Fig.  6.4.6)  indicates 
that  the  primary  difference  between  the  underground  cavity  explosions  and  the  surface  explo¬ 
sions  is  relatively  low  P  energy  at  low  frequencies  for  the  decoupled  explosions.  We  applied  the 
following  step-wise  procedure  to  find  a  quantity  which  summarises  the  visual  discriminant  in 
the  vertical  component  spectra  of  Fig.  6.4.6,  which  can  then  be  applied  to  all  of  the  events  in 
the  database  in  order  to  find  a  systematic  difference  between  the  cavity  (C)  events  and  the  sur¬ 
face  (S)  events. 
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For  each  event  in  the  data-base,  the  following  steps  were  carried  out! 

•  Waveform  data  was  obtained  for  every  available  vertical  component  station  in  the  HFS 
NRS  and  NOA  arrays. 

•  For  each  trace,  a  P-window  (5  seconds  following  the  P-arrival),  an  S-window  (10  seconds 
following  the  S-arrival)  and  a  noise  window,  N  (10  seconds  ending  2  seconds  before  the  P- 
wave  arrival),  were  defined. 

•  Power  density  spectra  were  calculated  over  all  time  windows  for  all  traces. 

•  Two  frequency  bands  were  defined:  band  1  (3  -  8  Hz)  and  band  2  (8  -  14  Hz).  Median  val¬ 
ues  of  the  P,  S  and  noise  power  density  spectra  were  calculated  for  frequencies  within 
these  bands.  For  a  given  trace,  we  denote  the  median  value  of  the  power  density  spectra 
between  frequencies  f  I  and  f2  P(fl-f2),  S(fl-f2)  and  N(fl-f2)  for  the  P  time  window,  S 
time  window  and  noi.se  time  window  respectively. 

•  For  each  trace,  the  ratios  P(8-14):N(8-14),  S(8-14):N(8-14),  P(3-8):N(3-8)  and  S(3- 
8):N(3-8)  were  calculated.  If  one  or  more  ofthe.se  ratios  was  below  2.0,  then  the  trace  was 
excluded.  For  the  remaining  traces,  the  ratios  P(8-14):P(3-8)  and  S(8-14):S(3-8)  were  eval¬ 
uated  and  the  mean  and  standard  deviations  calculated.  The  results  are  displayed  in  Fig. 

Fig.  6.4.7(a)  shows  a  clear  distinction  between  the  decoupled  explosions  (blue  bars)  and  the 
surface  detonations  (green  bars)  with  the  former  having  far  higher  ratios  of  P  energy  in  the 
higher  frequencies  to  P  energy  in  the  lower  frequencies.  This  ratio  is  highest  for  three  100(K)kg 
explosions  suggesting  that  they  were  the  mo.st  effectively  decoupled.  The  2500  kg  explosion 
(2001  Cl 50)  was  probably  even  more  effectively  decoupled,  but  was  excluded  from  this  analy¬ 
sis  on  the  grounds  of  poor  signal  to  noi.se  ratio.  This  ratio  is  approximately  an  order  of  magni¬ 
tude  greater  for  the  decoupled  explosions  than  for  the  surface  explosions. 

Fig.  6.4.7(b)  shows  a  similar  pattern  for  the  S  waves  but  with  a  much  smaller  difference  in  the 
ratios.  The  P(high):P(low)  ratio  appears  therefore  to  be  a  far  better  discriminant  than  the 
S(high):S(low)  ratio.  This  is  in  accordance  with  expectations,  since  only  P-energy  is  transmit¬ 
ted  through  the  air  surrounding  an  explosion  and  all  .shear  energy  is  generated  by  conversion 
from  P-energy.  Given  that  S  and  P  exhibit  similar  trends  (i.e.  decoupling  produces  both  less  P 
and  less  S  at  low  frequencies),  the  P(high):P(low)  ratio  is  probably  a  better  di.scriminant  than 
the  S:P  ratio. 

The  simplicity  of  this  single  ratio  discriminant  may  lead  to  spurious  conclusions  being  drawn 
and  it  therefore  provides  no  substitute  for  a  thorough  examination  of  the  spectra.  For  example, 
on  the  basis  of  the  P(8-14):P(3-8)  di.scriminant  (Fig.  6.4.7a),  event  1987C259  would  appear  to 
be  more  effectively  coupled  than  events  1987C146  and  1989C263:  two  events  which  involved 
equal  quantities  of  explosives  but  in  a  larger  chamber.  The  low  P(3-8)  value  for  event 
1987C259  is  the  result  of  a  decrea.se  in  power  between  5  and  8  Hz  (Fig.  6.4  4)  which  is  not 
shown  by  events  1987C146  and  1989C263.  Below  5  Hz,  event  1987C259  displays  more  P 
energy  than  the  other  two  events. 
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Event  identification  number 


Fig.  6.4.7.  Mean  values  of  (a)  log(  P(8-14)/P(3-8) )  and  (b)  log(  S(8-14)/S(3-8)  )for  all  vertical 

component  stations  in  the  three  arrays  HFS,  NRS  and  NOAfor  which  the  ratios  P :N  and  S:N 
were  greater  than  2.0  in  both  frequency  intervals.  Error  bars  surrounding  the  mean  values 
are  standard  deviations.  Decoupled  cavity  explosions  are  denoted  with  blue  bars  and  surface 
explosions  from  Mossibrdnden  with  green  bars. 
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Conclusions 

Throughout  the  analysis,  local  site  effects  have  been  reduced  as  a  result  of  averaging  of  values 
from  many  array  stations.  The  synthesis  below  and  the  main  conclusions  of  Fig.  6.4.7  are  based 
on  different  wavepaths  and  distances,  and  it  is  therefore  reasonable  to  assume  that  the  findings 
are  largely  source  related. 

•  The  recordings  from  explosions  in  the  cavities  at  Alvdalen  are  characterized  by  a  large 
amount  of  high-frequency  energy  (8-16  Hz)  compared  with  the  energy  at  lower  frequencies 
(2"5  Hz).  Such  spectral  characteristics  are  not  observed  for  any  of  the  other  events  included 
in  the  database.  The  effect  is  more  pronounced  for  P  than  for  S  waves. 

•  The  three  10000kg  explosions  in  a  lOOOm^  chamber  appear  to  have  been  more  efficiently 
decoupled  than  the  5000kg  explosions  in  chambers  of  200m^  and  300m^.  The  signals  from 
the  explosions  in  the  smaller  chambers  were  far  stronger,  in  spite  of  significantly  lower 
yields,  and  displayed  a  less  pronounced  high-to-low  frequency  energy  ratio  than  the  explo¬ 
sions  in  the  larger  chambers. 

•  The  2500kg  explosion  in  the  lOOOm^  chamber  appeared  to  be  very  effectively  decoupled, 
and  was  visible  on  our  recordings  only  at  frequencies  above  8  Hz. 

The  Office  of  Technology  Assessment  report  (OTA,  1 988)  estimates  that  a  one  kiloton  explo¬ 
sion  in  granite  at  a  depth  of  828  meters  would  be  decoupled  by  a  spherical  cavity  with  a  radius 
of  approximately  20  metres.  It  might  be  po.ssible  to  scale  this  to  a  yield  of  10  tons  at  a  depth  of 
100  metres,  as  in  Alvdalen,  to  obtain  the  cavity  volume  required  for  full  decoupling.  There  are, 
however,  large  uncertainties  associated  with  such  scaling,  in  particular  when  it  comes  to  the  use 
of  high-explosive  materials  and  non-spherical  cavities.  The  10  ton  explosions  in  Alvdalen  were 
detonated  in  a  cavity  5  meters  high  by  25  meters  long,  by  10  meters  wide  (-1000  m^),  and  the 
explosive  charge  was  made  up  of  10  simultaneous  one  ton  charges  detonated  on  the  floor  of  the 
cavity  (lifted  from  the  floor  on  styropor  ‘mattresses’).  It  is  therefore  uncertain  whether  the  10 
ton  shots  were  fully  or  only  partially  decoupled.  As  a  qualitative  indicator  for  the  degree  of 
decoupling  it  can  be  noted  that  the  damages  to  the  cavity  are  minor  and  limited  to  ‘surface 
scratches’,  and  the  repeated  use  of  the  same  cavity  is  warranted.  On  the  other  hand,  both  the 
200m'  and  300m  cavities  have  been  closed  for  further  use  due  to  structural  damage  caused  by 
the  explosions.  From  the  P-wave  spectra  shown  in  Figure  6.4.4  we  make  the  following  observa¬ 
tions: 

•  The  power  spectral  levels  of  the  2000C348  and  200 1 C 1 86  1 0  ton  explosions  detonated  in 
the  1000  m  cavity  differ  by  a  factor  2{j2  m  amplitude)  over  the  entire  observable  fre¬ 
quency  band.  At  high  frequencies,  the  2002C164  10  ton  explosion  has  a  similar  spectral 
level  to  the  (stronger)  2000C348  event,  although  it  shows  slightly  lower  power  at  lower  fre¬ 
quencies. 

•  Between  8  and  16  Hz  we  find  a  factor  6  amplitude  difference  between  the  2.5  ton  and  the 
weakest  10  ton  explosion,  both  detonated  in  the  1000  m'^  cavity.  Below  8  Hz,  the  2.5  ton 
explosion  had  too  low  SNR  to  enable  a  comparison. 
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•  The  two  5  ton  charges  detonated  in  the  300  cavity  had  very  similar  spectra.  Below  5  Ife, 
the  5  ton  charge  detonated  in  the  200  cavity  had  higher  spectral  levels  than  the  300  m 
cavity  shots.  Between  5  and  8  Hz,  the  200  cavity  shot  displays  a  reduction  in  power  spec¬ 
trum  whereas  the  300m^  shots  display  an  increase.  Above  8  Hz  the  spectral  levels  were  all 
similar. 

•  When  comparing  the  spectral  levels  of  the  “smallest”  10  ton  explosion  detonated  in  the  1000 
m^  cavity  with  the  5  ton  explosions  detonated  in  the  300  m^  cavity,  we  find  the  amphtudes 
of  the  5  ton  explosion  exceeds  the  amplitudes  of  the  10  ton  explosion  by  a  factor  3-3.5  in  the 
band  10-16  Hz.  At  5  Hz,  the  amphtude  difference  is  increased  to  a  factor  7. 

The  spectral  characteristics  found  for  the  decoupled  explosions  as  shown  in  Fig.  6.4.7  may  be 
influenced  by  the  open  tunnel  during  the  explosions.  This  entails  that  it  is  presently  unknown  if 
these  characteristics  are  reproducible  in  a  sealed  decoupled  explosion. 

These  observations  of  the  cavity  explosions  provide  information  that  can  be  used  to  assess  the 
full  and  partial  decoupling  effects.  This  is  a  topic  that  we  plan  to  look  into  in  the  continuation 
of  this  project.  For  the  newer  events  we  have  also  collected  high-frequency  data  that  we  will 
investigate  in  more  detail.  The  high-frequency  data  was  obtained  using  portable  stations  which 
were  placed  both  alongside  fixed  instruments  within  the  arrays  and  at  selected  sites  closer  to 
the  source  locations. 

We  also  plan  to  assemble  NOA,  HFS  and  NRS  recordings  of  earthquakes  located  at  distances 
similar  to  Alvdalen  and  Mossibranden,  but  at  different  azimuths,  in  order  to  make  comparative 
studies  of  their  spectral  characteristics. 

Comparisons  between  records  from  decoupled  explosions  and  ripple  fired  mining  explosions 
are  recognized  as  an  important  future  target. 
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6.5  Site-Specific  Generalized  Beamforming  (SSGBF)  applied  to  the  Lop 
Nor  test  site 

Introduction 

In  the  preceding  NORSAR  Semiannual  Technical  Summary  we  reported  on  our  initial  efforts 
to  develop  an  optimized  site-specific  threshold  monitoring  (SSTM)  system  for  the  Lop  Nor  test 
site  in  China  (Lindholm  et.  al.,  2002).  Our  aim  was  to  calibrate  all  available  high-quality  seis¬ 
mic  stations  within  regional  distance,  and  supplement  with  the  best  IMS  arrays  at  teleseismic 
distances.  As  a  calibration  data  base,  we  used  several  past  nuclear  explosions  at  the  test  site, 
combined  with  nearby,  well-recorded  earthquakes. 

We  have  now  developed  a  site-specific  Generalized  Beamforming  (SSGBF)  approach  to  sup¬ 
plement  the  SSTM  system  for  monitoring  the  Lop  Nor  test  site.  This  paper  describes  the  meth¬ 
odology  and  initial  results  using  the  SSGBF  approach.  We  also  comment  briefly  on  the  relative 
merits  of  the  two  approaches  (SSTM  and  SSGBF),  and  the  extent  to  which  they  complement 
each  other  in  this  particular  case  study. 

Development  of  site-specific  GBF 

The  Generalized  Beamforming  (GBF)  technique,  originally  developed  by  Ringdal  and  Kv^ma 
(1989),  is  now  widely  accepted  as  the  most  efficient  method  for  associating  seismic  phases 
from  a  global  or  regional  network.  In  a  typical  implementation,  a  large  number  of  generalized 
beams  are  steered  to  the  points  in  a  global  or  regional  grid.  An  automatic  detector  is  applied 
to  each  station  or  array  in  the  network,  and  a  set  of  “box-car”  functions  is  generated  for  each 
station,  so  that  the  non-zero  parts  of  these  functions  correspond  to  a  time  interval  around  a 
detection.  By  adding  these  functions  with  appropriate  weights  and  with  time  delays  corre¬ 
sponding  to  the  particular  phase-station-grid  point  combination,  one  obtains  a  “beam”  that  may 
then  be  subjected  to  a  detector  algorithm. 

A  common  feature  of  such  implementations  is  the  need  to  group  the  large  number  of  individual 
beam  detections”  so  as  to  eliminate  side-lobes  and  associate  the  detection  group  with  an  event 
at  the  correct  hypocenter.  This  grouping,  which  can  become  quite  complex,  is  clearly  unavoid¬ 
able  when  a  large  geographical  region  is  to  be  monitored.  However,  if  one  is  monitoring  a  par¬ 
ticular  site,  such  as  a  suspected  nuclear  test  site,  it  is  possible  to  simplify  this  procedure 
considerably,  and  at  the  same  time  optimize  the  parameter  settings  to  ensure  the  best  possible 
detection  probability  for  the  target  site.  This  idea  was  first  tested  by  Ringdal  and  Kv£ema 
(1993)  to  monitor  the  aftershocks  of  a  large  earthquake  sequence  occurring  in  Western  Cauca¬ 
sus  during  the  GSETT-2  experiment.  They  concluded  that  the  approach  showed  a  superior  per¬ 
formance  compared  to  the  association  procedures  being  employed  at  the  four  experimental 
international  data  centers  operating  during  GSETT-2.  In  the  present  paper  we  elaborate  further 
on  this  site-specific  GBF  (SSGBF)  approach  to  monitor  the  Lop  Nor  test  site. 

Data  base  and  analysis  procedure 

The  data  base  for  this  study  has  consisted  of  recordings  by  an  experimental  seismic  network 
composed  of  stations  and  arrays  at  teleseismic  as  well  as  regional  distances  from  the  Lop  Nor 
test  site.  We  develop  and  test  a  parameter  set  calibrated  towards  the  test  site,  and  discuss  the 
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SSGBF  results  in  conjunction  with  the  results  from  the  site-specific  threshold  monitoring 
(SSTM)  technique  applied  to  the  same  site. 

For  the  purposes  of  this  study,  18  seismic  stations  (including  9  arrays)  were  initially  chosen  as 
shown  in  Fig.  6.5.1.  These  stations  are  a  subset  of  the  23  stations  used  in  the  study  by  Lindholm 
et.  al.  (2002).  Table  6.5.1  summarizes  the  station  information  and  also  lists  the  range  of  the  azi¬ 
muth  and  slowness  parameters  used  for  forming  the  heam  towards  the  test  site.  Only  P-phases 
are  used  in  this  study.  We  note  that  the  slowness  restriction  is  only  defined  for  the  array  sta¬ 
tions,  whereas  for  the  3-component  stations  we  do  not  impose  any  restriction  on  the  slowness 
value.  The  azimuth  restrictions  are  defined  for  all  stations,  and  have  been  set  to  +!-  30  degrees 
of  true  azimuth. 

The  limits  on  arrival  time  tolerances  are  not  defined  in  Table  6.5.1,  since  these  tolerances  are  to 
a  large  extent  dependent  on  the  desired  sharpness  of  the  generalized  beam.  We  have  used  the 
formula: 

dT=S*R/111.13 

where  dT  is  the  time  tolerance  in  seconds,  R  is  the  desired  beam  radius  in  km  and  S  is  the  P- 
phase  slowness  (in  s/km)  corresponding  to  the  distance  from  the  target  site  to  the  station  being 
processed.  If  a  beam  radius  of  e.g.  50  km  is  desired,  and  the  expected  P-wave  phase  velocity  at 
a  given  station  is  8  km/s,  then  the  assumed  range  of  the  time  tolerance  is  approximately  +!-  6 
seconds. 

In  general,  there  is  a  trade-off  between  the  sharpness  of  the  generalized  beam  and  the  size  of 
the  tolerance  intervals  of  the  parameters.  We  have  chosen  to  use  fairly  large  tolerances,  to 
ensure  that  there  is  only  a  small  probability  of  missing  real  detections  corresponding  to  an 
event  at  the  site.  In  future  studies,  these  parameters  may  be  fine  tuned  to  take  advantage  of 
those  arrays  and  3-component  stations  that  are  particularly  well  calibrated.  For  example,  for 
well-calibrated  stations  the  azimuth  tolerances  may  be  greatly  narrowed,  resulting  in  a  reduced 
false  association  rate  at  no  loss  in  detectability  for  real  on-site  events. 

The  beamforming  procedure  follows  the  GBF  standard,  except  that  only  one  generalized  beam 
is  formed  in  the  site-specific  case.  The  main  steps  are: 

•  Applying  an  automatic  detector  at  each  of  the  stations/arrays  in  the  network 

•  Adding  together  weighted  “box-car  functions”  representing  the  detector  outputs  with  the 
appropriate  restrictions  on  travel  time,  azimuth  and  slowness 

•  Applying  a  thresholding  procedure  on  the  resulting  generalized  beam 

We  have  used  0-1  weights  for  the  beamforming  in  our  initial  analysis,  with  standard  (rectangu¬ 
lar)  box-car  functions.  In  addition,  we  have  made  some  experiment  with  using  alternative  func¬ 
tions  (e.g.  triangular  functions  centered  at  the  expected  arrival  time).  We  will  return  to  the 
results  of  these  experiments  in  a  later  paper.  For  now,  we  proceed  to  present  preliminary  results 
using  the  standard  approach. 
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Table  6.5.1.  The  stations  used  for  Site-Specific  Generalized  Beamforming  processing  of  the  Lop  Nor  test 
site.  In  order  for  a  signal  detection  to  be  used  as  a  candidates  an  event  at  the  Lop  Nor  test  site,  its  azimuth 
and  slowness  ^timates  has  to  be  within  the  ranges  given  in  the  table.  Notice  that  for  three-component 
stations  only  the  azimuth  range  is  used. 


Station 

Type 

Delta 

P 

travel 

time 

Azimuth 

Range 

(deg) 

Slowness 

Range 

(sec/deg) 

ARCES 

Array 

42.47 

478.2 

60-124 

5.5  -  12.5 

ASAR 

Array 

77.16 

719.2 

300  -  359 

3.0  -  9.0 

CMAR 

Array 

24.64 

326.8 

305  -  359 

6.3  -  12.3 

FINES 

Array 

41.75 

472.4 

58-118 

5.2-11.2 

GERES 

Array 

51.41 

548.6 

42  - 102 

6.3  -  12.3 

HFS 

Array 

47.93 

521.0 

45  - 105 

3.2  -  9.2 

ILAR 

Array 

65.57 

645.9 

283  -  343 

2.1  -8.1 

NORES 

Array 

48.84 

528.5 

50-110 

4.9  - 10.9 

MKAR 

SC/Array 

6.88 

103.3 

120  - 170 

10.1  -  16.1 

ARU 

3C 

24.43 

322.8 

r  84-144 

NA 

BRVK 

3C 

16.87 

239.6 

95  - 155 

NA 

KURK 

3C 

11.47 

165.7 

109-169 

NA 

KZA 

3C 

9.92 

149.8 

58-118 

NA 

NIL 

3C 

14.38 

209.0 

22-82 

NA 

TKM2 

3C 

9.68 

146.0 

63  - 123 

NA 

ULHL 

3C 

9.18 

139.8 

60-120 

NA 

ULN 

3C 

14.64 

216.4 

221  -281 

NA 

USP 

3C 

10.51 

156.6 

65  - 125 

NA 

Preliminary  results 

We  present  below  some  preliminary  results  in  applying  site-specific  generalized  beamforming 
to  the  Lop  Nor  test  site  in  China.Two  data  sets  for  performance  testing  were  compiled.  The  first 
data  set  covered  one  day  (10  September  2001),  and  in  these  data  the  recordings  of  4  explosions 
were  scaled  and  embedded  at  some  (but  not  all  of)  the  stations  (see  Lindholm  et.  al.,  2002). 
Secondly,  a  10  day  test  period  with  data  from  August  2  through  1 1, 2001  was  selected  and  data 
were  collected  for  the  stations. 

Example  1: 10  August  2001 

Our  first  example  is  the  day  10  August  2001.  Figs.  6.5.2  and  6.5.3  show  the  SSGBF  traces  for 
arrays  and  single  stations  respectively  for  that  day.  In  Fig.  6.5.2,  the  peaks  of  the  panels  for 
each  array  represent  detections  that  fall  within  the  Lop  Nor  azimuth  and  slowness  ranges  given 
m  Table  6.5.1.  To  align  the  detections  we  have  subtracted  the  P  travel-time  from  Lop  Nor  to  the 
respective  arrays.  The  network  trace  on  top  is  calculated  by  adding  boxcar  functions  surround¬ 
ing  each  detection,  using  both  the  arrays  of  this  figure  and  the  three-component  stations  of  Fig¬ 
ure  6.5.3.  In  Fig.  6.5.3,  the  same  type  of  display  is  shown,  but  only  for  the  3-component 
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stations.  Furthermore,  in  Fig.  6.5.3,  the  top  trace  uses  only  detector  information  from  the  3- 
component  stations  in  the  network. 

Alerts  (red  diamonds)  are  declared  when  the  network  trace  (the  full  network)  has  at  least  three 
matching  detections,  where  at  least  one  detections  has  to  be  at  an  array.  The  red  diamonds  at 
the  individual  arrays  or  3-component  stations  on  the  two  figures  show  detections  that  fall 
within  the  alert  interval  defined  on  the  network  trace.  The  event  at  03:15  is  located  in  northern 
Xinjang,  about  200  km  north  of  the  Lop  Nor  test  site.  It  is  interesting  to  note  that  this  event  trig¬ 
gers  a  clear  detection  on  the  generahzed  beam  steered  towards  Lop  Nor,  even  if  the  distance  to 
the  test  site  exceeds  the  beam  coverage  radius  of  50  km.  In  practice,  it  is  difficult  to  avoid  such 
side  lobe  detections  at  nearby  sites,  and  in  any  case,  it  may  be  useful  to  have  indications  on  the 
beam  trace  that  such  nearby  events  may  have  occurred.  The  two  alerts  between  19:00  and 
21:00  are  caused  by  signals  of  unknown  origin.  Most  likely  these  alerts  are  false  alarms  caused 
by  coinciding  detections  at  the  different  stations. 

Example  2:  5  August  2001 

As  another  example.  Fig.  6.5.4  shows  (in  the  top  panels)  SSGBF  processing  results  for  5 
August  2001.  Only  the  arrays  and  the  full  network  trace  are  shown.  Two  alerts  are  declared, 
and  we  note  that  several  detections  fall  within  the  alert  interval.  However,  these  detections  do 
not  add  constructively  to  the  network  trace,  indicating  that  they  were  associated  with  an  event 
outside  the  source  region.  At  the  bottom  panel  of  Fig.  6.5.4  we  show  the  corresponding  SSTM 
results  (presented  by  Lindholm  et.  al.,  2002),  having  threshold  peaks  at  the  same  times  as  the 
SSGBF  network  trace.  These  peaks  are  verified  to  be  caused  by  events  far  from  the  Lop  Nor 
test  site,  (Andaman  Islands  and  Afghanistan,  respectively). 

Example  3:  2  August  2001 

Our  third  example  is  similar  to  Example  2,  but  covers  a  “quiet”  day  (2  August  2001)  during 
which  there  were  no  large  seismic  events  globally,  and  no  significant  event  near  Lop  Nor.  The 
top  panels  of  Fig.  6.5.5  show  the  SSGBF  processing  results.  No  alerts  are  declared.  The  corre¬ 
sponding  SSTM,  shown  below  the  red  fine,  has  no  significant  peaks  and  a  background  level 
varying  between  2.6  and  3.0. 

Example  4: 10  September  2001 

Our  final  example  (Figs.  6.5.6  and  6.5.7)  is  the  one  day  (10  September  2001)  which  contained 
the  scaled,  embedded  events,  and  the  combined  network  SSGBF  trace  (using  aU  stations)  is 
shown  on  top  of  each  figure.  In  Fig.  6.5.6,  which  gives  the  standard  display  of  the  arrays  and 
the  combined  network  SSGBF  traces,  we  note  the  large  number  of  alerts  declared  automati¬ 
cally.  A  total  of  1 1  alerts  are  declared  for  that  day  by  the  SSGBF  method.  These  alerts  corre¬ 
spond  to  the  eight  embedded  explosions  (scaled  magnitudes  3.0  to  3.5),  one  additional  trigger 
that  turned  out  to  be  a  real  event  near  Lop  Nor,  and  two  alerts  of  unknown  origin. 

Fig.  6.5.7  shows  a  combination  of  SSGBF  and  SSTM  results  for  the  same  day  (10  September 
2001)  with  scaled,  embedded  events.  The  top  panel  shows  the  SSGBF  network  trace  for  10 
September  2001.  Below  is  shown  the  SSTM  network  trace  for  the  same  day.  Red  circles  indi¬ 
cate  events  for  which  the  SSTM  processing  issued  an  alert,  whereas  blue  circles  indicate  events 
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for  which  an  alert  was  not  issued.  Notice  the  good  correspondence  between  the  SSGBF  alerts 
(red  diamonds)  and  the  SSTM  alerts. 

Conclusions 

From  these  tests  the  following  preliminary  observations  can  be  stated: 

•  Out  of  the  eight  embedded  explosions  on  10  September  2001,  six  were  flagged  by  the  auto¬ 
matic  SSTM  detector,  whereas  all  eight  were  flagged  by  the  automatic  SSGBF  detector. 

The  quiet  day  (2  August  200 1 )  did  not  have  any  events  flagged,  neither  by  the  SSTM  nor  the 
SSGBF  detector. 

•  On  5  August  2001  the  SSGBF  as  well  as  the  SSTM  method  triggered  on  two  large  teleseis- 
mic  events. 

The  above  observations  should  be  evaluated  under  the  perspective  that  the  study  is  of  a  very 
preliminary  nature,  and  that  further  optimization  and  tuning  of  the  SSGBF  method  is  clearly 
possible.  Furthermore,  the  addition  of  newly  developed  and  planned  IMS  stations  would  con¬ 
tribute  to  a  significant  improvement  in  detector  performance. 

The  false  association  problem  of  the  SSGBF  approach  has  not  yet  been  evaluated  in  detail,  but 
preliminary  indication  are  that  the  likelihood  of  spurious  associations  is  low,  even  with  the 
fairly  relaxed  set  of  parameter  tolerances  used  in  this  study.  A  more  extensive  study  of  this 
topic  is  being  planned. 

The  combination  of  the  SSTM  and  the  SSGBF  techniques  appears  to  be  a  very  promising  new 
development.  Each  of  these  two  approaches  has  its  specific  strengths.  The  SSTM  technique  has 
as  its  main  strength  the  ability  to  display  the  real  seismic  field,  without  regard  to  “station  detec¬ 
tor  performance”.  On  the  other  hand,  the  SSGBF  technique  takes  advantage  of  the  individual 
station  detector  outputs,  and  uses  this  combined  information  to  narrow  down  the  number  of 
possible  candidates  for  events  in  the  target  area. 

From  the  above  we  conclude  that  the  Site-Specific  Generalized  Beamforming  performance 
tests  for  the  Lop  Nor  test  site  were  successful.  It  is  expected  that  these  initial  tests  will  be  fol¬ 
lowed  by  more  detailed  studies,  where  in  particular  the  calibration  parameters  will  be  more 
firmly  established  and  the  combined  benefits  of  the  SSTM  and  the  SSGBF  methods  will  be  fur¬ 
ther  explored. 


Torinod  Kvaema 
Erik  Hicks 
Erode  Ringdal 
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Fig.  6.5.1  Mapsjhowing  arrays  and  stations  used  for  Site-Specific  GBF  processing  of  the  Lop 
Nor  Test  Site  (indicated  with  red  star). 
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Site  specitic  GBF  -  Lop  Nor  test  site 
10  August  2001  (day  222) 
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Fig.  6.5.2.  The  peaks  of  the  panels  for  each  array  represent  detections  that  fall  within  the  Lop  Nor 
azimuth  and  slowness  ranges  given  in  Table  6.5.1.  To  align  the  detections  we  have  sub¬ 
tracted  the  P  travel-time  from  Lop  Nor  to  the  respective  arrays.  The  network  trace  on  top  is 
calculated  by  adding  boxcar  Junctions  surrounding  each  detection,  using  both  the  arrays  of 
this  figure  and  the  three-component  stations  of  Figure  6.5.3.  Alerts  (red  diamonds)  are 
declared  when  the  network  trace  has  at  least  three  matching  detections,  where  at  least  one 
detections  has  to  be  at  an  array.  The  red  diamonds  at  the  individual  arrays  show  detections 
that  fall  within  the  alert  interval  defined  on  the  network  trace. 

The  event  at  03:15  is  located  in  northern  Xinjang,  about  200  km  north  of  the  Lop  Nor  test 
site.  The  two  alerts  between  19:00  and  21:00  are  caused  by  signals  of  unknown  origin. 
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Site  specific  GBF  -  Lop  Nor  test  site 
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Fig.  6.5.3.  The  peaks  of  the  panels  for  each  3C  station  represent  detections  that  fall  within  the  Lop 
Nor  azimuth  ranges  given  in  Table  6.5. 1 .  To  align  the  detections  we  have  subtracted  the  P 
travel-tune  from  Lop  Nor  to  the  respective  stations.  The  3C  network  trace  on  top  is  calcu¬ 
lated  by  adding  boxcar  functions  surrounding  each  detection,  using  the  three-component 
stations  only.  The  red  diamonds  at  the  3C  network  trace  and  the  individual  stations  show 
l^eak^  that  fall  within  the  alert  interval  defined  on  the  full  network  trace  shown  in  Figure 
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Large  earthquakes  far  from  Lop  Nor 


Fig.  6.5.4.  The  top  panels  show  SSGBF  processing  results  for  5  August  2001.  Two  alerts  are 

declared,  and  we  find  that  several  detections  fall  within  the  alert  interval.  However,  these 
detections  do  not  add  constructively  to  the  network  trace,  indicating  an  event  outside  the 
source  region.  Below  the  red  line  we  show  the  corresponding  SSTM  results  (presented  in  the 
previous  NORSAR  Semiannual  Report),  having  threshold  peaks  at  the  same  time  as  the 
SSGBF  network  trace.  These  peaks  are  verified  to  be  caused  by  events  far  from  the  Lop  Nor 
test  site,  respectively,  Andaman  Islands  and  Afghanistan. 
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Fig.  6.5.5.  The  top  panels  show  SSGBF  processing  results  for2  August  2001.  No  alerts  are 

aeclared.  The  conesponding  SSTM,  shown  below  the  red  liney  has  no  significant  peaks  and 
a  background  level  varying  between  2.6  and  3.0. 
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Site  specitic  GBF  -  Lop  Nor  test  site 
10  September  2001  (day  253) 
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Fig.  6.5.6.  The  top  panels  show  SSGBF  processing  results  forlO  September  2001,  with  several 
embedded  Lop  Nor  event  recordings  scaled  to  magnitudes  3.0  and  3.5.  Notice  the  large 
number  of  alerts  declared. 
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Fig.  6.5. 7.  The  top  panels  show  SSGBF  network  trace  for  10  September  200 1. 

Below  is  shown  the  Network  SSTM  trace  for  the  same  day.  Red  circles  indicate  events  for 
winch  the  SSTM  processing  issued  an  alert,  whereas  blue  circles  indicate  events  for  which 
an  alert  was  not  issued. 

Notice  the  good  corre.spondence  between  the  SSGBF  alerts  (red  diamonds)  and  the  SSTM 
alerts.  The  event  at  03: 1 5  is  located  in  northern  Xinjang,  about  200  km  north  of  the  Lop  Nor 
test  site.  ‘ 
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6.6  Regional  Seismic  Threshold  Monitoring 

Sponsored  by  Defense  Threat  Reduction  Agency,  Contract  No.  DTRAOl-OO-C-0107 


Abstract 

A  database  comprising  a  total  of  45  events,  selected  to  provide  the  best  possible  ray  path  cover¬ 
age  of  the  Barents  Sea  and  adjacent  areas,  was  compiled  and  reanalyzed  in  a  consistent  manner. 
This  resulted  in  new  regional  attenuation  relations  for  Pn  and  Sn,  together  with  a  preferred 
average  velocity  model  to  be  used  for  predicting  the  travel  times  of  regional  phases.  We  have 
now  applied  these  attenuation  relations  to  investigate  a  regional  threshold  monitoring  scheme 
for  the  Barents  Sea  area. 

A  grid  system  with  an  approximately  100-km  grid  spacing  was  deployed  for  the  Barents  Sea 
region,  and  the  observations  at  the  arrays  ARCES,  SPITS,  FINES  and  NORES  were  then  used 
for  calculating  threshold  magnitudes  for  each  of  the  grid  points.  During  an  interval  without 
seismic  signals,  the  threshold  magnitudes  showed  large  variations  over  the  region,  and,  in  par¬ 
ticular,  in  the  vicinity  of  each  array.  However,  for  the  region  around  the  island  of  Novaya 
Zemlya  the  variations  are  modest,  varying  around  a  mean  of  magnitude  2. 1-2.2. 

In  order  to  investigate  in  more  detail  the  variations  in  threshold  magnitudes  for  the  Novaya 
Zemlya  region,  we  deployed  a  dense  grid  with  an  areal  extent  of  about  500  x  500  km  around 
the  former  Novaya  Zemlya  nuclear  test  site.  For  each  of  the  grid  nodes,  we  calculated  magni¬ 
tude  thresholds  for  the  two-hour  time  interval  00:00  -  02:00  on  23.  February  2002.  At 
01:21:12.1  there  was  an  event  with  a  magnitude  of  about  3,  located  about  100  km  northeast  of 
the  former  nuclear  test  site. 

Regions  of  different  sizes  were  constructed  by  selecting  grid  points  within  different  radii  from 
the  former  nuclear  test  site.  Average,  minimum  and  maximum  threshold  magnitudes  were  cal¬ 
culated  for  circular  regions  with  radii  of  20,  50, 100  and  200  km  respectively. 

The  most  important  result  is  that  even  for  a  target  region  with  radius  as  large  as  100  km,  the 
variations  in  threshold  magnitudes  are  all  within  0.2  magnitude  units.  This  applies  both  for  the 
time  interval  with  the  event  and  for  background  noise  conditions.  For  the  investigated  station 
geometry,  it  will  therefore  be  meaningful  to  represent  the  monitoring  threshold  of  the  entire 
Novaya  Zemlya  region  with  the  values  of  a  single  target  point,  together  with  the  d  priori  deter¬ 
mined  uncertainty  bounds. 

We  can  therefore  conclude  that  the  experimental  site-specific  threshold  monitoring  which  has 
been  run  daily  by  NORSAR  for  the  past  5  years,  aiming  at  the  Novaya  Zemlya  test  site,  can  be 
used  with  only  minor  adjustments  to  assess  the  threshold  within  100  km  of  the  site.  This  moni¬ 
toring  has  shown  that  between  16  August  1997  and  23  February  2002,  the  threshold  level  has 
been  consistently  below  2.5,  except  during  “interfering”  large  regional  or  teleseismic  events 
located  well  outside  the  target  region. 

For  areas  with  larger  variations  in  threshold  magnitudes,  like  in  the  vicinity  of  the  arrays,  a 
100-km  radius  target  region  will  obviously  show  larger  differences  between  the  maximum  and 
minimum  values.  Examples  illustrating  this  point  will  be  shown. 
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Objective 

The  main  objective  of  this  research  is  to  develop  and  test  a  new,  advanced  method  for  applying 
regional  seismic  array  technology  to  the  field  of  nuclear  test  ban  monitoring.  To  that  end,  we 
address  the  development  and  testing  of  a  method  for  optimized  seismic  monitoring  of  an 
extended  geographical  region,  using  a  sparse  network  of  regional  arrays  and  three-component 
stations.  Our  earlier  work  on  optimized  site-specific  threshold  monitoring  serves  as  a  basis  for 
the  development  of  this  new  method.  Emphasis  of  the  research  is  on  algorithms  that  can  be  effi¬ 
ciently  applied  in  a  real-time  monitoring  environment,  which  are  using  primarily  automated 
processing,  and  which  can  be  readily  implemented  in  an  operational  CTBT  monitoring  system. 

Research  Accomplished 

Travel-times  and  attenuation  relations  for  regional  phases  in  the  Barents  Sea  region 

A  database  containing  45  events  in  the  Barents  Sea  region  (Fig.  6.6.1)  has  been  compiled  and 
analyzed  with  the  aim  of  evaluating  crustal  models,  travel-times  and  attenuation  relations  in  the 
context  of  performing  regional  detection  threshold  monitoring  of  this  region.  The  45  events  are 
mostly  located  around  the  circumference  of  the  study  area  due  to  the  virtually  aseismic  nature 
of  the  Barents  Sea  itself.  Regional  Pn  and  Sn  phases  were  observable  for  most  events  in  the 
database,  while  Pg  and  Lg  phases  were  only  observable  for  events  with  ray  paths  that  do  not 
cross  the  tectonic  structures  in  the  Barents  Sea.  This  corroborates  a  number  of  previous  obser¬ 
vations  of  Lg-wave  blockage  within  the  Barents  Sea. 


lg.  6.6.1.  Events  (circles)  and  seismic  stations  used  for  deriving  wave  propagations  characteris¬ 
tics  of  the  Barents  Sea  and  surrounding  areas.  Array  stations  are  shown  as  squares,  while  3C 
stations  are  shown  as  triangles.  The  symbol  size  for  the  events  are  proportional  to  the  net¬ 
work  magnitudes. 
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In  order  to  estimate  magnitudes,  short  term  average  (STA)  and  spectral  amplitude  values  were 
calculated  in  several  frequency  bands  for  all  phase  arrivals  in  the  data  base.  There  were  no  sig¬ 
nificant  differences  between  spectral  and  STA  amplitudes,  so  the  latter  were  used  as  this 
parameter  is  more  efficient  to  calculate  in  real-time  processing.  A  joint  inversion  of  the  regional 
phases  Pn,  Pg,  Sn  and  Lg  was  performed  in  order  to  determine  attenuation  relations  specific  for 
this  region  according  to  the  relation  given  in  equation  1 . 

=  \ogA  -  loge  ■  •  /  +  ia.f  +  +  ^ik  + 

where  A  is  observed  STA  amplitude  within  a  frequency  band  with  logarithmic  center  frequency 
/,  A  is  epicentral  distance,  the  a^j  term  represents  the  total  attenuation  for  phase  i  out  to  the  ref¬ 
erence  distance  of  200  km,  and  a,-  and  bi  are  phase-dependent  attenuation  constants.  The  inver¬ 
sion  results  for  the  parameters  a  ,  a  and  b  are  given  in  Table  6.6.1. 

Table  6.6.1.  The  inversion  results  for  the  a,  b  and  a”  coefficients  (Is)  for  Pn,  Sn,  Pg  and  Lg  phases 
used  in  the  attenuation  relation  (equation  1). 


Phase 

a 

b 

- ^ 

Pn 

-0.002  +  0.023 

2.340  ±  0.099 

0.584  ±  0.030 

s„ 

0.141  ±0.028 

2.021  ±0.110 

0.419  ±0.037 

Pg 

0.091  ±0.084 

0.851  ±0.366 

-0.538  ±  0.035 

0.534  ±  0.062 

-0.1 86  ±0.123 

-0.609  ±  0.063 

Fig.  6.6.2  shows  a  comparison  between  the  phase  magnitude  residuals  calculated  using  the 
relations  and  parameters  of  Jenkins  et  al.  (1998)  and  our  inversion  results.  Our  results  using  the 
relations  and  parameters  of  Jenkins  et  al.  (1998)  revealed  a  relatively  high  scatter  between  indi¬ 
vidual  station  and  phase  magnitudes,  and  also  some  systematic  inconsistencies,  most  notably 
magnitudes  calculated  from  different  frequency  bands  at  the  same  station.  Magnitudes  calcu¬ 
lated  from  STA  values  in  the  2-4  Hz  passband  are  mostly  higher  than  magnitudes  calculated  in 
the  3-6  Hz  passband,  which  again  are  generally  higher  than  magnitudes  calculated  from  the  4-8 
Hz  passband.  The  coefficients  used  in  this  case  were  determined  using  data  from  eastern  North 
America,  central  Asia  and  Australia.  However,  this  relation  is  not  primarily  intended  for  local 
magnitude  calculation,  and  some  of  the  scatter  in  the  magnitudes  from  Pg  and  Lg  arrivals  in 
particular  may  be  due  to  the  small  distance  for  some  of  these  observations,  below  the  lower  dis¬ 
tance  limit  of  1.8°  used  by  Jenkins  et  al.  (1998). 

Phase  magnitudes  calculated  using  equation  1  and  the  parameters  of  Table  6.6.1  are  shown  to 
the  right  of  Fig.  6.6.2.  These  results  show  that  the  scatter  (expressed  as  standard  deviation)  was 
significantly  reduced  compared  to  the  original  calculations.  There  is  also  no  apparent  fre¬ 
quency  dependency  in  the  magnitude  residuals. 
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Fig.  6.6.2.  Magnitudes  calculated  using  the  relation  from  Jenkins  et  al.  (1998)  (left)  and  this  study 
(right)  from  individual  amplitude  readings,  plotted  V5.  network  magnitudes  for  the  45  event's 
Note  the  significant  reduction  in  scatter  (St.  dev.)  and  also  the  absence  of  frequency-depen- 
dent  effects  when  the  relation  from  this  study  is  used 


As  an  example  we  show  in  Table  6.6.2  the  individual  phase  magnitude  estimates  of  the  23  Feb¬ 
ruary  2002  event  located  on  the  northeastern  coast  of  Novaya  Zemlya.  The  consistency  of  these 
phase  magnitudes  is  remarkably  high. 
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Table  6.6.2.  Phase  magnitudes  and  network  magnitudes  for  the  23  February  2002  event  located  on  the 
northeastern  coast  of  Novaya  Zemlya,  using  the  attenuation  relations  developed  in  this  study. 


Magnitude  Type 

Magnitude 

Network  Pn 

3.19 

Network  Sn 

3.11 

Network  All 

3.15 

Station 

Phase 

Distance 

Magnitude 

AMD 

Pn 

509 

3.19 

AMD 

Sn 

509 

3.15 

LVZ 

Pn 

1055 

3.22 

LVZ 

Sn 

1055 

3.01 

SPITS 

Pn 

1095 

3.44 

SPITS 

Sn 

1095 

3.11 

ARCES 

Pn 

1144 

2.97 

ARCES 

Sn 

1144 

3.08 

KBS 

Pn 

1197 

3.16 

KBS 

Sn 

1197 

3.19 

FINES 

Pn 

1850 

3.17 

Fig.  6.6.3  shows  comparisons  of  corrected  network  magnitudes  compared  to  magnitudes  calcu¬ 
lated  from  individual  phases.  Although  relative  P  and  S  magnitudes  could  be  used  as  an  aid  in 
discriminating  between  earthquakes  and  explosions,  Fig.  6.6.3  shows  that  regional  path  effects 
in  this  area  also  give  rise  to  substantial  differences  in  magnitude.  This  is  particularly  visible  for 
Pn  and  Sn  magnitudes,  as  they  are  available  for  events  covering  the  entire  region.  Events  that 
predominantly  have  ray  paths  within  Fennoscandia  have  larger  Sn  magnitudes,  while  the  oppo¬ 
site  is  true  for  events  that  have  ray  paths  crossing  the  sediment  basins  of  the  Barents  Sea 
(Novaya  Zemlya/Kara  Sea  and  the  western  Barents  Sea/Mid-Atlantic  ridge  areas). 

From  this  study,  it  is  clear  that  Pn  and  Sn  are  the  most  useful  phases  for  calculating  magnitudes 
for  events  in  the  Barents  Sea.  In  fact,  Fig.  6.6.3  shows  that  Pg  and  Lg  are  mainly  observed  at 
close  epicentral  distances  (<  300  km).  This  situation  is  quite  different  from  what  we  have  previ¬ 
ously  found  for  the  Seandinavian  Peninsula  and  the  Baltic  Shield,  where  Lg  is  the  dominant 
phase  on  the  seismogram  out  to  at  least  1000  km.  Thus,  even  for  a  stable  continental  region, 
one  may  expect  quite  significant  regional  variations  in  the  magnitude  correction  factors. 


107 


NORSAR  Sci.  Rep.  2-2002 


August  2002 


-0.6  -0.4  -0.2  0.0  0.2  0.4  0.6 

Magnitude  diff. 


Fig.  6.6.3.  Nenvork  event  magnitude  comparisons  and  maps  of  the  geographical  distribution  of  the 
magnitude  differences  for  individual  phase  magnitudes  compared  to  network  magnitudes. 
Note  that  Sn  magnitudes  are  overestimated  for  events  that  have  paths  predominantly  within 
the  Baltic  shield,  while  events  with  paths  that  cross  the  Barents  Sea  have  lower  Sn  magni¬ 
tudes  Pg  and  Lg  magnitudes  appear  to  he  quite  stable  within  the  limited  distance  range  from 
which  readings  are  available. 
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Regional  Threshold  Magnitudes 

Using  the  developed  attenuation  relations  described  above,  observations  at  the  arrays  ARCES, 
SPITS,  FINES  and  NORES  were  used  for  calculating  threshold  magnitudes  for  a  grid  system 
covering  the  entire  Barents  Sea  region.  The  grid  spacing  was  approximately  100  km.  Fig.  6.6.4 
shows  the  threshold  magnitudes  during  a  time  instant  without  seismic  signals.  We  find  large 
variations  over  the  region,  and  in  particular  when  approaching  each  of  the  arrays.  However,  for 
the  region  around  the  island  of  Novaya  Zemlya,  the  variations  are  modest,  varying  around  a 
mean  of  magnitude  2.2. 


Fig.  6.6.4.  Threshold  magnitudes  for  the  time  instant  2002-054:01.11.20.0.  Notice  the  improved 

monitoring  capability  in  the  vicinity  of  each  station.  For  distances  above  1.5  degrees  of  each 
station  we  have  considered  the  Pn  and  Sn  phases,  whereas  P  g  and  Lg  have  been  used  for  dis¬ 
tances  less  than  1.5  degrees. 


In  order  to  investigate  in  more  detail  the  variations  in  threshold  magnitudes  for  the  Novaya 
Zemlya  region,  we  deployed  a  dense  grid  with  an  areal  extent  of  about  500  x  500  km  as  shown 
in  Fig.  6.6.5.  For  each  of  the  grid  nodes  we  calculated  magnitude  thresholds  for  the  2  hour  time 
interval  00:00  -  02:00  on  23  February  2002.  At  01:21:12.1  there  was  an  event  with  a  magnitude 
of  about  3.2  (see  Table  2),  located  about  100  km  north-east  of  the  former  nuclear  test  site. 

Regions  of  different  sizes  were  constructed  by  selecting  grid  points  within  different  radii  from 
the  former  nuclear  test  site.  Fig.  6.6.6  (left)  shows  the  variations  in  threshold  magnitudes  for  a 
circular  region  with  a  radius  of  20  km  around  the  test  site.  The  blue  line  shows  the  average 
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threshold,  whereas  the  red  lines  represent  the  minimum  and  maximum  values.  Fig.  6.6.6  (right) 
show  similar  curves  for  a  region  with  radius  100  km. 


?  S  a  S  f:. 


Fig.  6.6.5.  Dense  grid  deployment  around  Novaya  Zemlya  (grid  .^pacing  1 1  km ).  The  red  star  shows 
the  location  of  the  former  nuclear  test  site,  whereas  the  red  diamond  shows  the  location  of 
the  event  on  23.  February  2002  with  origin  time  01:21:12. 1. 


Fig.  6.6.6.  Threshold  magnitudes  for  the  time  interval  00:00  -  02:00  on  23.  February  2002 for 

20km  (left)  and  1 00  km  ( right)  radius  target  regions  centered  around  the  former  nuclear  test 
site^  The  peak  at  about  5000  seconds  corresponds  to  signals  for  the  3.2  event  located  about 
100  km  north-ea.st  of  the  te.st  site. The  blue  line  shows  the  average  threshold,  whereas  the  red 
lines  represent  the  minimum  and  maximum  values. 
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It  is  interesting  to  notice  that  even  for  a  region  with  100  km  radius,  the  variations  in  threshold 
magnitudes  are  all  within  0.2  magnitude  units.  For  this  particular  configuration  of  the  monitor¬ 
ing  network  relative  to  the  target  area,  it  will  therefore  be  meaningful  to  represent  the  monitor¬ 
ing  threshold  of  the  entire  region  with  the  values  of  a  single  target  point,  together  with 
uncertainty  bounds  as  shown  in  Fig.  6.6.6.  For  areas  with  larger  variations  in  threshold  magni¬ 
tudes,  like  in  the  vicinity  of  the  arrays,  a  100  km  radius  target  region  will  obviously  show  larger 
differences  between  the  maximum  and  minimum  values. 

We  would  like  to  comment  on  the  threshold  magnitude  of  the  peak  corresponding  to  the  event 
located  north-east  of  the  test  site.  In  cases  where  an  event  actually  occur  in  the  target  region,  the 
magnitude  thresholds  will  often  be  biased  slightly  low.  In  Fig.  6.6.6,  we  find  a  maximum  value 
of  about  2.9,  whereas  the  event  magnitude  is  estimated  to  3.15.  In  such  cases,  a  maximum-like¬ 
lihood  magnitude  estimation  algorithm  should  be  activated.  However,  for  small  events  with  a 
size  close  to  the  threshold  magnitudes,  this  bias  will  not  be  significant. 

Regional  Threshold  Monitoring  including  the  Amderma  station 

The  KRSC  group  in  Apatity,  Russia  has  provided  us  with  about  3  days  of  continuous  data  from 
the  station  in  Amderma,  located  on  mainland  Russia,  just  south  of  the  island  of  Novaya 
Zemlya.  The  data  interval  is  centered  around  the  origin  time  of  the  23  February  2002  event 
located  on  the  northeastern  coast  of  Novaya  Zemlya.  In  addition,  data  from  the  array  in  Apatity 
on  the  Kola  Peninsula  were  included. 

Fig.  6.6.7  shows  the  threshold  magnitudes  during  an  instant  without  seismic  signals,  using  the 
developed  attenuation  relations  for  Pn,  Pg,  Sn  and  Lg.  We  find  large  variations  over  the  region, 
and  in  particular  when  approaching  each  of  the  arrays.  With  the  Amderma  station  included,  we 
also  find  significant  variations  for  the  island  of  Novaya  Zemlya,  ranging  from  1.4  at  the  south¬ 
ern  tip  to  2.2  at  the  northern  tip.  This  implies  that  a  regional  threshold  monitoring  scheme  for 
the  NZ  region  has  to  be  divided  into  geographical  sub-regions  having  similar  threshold  magni¬ 
tudes  during  background  noise  conditions. 
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Fig.  6.6.7.  Threshold  magnitudes  for  the  time  instant  2002-054:01. 1 1.20.0,  including  data  from  the 
Amderma  station  and  the  Apatity  array.  Notice  the  improved  monitoring  capability  in  the 
vicinity  of  the  Amderma  station  as  compared  to  the  thresholds  shown  in  Fig.  6.6.4. 


Conclusions  and  Recommendations 

The  pattern  of  Lg  arrivals  and  associated  amplitudes  supports  the  previously  published  indica¬ 
tions  that  the  deep  sediment  basins  and  Moho  topography  under  the  Barents  Sea  efficiently 
block  Lg  wave  energy  from  crossing.  From  this,  it  is  clear  that  Pn  and  Sn  are  the  most  useful 
phases  for  calculating  stable  and  consistent  magnitudes  for  events  in  the  Barents  Sea. 

The  ‘BAREY’  model  from  Schweitzer  and  Kennett  (2002),  based  on  a  model  for  the  Barents 
Sea  area  from  Kremenetskaya  et  al.  (2001),  provides  the  smallest  overall  travel-time  residuals 
when  locating  events  within  the  vicinity  of  the  Barents-  and  Kara  Seas. 

The  attenuafion  in  the  Barents  Sea  region  differs  somewhat  from  that  observed  in  other  stable 
tectonic  regions,  as  evidenced  by  the  fact  that  the  coefficients  given  by  Jenkins  et  al.  (1998)  for 
such  regions  do  not  give  consistent  magnitudes  across  frequencies,  phases  and  stations  for  our 
amplitude  observations  from  the  events  in  the  Barents  Sea  region. 

Amplitude  inversion  has  been  used  in  this  study  to  resolve  new  attenuation  coefficients  and  sta¬ 
tion  corrections  for  estimating  magnitudes  from  STA  amplitude  observations  for  Pn,  Pg,  Sn 
and  Lg  phases  in  the  Barents  Sea  region.  The  distance  range  of  observations  on  which  the  Pg 
and  Lg  relations  are  based  is  limited;  a  future  study  using  a  greater  number  of  continental 
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events  could  most  likely  provide  a  relation  for  STA  based  Lg  magnitudes  that  is  applicable  at 
larger  distances,  albeit  limited  to  paths  within  Fennoscandia. 

The  seismic  station  in  Amderma  can  be  incorporated  into  the  regional  network  in  Fennoscandia 
and  on  the  Svalbard  archipelago  using  an  appropriate  crustal  model,  and  is  able  to  provide 
important  information  regarding  the  location  of  events  in  the  eastern  parts  of  the  Barents  Sea 
and  the  Kara  Sea  (Schweitzer  and  Kennett,  2002).  Magnitudes  calculated  at  this  station  are  on 
the  whole  consistent  with  the  other  observations. 

For  the  time  interval  under  study,  the  seismic  arrays  ARCES,  SPITS,  FINES  and  NORES  pro¬ 
vide  an  average  monitoring  capability  of  about  magnitude  2.2  for  the  island  of  Novaya  Zemlya. 
For  a  region  with  a  100  km  radius  around  the  former  nuclear  test  site  the  variations  in  threshold 
magnitudes  are  aU  within  0.2  magnitude  units.  It  will  therefore  be  meaningful  to  represent  the 
monitoring  threshold  of  the  entire  region  with  the  values  of  a  single  target  point,  together  with 
uncertainty  bounds  as  shown  in  Fig.  6.6.6. 

In  cases  when  data  from  the  Amderma  station  can  be  retrieved,  we  find  significant  variations  in 
threshold  magnitudes  over  the  island  of  Novaya  Zemlya,  ranging  from  1.4  at  the  southern  tip  to 
2.2  at  the  northern  tip.  For  the  actual  time  interval,  the  monitoring  capability  for  the  former 
nuclear  test  site  is  lowered  by  about  0.3  magnitude  units  to  about  1.9.  This  implies  that  a 
regional  threshold  monitoring  scheme  for  the  NZ  region  has  to  be  divided  into  geographical 
sub-regions  having  similar  threshold  magnitudes  during  background  noise  conditions. 

T.  Kvaema 

E.  Hicks 

J.  Schweitzer 

F.  Ringdal 
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